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ABSTRACT 


This thesis develops a Lanchester-type force-on-force combat model 
simulating small-unit amphibious operations. The model commences with 
a ship-to-shore assault of aggressor forces mounted onboard Landing 
Vehicle Assault craft moving against a defensive force ashore. Once 
the ship-to-shore phase of combat is completed, the mcdel continues 
to simulate land combat further inland between the assaulting aggressor 
forces and other defensive forces occupying key terrain. 

The main thrust of the thesis is to alleviate some of the problems 
associated with the inherent abstractness of Lanchester-type combat 
models; specifically, to develop "user-friendly" input-data and output 
structure, and more thorough documentation of the model's algorithms 
to provide a model which would be more easily understcod and utilized 
by students of combat modeling. 
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I. INTRODUCTION 


A. OVERVIEW 

This thesis develops a Lanchester-type force-on-force combat model 
simulating small-unit amphibious operations. The model commences with 
a ship-to-shore assault of aggressor forces (e.g,, a U.S. Marine Infantry 
Battalion), mounted onboard Landing Vehicle Assault craft (LVA) moving 
against a defensive force ashore located in fixed positions along the 
coast the aggressor force is attempting to occupy. Once the ship-to- 
shore phase of combat is completed, the model continues to simulate 
land combat further inland between the assaulting aggressor forces and 
other defensive forces occupying key terrain. 

The main thrust of the thesis is to alleviate some of the problems 
associated with the inherent abstractness of Lanchester-type combat 
models (see [Ref. 1]), and specifically to integrate and enhance work 
done in previous models, to develop "user-friendly" enhancements, and 
more thorough documentation of algorithms to provide a model which would 
be more easily understood and utilized by students of combat modeling. 

B. BACKGROUND AND GENERAL MODEL 
1. Overview 

The small-unit'amphibious operation combat model presented in 
this thesis is the result of the integration and enhancement of two 
independent combat models. The first model is a ship-to-shore combat 
model which models a ship-to-shore assault conducted by landing vehicle 
assault craft against fixed enemy positions ashore. The second model 
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is a land combat model which models a land assault conducted by LVA 
forces on a beach, against fixed enemy positions further inland, 

2. Original $hip-to-Shore Combat Model 

The ship-to-shore combat model used as a basis for this thesis 
was presented in a thesis by David L. Chadwick [Ref. 2], It modeled 
the amphibious assault of five waves of LVA against a defensive force 
composed of tanks and antitank guided missiles (ATGM) in fixed positions 
ashore. Attrition was modeled using Lanchester area-fire and aimed-fire 
equations. The purpose of developing such a model was to determine tne 
optimal design characteristics of LVA in an amphibious assault for a 
given combat scenario. The optimal design of an LVA was considered 
to be that design which produced the lowest level of LVA attrition for 
the given combat scenario. 

3. Original Land Combat Model 

The land combat model used as a basis for this thesis was 
developed in Joseph Smoler's thesis [Ref. 3]. It modeled land combat 
conducted by three aggressor force units utilizing tanks assaulting 
three defensive force units armed with Tube-Launched, Optical-Guiaed, 
Wire-Controlled missiles (TOW's) in fixed positions. The location of 
the land combat was the Fulda Gap region in West Germany. Attrition 
was modeled using Lanchester aimed-fire equations. The purpose of 
Smoler's thesis was to develop a basic small-unit land combat model for 
determining optimal defensive unit locations for a given combat scenario. 
The optimal locations of the defensive units were considered to be those 
locations which provided the lowest level of attrition of the defensive 
units for the given combat scenario. 
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4. The Enhanced Land Combat Model 


An enhanced version of Smoler's land combat model was developed 
by Glenn Mills in his thesis [Ref. 4]. The enhancements developed by 
Mills added flexibility to Smoler's land combat model by providing 
user selected options which could be employed depending upon the abili¬ 
ties and desires of the model's user. The enhancements included the 
option of altering the aggressor force's attack routes enabling the user 
to study not only the optimal defensive unit locations, but the optimal 
aggressor force attack routes for the given defensive unit locations as 
wel 1. 

A second enhancement was the option of selecting a stochastic 
attrition-rate coefficient. This introduced the element of randomness 
into the model's attrition algorithm providing a more realistic approach 
to modeling a unit's fighting effectiveness. 

The third enhancement is the option of providing alternate 
defensive positions so that the defensive units could move to more 
defensible terrain once their original positions had become untenable. 

5. The Original Small-Unit Amphibious Warfare Model 

The original small-unit amphibious warfare model used as a 
basis for this thesis was developed by Soon Dae Park in his thesis 
[Ref. 5]. Park's model attempted to conceptualize the flow of events 
of an amphibious assault by first running the ship-to-shore model, 
followed immediately by running the land combat model. The analysis 
of this model as a class project served as the catalyst for the develop¬ 
ment of the small-unit amphibious operation combat model presented in 
this thesis. 
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6. The Analysis of Park's Model 


The class project conducted by Clay Grubb, Robert Larson, and 
this author had as its purpose the analysis of Soon Dae Park's small- 
unit amphibious operation combat model. The results of the analysis 
revealed the value of Park's thesis in providing a general scheme of 
events for the modeling of small-unit amphibious operations. The results 
also identified some enhancements that could be applied to his conceptua¬ 
lized model that would integrate the ship-to-shore and land combat models 
into a singular small-unit amphibious operation combat model. The develop¬ 
ment and application of these enhancements to Park's model served as the 
foundation for this thesis, and the development of the model presented. 

C. MAJOR GOAL AND OBJECTIVES OF THE THESIS 

1. Major Goal of the Thesis 

The overall goal of the thesis is the development of a small-unit 
amphibious operation combat model. It will be based on the integration 
and enhancement of the two combat models discussed in the previous section 
of this chapter. There are three underlying objectives of the thesis 
which will guide the development of the model toward the accomplishment 
of this goal. 

2. Objectives of the Thesis 

a. Integration of Independent Combat Models 

The first objective in the development of the model was to 
integrate two initially independent combat models into a singular con¬ 
tinuous flow combat model. This was accomplished by first allowing 
force levels at the completion of the ship-to-shore phase of combat to 
be used as the initial force levels in the land phase of combat. 
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Secondly, it was recognized that four combat modelers con¬ 
tributed to the resulting model presented by Park in his thesis. As 
such, four individualized FORTRAN coding techniques were reformulated 
into one style to provide a more tractible small-unit amphibious oper¬ 
ation combat model, 

b. User-Friendly Input-Data and Output Structure 

The second objective of the thesis was to provide a user- 
friendly combat model. It is a major contention of this thesis that 
combat modelers have not adhered closely to this principle when pro¬ 
viding combat models for the United States military. Furthermore, it 
is believed that the lack of concern given to this approach of combat 
modeling is a major reason for the less than unanimous reception that 
combat models have received by the United States military as tools for 
training its commanders aid staffs. Therefore, the model presented in 
this thesis was designed and documented with the user's capabilities 
and needs in mind as opposed to those of the programmer. 

c. Student-Oriented Combat Model 

The third objective of the thesis was to provide the student 
of combat modeling with a combat model which was easily understood and 
studied. As a result, the model presented in this thesis was designed 
with a low level of complexity to allow the student with little or no 
experience in combat modeling to understand more easily the combat 
modeling theory and its application. 
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II. NK)DEI. ENHANCEMENTS 


A. OVERVIEW 

This thesis had as its goal the development of a small-unit amphibious 
operation combat model. Guided by the three objectives discussed in 
Chapter One, five modeling enhancements were applied to the two original 
combat models serving as the foundation for the resulting small-unit 
amphibious operation combat model presented in this thesis. The enhance¬ 
ments provide for the proper integration of the ship-to-shore and land 
combat models. In addition, they have contribjted to the development 
of a more user-friendly combat model which can be used to assist combat 
modeling students in their understanding the theory of combat modeling 
and its application. 

B. INTEGRATION OF SHIP-TO-SHORE AND LAND COMBAT MODELS 

The intent of the model presented here is to view the amphibious 
assault as a continuous process made up of two phases of combat (ship- 
to-shore, and land combat) where the land combat phase is dependent 
upon the outcome of the ship-to-shore combat pnase of the model. 

Implementation of this enhancement called for the creation of a 
new variable. Total Landing Force Ashore (TLF), which would accumulate 
the surviving landing force of each assault wave as it reached the beach. 
This total landing force ashore would than be redistributed into three 
main assault units for the land combat phase of the model. The 
rationale for the redistribution of forces is based on realistic 
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military doctrine which is to maintain a well-balanced force when the 
strength and location of the enemy is unknown to the assaulting forces 
(as is assumed in the model). 

Since the manner in which defender force levels are determined by 
the ship-to-shore and land combat models appears to be quits realistic, 
the defending force level as modeled by Soon Dae Park was used as input 
to the land combat phase. In particular, if the aggressor force had been 
successful in routing the defending forces situated on the beach, 
defending forces situated further inland naturally would be impelled 
to defend the remaining terrain still in their possession. It should 
be noted that the size of these defending forces further inland is an 
option of the user which in itself can be varied for analysis of variant 
battle scenarios. 

C. AGGRESSOR FORCE ATTRITION—SHIP-TO-SHORE PHASE 

Attrition in Lanchester-type combat modeling is based upon the 
expected percentage of the original force remaining at a given point 
in time. The expected percentage of forces remaining then car be 
restated in terms of a real number to represent the expected number of 
forces remaining. This method of computing reduced force levels is 
considered to be quite appropriate when modeling land combat, and was 
implemented by Chadwick in his ship-to-shore combat model to simulate 
LVA attrition. However, use of Lanchester equations to model such 
vehicular attrition of a vehicle at sea was determined to be inappro¬ 
priate. Where it is a reasonable assumption that a disabled vehicle 
on land still can contribute something toward the final outcome of 
the battle if any of its weapons systems or onboard troops survive. 
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an LVA that is disabled at sea is of no use to the amphibious assault 
and subsequent land combat phase. The LVA will be recovered, and on¬ 
board troops brought to the landing site after the assault has caken 
place. 

Chadwick was not concerned with this distinction due to his model's 
purpose of modeling LVA attrition in terms of ship-to-shore movement 
only. Therefore, he simply utilized Lanchester equations in modeling 
LVA attrition resulting in fractionalized losses of assaulting LVA's. 
However, if a ship-to-shore combat model is to be properly integrated 
with a land combat model, only whole numbers of LVA's ashore should be 
used as input. Hence, an enhancement was made to the model. 

The approach was to find the integer value of the number of surviv¬ 
ing LVA's in each assault wave, and then sum these values resulting 
in the total landing force ashore (TLF). The fractional portion remain¬ 
ing was considered to be those LVA's disabled at sea and unable to 
participate in the land combat phase of the operation. 

D. STOCHASTIC ATTRITION-RATE COEFFICIENT MODIFICATION 

Mill's land combat model allowed the user the option of selecting 
either deterministic or stochastic attrition-rate coefficients to be 
used in assessing the attrition of opposing forces. The justification 
for utilizing stochastic attrition-rate coefficients to model force-on- 
force attrition rates was based upon the assumption that the attrition- 
rate coefficient is a random quantity measuring a unit's fighting 
ability, and can be estimated before any given battle. 

This can be illustrated by considering the expected value of a 
random variable. For example, assume a probability distribution is 
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selected for the random variable such that the expected value of the 
random variable is equal to the deterministic attrition-rate coefficient 
set for all units. When a random sample is taken from this distribution, 
the individual values assigned the random variable will serve as 
individual unit attrition-rate coefficients, where the sample mean 
will serve as the overall force attrition-rate coefficient. The result 
is that the overall force attrition-rate is equal to the sample mean, 
which is approximately equal to the population mean of the random 
variable. Recalling that this population was selected with a mean 
that equalled the deterministic attrition-rate coefficient, units now 
have their own individual attrition-rate coefficients, while the force 
attrition-rate coefficient has remained close to the intended value of 
the deterministic attrition-rate coefficient. This is more realistic 
than the deterministic option since each unit would be expected to 
have a different level of effectiveness, which necessarily would imply 
different attrition rates while maintaining one overall force level 
attrition rate. 

The attrition-rate coefficient, A^-j, is used as the measure of the 
rate a firer in Unit i attrits a target in Unit j. This has been 
likened to the fighting effectiveness of a particular Unit i. Obviously, 
this is a variable quantity influenced by a myriad of factors to include 
esprit de corps, past history of success or failure, prior exposure to 
combat, weather, quality of leadership, etc. The intent of such a 
basic model as this is to attempt to capture the overall effect of 
these factors by developing a distribution of a unit's initial fighting 
capabilities (specifically, to develop a distribution of for the 

unit). 
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Mills proposed a distribution based upon a quadradic function which 
would produce a symmetric distribution with a mean value of approximately 
0.55. This distribution restricted a unit's maximum effectiveness to 
only 80 percent of its maximum capable effectiveness level. It also 
implied that the civerage unit in combat will only perform at 55 percent 
of its maximum effectiveness level at any given time. 

A more plausible way of assigning a distribution to the A--'s might 

• J 

be a truncated Nor'mal Distribution limited to values between 0.00 and 
1.00. However, this approach would leave little flexibility in terms 
of modeling variant scenarios since the opposing forces always would 
have attrition-rate coefficients associated with that particular dis¬ 
tribution whenever the stochastic option was selected. This restriction 
is due to the programming constraints encountered in attempting to 
implement variant truncated Normal Distributions in the model. There¬ 
fore, a Beta Distribution was selected for use in the model. 

The natural range of the Beta Distribution is from 0.00 to 1.00 
thereby alleviating the burden of constructing a truncated distribution. 
Furthermore, its two scaling parameters, P and Q, can be selected readily 
and input by the user to construct virtually any variant of the Beta 
Distribution so desired. The specific values selected for P and Q would 
parameterize the distribution of the A^j's according to the user's par¬ 
ticular combat scenario without the burden of reprogramming the distri¬ 
bution on each successive run of the model. 

The density function for the Beta Distribution is as follows: 
f(x) = X^’^{l-X)^“^ for 0.04 x< 1.0 






























Therefore, a P=21 and Q=7 would yield a distribution of A..'s with a 
mean of 0.75. This says that a unit with an A., of 0.75 is operating 
at 75 percent of its potential effectiveness. Whereas, a P=7 and Q=21 
would yield a distribution of A^.^'s wi'ch a mean of 0.25, indicating that 
a unit is operating at 25 percent of its potential effectiveness. 

To illustrate the flexibility of tin's approach in determining sto¬ 
chastic attrition-rate coefficients. Figure 2-1 is provided displaying 
the distribution of A^.^'s that would be obtained when the user alters 
the parameters of the Beta Distribution. The user now can model a 
strong elite force using, for example, parameter values P=21, Q=7, or 
model a weak and poorly lead force using parameter values P=7 and Q=21, 
depending upon the particular battle scenario the user is analyzing. 

While the Beta Distribution used in this thesis is different than 
the Quadradic Distribution used by Glenn Mills, the implementation of 
this distribution for the attrition-rate coefficients is exactly the 
same as originally modeled. Since it was assumed earlier that the fight¬ 
ing effectiveness of each unit is a random quantity prior to a given 
battle, it is only necessary to obtain a realization of the random 
variable for each unit prior to the initialization of the battle. This 
realization. A?,., is determined by the user-supplied inputs P and Q, 
and subsequent calls to a Beta Distribution Random Deviate Generator 
[Ref. 6]. Therefore, an attrition-rate coefficient is computed for 
each unit using the following equation: 

' A?. X (1 - r/r for 0 ± r ^ r 

1 j e c 

^ij 

, 0 for ^ r 
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Figure 2-1. Beta Density 
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Where: A°j. = Realization of unit's fighting effectiveness 
r = Current range between firer and target 
r^ = Maximum effective range of a firer's weapon 
This function was utilized because it is a function of both range and 
A?j thus creating a different attrition-rate curve for each unit, 
depending on that unit's effectiveness level prior to the battle. A 
graphic illustration of an attrition-rate coefficient curve for an 
equal to a mean of 0.75 from the Beta Distribution where P=21 and 
Q=7, and the maximum effective range r^, of 3000 meters would look 
1 ike Figure 2-2. 



Figure 2-2. Attrition-Rate Coefficient Curve for A?. = 0.75 

and r = 3000m 
e 

To illustrate the effect that this stochastic attrition option has 
on the outcome of the model, two runs of the model were made using this 
option, while varying the Beta Distribution parameter values for both 
forces on each run. All other characteristics of both forces were left 
unaltered. In the first run, the aggressor forces were modeled to 
operate at 75 percent of their potential effectiveness, and the defending 
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forces were modeled to operate at 25 percent of their potential effec¬ 
tiveness. The battle outcome, as listed in Table 2-1, indicates that 
the aggressor forces won the battle. In the second run, the potential 
effectiveness of the opposing forces was reversed. The aggressor forces 
were now modeled to operate at 25 percent of their potential effective¬ 
ness, and the defending forces were modeled to operate at 75 percent of 
their potential effectiveness. The battle outcome, as listed in Table 
2-1, indicates that the battle was terminated due to the opposing forces 
being too close. The aggressor forces were unable to overrun the defend¬ 
ing forces, as was the case in the first run, which was due solely to the 
change in the potential fighting effectiveness of the opposing forces. 

A change in the battle outcome was expected; however, to what degree 
that change would be was unknown. The fact that the defender forces 
were unable to win the battle on the second run, while having a much 
higher level of effectiveness, indicates that other characteristics of 
the opposing forces were also playing an important role in the battle 
(e.g., types of weapons employed, original force levels, speed of attack, 
etc.). 

Through the use of the stochastic attrition option, the user now 
has the capability of studying one more facet of combat (i.e., potential 
fighting effectiveness), and can analyze to what degree different fight¬ 
ing effectiveness levels will have on final battle outcome. 

E. USER-FRIENDLY I/O STRUCTURE 

A significant and important part of writing a computer program for 
a combat model is to provide for the input and output of data to and 
from the program. It is my belief that one of the major factors 
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contributing to the lukewarm reception, in general, that combat models 
have received by the United States military is due, in part, to the 
poorly designed input-data and output structure of the combat models. 

The primary user of those models, the military commander, normally 
finds it difficult to decipher the myriad of input-data requirements, 
or the voluminous output from combat models that supposedly were de¬ 
signed for the commander's use. It is a contention of this thesis that 
if more attention was given to the development of user-friendly input- 
data/output requirements, that more interest would be generated toward 
the use of such models in training military commanders. Therefore, an 
enhancement was made to the input-data and output requirements of the 
model to demonstrate a method of alleviating this problem. 

1. User-Friendly Input Structure 

A readymade input data file was provided with the model to serve 
as a guide for entering all of the required data in the correct format 
required by the model (see Appendix C). Each variable requiring input 
for the model has been listed in the sample input file with sufficient 
space provided for ensuririg that data is entered in the correct format. 
This file, therefore, provides the unfamiliar user of the model with the 
opportunity to utilize the model with only a limited knowledge of the 
model's algorithm and input requirements. This type of user-oriented 
input requirement will alleviate some of the apprehension that an unfamil¬ 
iar user of the program might have, and might actually act as a catalyst 
in increasing the amount of use the model receives. 
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2. User-FrieiTdly Output Structure 

Indecipherable output, or too much output from a model, can be 
just as much of a deterrent to a model's use as complex input require¬ 
ments can be. This point was brought out by Ye S. Venttsal in her dis¬ 
cussion of good combat models: 

It is advisable in such "training" modeling of combat actions 
that the commander receive information from the computer not in 
the form of mean characteristics averaged over a set of realiza¬ 
tions, but rather in the form of only one specific realization, 
on the basis of which a decision is in fact made. [Ref. 7] 

To paraphrase Venttsal, the combat model output must be clear, 
concise, and identifiable to the military commander. Furthermore, it 
must answer the questions that were originally asked by the user--speci- 
fically, who won and why? 

The model output was therefore restructured to provide a concise 
listing of what input parameters were entered into the program for pro¬ 
cessing, and a concise and understandable output summary of what occurred 
throughout the battle (see Appendix F). Additional'y, a new feature was 
introduced into the model which gives the user the option of viewing 
either a detailed time-step battle summary, or just a final battle summary 
of what occurred in the running of the model. 

F. DOCUMENTATION AND PROGRAM FORMAT 

Two of the objectives of the model presented in this thesis were 
first to serve as an example of the way in which combat models should 
be designed to be user-friendly to ensure their acceptance and use in 
training military commanders; and secondly, to serve as a model for 
combat-modeling students so that they might acquire a better understand¬ 
ing of how combat models ought to be programmed into a computer. 
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It already has been discussed how the user-friendly I/O structure assists 
the user of the model. However, proper structuring of programs for read¬ 
ability and good documentation is equally necessary to ensure readability 
and understanding by students and analysts. 

The FORTRAN program presented in this thesis which integrated and 
enhanced the ship-to-shore and land combat models is an amalgamation of 
subroutines originally written by different people, with their own unique 
style of programming. The interweaving of these four styles of program¬ 
ming throughout the program seriously detracted from the smooth flew of 
program structure and readability desired when analyzing the computer 
program. Therefore, an enhancement was made to the model: the program 
was restructured so that it would follow one basic style of programming 
(see Appendix B). New labeling and structuring of formatted statements 
and nested FORTRAN functions were provided to make the computer program 
more readable. This restructuring should assist the student in under¬ 
standing the program flow, and provide an incentive to those interested 
students to develop future enhancements to the model. 

In addition to developing one style of programming, more detailed 
documentation of variable definitions and descriptions of program flow 
were added to the program. The purpose of this documentation was to 
have the program serve as a reference to itself in order that the reader 
would not be forced to refer to various manuals outside of the program 
each time an explanation of the functioning of a particular aspect of 
the program is desired. 
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III. CURRENT MODEL DESIGN 


A. OVERVIEW 

The small-unit amphibious operation combat model presented in this 
thesis consists of the integration and enhancement of a ship-to-shore 
combat submodel, and a land combat submodel. Both of the original sub¬ 
models were similar in design, basing force attrition on Lanchester- 
type expected-value equations. As presented earlier, enhancements to 
both submodels reduced the differences in design of these submodels, 
molding them into what may be called a small-unit amphibious operation 
combat model. Figure 3-1 provides the scheme for the sequence and 
general flow of events in the overall model. 

It should be noted that although the ship-to-shore and land combat 
models are quite similar, they still have their own unique characteristics 
in modeling certain events that take place throughout the battle. There¬ 
fore, in discussing the model as a whole, the two phases of the battle 
will be addressed separately, and those events which are of particular 
interest in each phase of combat will be elaborated on in order that the 
reader might acquire an overall appreciation of the contributions each 
submodel makes to the overall model. 

B. SHIP-TO-SHORE PHASE 

1. Overview 

Since the objectives of the thesis are to provide a user-friendly 
tractible combat model, a number of broad assumptions have been made 
regarding the exact method of employment of the LVA in the ship-to-shore 
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phase of the amphlDious operation. First of all, it is envisioned that 
for command and control purposes, as well as for mine clearing operations, 
there will exist LVA approach lanes as depicted in Figure 3-2, along which 
columns of LVA will transit a 25-mile distance to shore from the Amphibious 
Task Force (ATF). The 25-mile distance is based upon recent requirements 
studies indicating that in future amphibious operations, due to the 
increased lethality of anti-ship missiles and long range artillery, it 
will be necessary to increase the Amphibious Task Force standoff distance 
to approximately 25 miles from shore to reduce the vulnerability of the 
amphibious shipping against this anticipated threat [Ref. 8]. Secondly, 
it is assumed that a maneuver area will exist within which the columns of 
waves of LVA will form into a conventional landing formation composed of 
waves of landing craft as prescribed by current doctrine. 

The two previous assumptions set the stage for the primary assump¬ 
tion used in computing LVA force level attrition: that is, direct fire 
weapons will be assumeii to be the primary anti-LVA threat,--specifically, 
modified versions of current tank and antitank guided missiles (ATGM) 
assets. Although in reality some attrition of LVA can be expected in 
the maneuver area, it will be assumed that the critical exposure period 
will be that portion of time in the ship-to-shore movement that the first 
assault wave comes within 5,000 meters of the shore defenses until, up to, 
and including the arrival of the last assault wave ashore. Figure 3-3 is 
a flowchart depicting the general sequence of events of the ship-to-shore 
phase model. 

2. LVA Movement Conceptualization 

Two tactical decision variables were utilized for modeling LVA 
ship-to-shore movement: 
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Figure 3-3. Generalized Flowchart for Ship-To-Shore Phase 
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— TBW is the decision variable for the time betv/een successive 


waves. As TBW is shortened, coordination problems will arise resulting in 
confusion on the beach due to insufficient time provided for an assault 
wave to move inland prior to the next wave's arrival. The level of con¬ 
fusion generated by a short TBW must be balanced against the cost of not 
having sufficiently rapid initial buildup of offensive forces ashore. 

—RD is the distance from the shoreline that each wave will 
commence the transition from planning model to displacement mode. This 
process will be termed a sequential wave transition since each of the 
assault waves sequentially performs the mode transition. This is illus¬ 
trated in Figure 3-4. The reason for this transition is due to engineer¬ 
ing stability requirements that this displacement configuration be 
achieved prior to crossing the surf line. The obvious effect of this 
transition is that exposure time to close-in direct-aimed fire will be 
created. 

3. Overall Force Structure 

The model aggregates the combat organizations involved in the 
ship-to-shore phase of the amphibious operation into several homogeneous 
combat units. Each unit is characterized by certain offensive and 
defensive capabilities in comparison to each of the other units. 

Table 3-1 illustrates the combat organizations which have been 
explicitly modeled. The force level of each unit was represented by 
state variables as indicated. The initial force level for each unit is 
input-data to the model. This, therefore, permits the user to investi¬ 
gate alternative wave composition options as well as various defensive 
scenarios without having to make modifications to the model algorithm. 
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Table 3-1. State Variables Representing Combat Organizations 


Combat Organization 

State Variable 

Shore Defenses - Tank Assets 

DT 

Shore Defenses - ATGM Assets 

DS 

Incoming assault Waves of LVA 
representing waves 1 thru 5 

WV(I) 1=1 thru 5 

A cumulative combat force 
comprised of those Marine 
grounc units which have 
arrived on the beach and 
debarked their LVA 

TLF 

P'ire support assets of the 

I amphibious task force 

!_ 

ATFFS 


The tactical combat interactions that exist between these nine combat 
units 'Within the overall force structure are illustrated in Figure 3-5. 

4. Shore-Oefense Scenario 

The defensive scenario utilized in the model includes a force 
comprised of both tank (DT) and antitank guided missiles (DS). Tank 
and ATGM units are emplaced 75 meters inland of the waterline at an 
elevation of approximately 5-10 meters. The model does not explicitly 
maneuver or emplace individual tanks or ATGM systems within each unit, 
but aggregates the cumulative effects of the individual vehicles and 
weapons within each category. 

a. Defensive Unit Force Levels 

The state variables DT and DS represent the total unit 
"strengths" in each of the defensive unit categories. Specifically, a 
DT=3 indicates that within the shore defenses there exists a unit of 
tanks having a total combat effectiveness equivalent to three continuously 
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firing independent weapon systems. A similar interpretation is applicable 
to the state variable D5. 

b. Defensive Fire Allocation 

The two categories of direct-fire weapons are assumed to 
engage targets (incoming LVA) according to a predetermined tactical scheme. 
The defensive "plan" was parameterized as follows: 

(1) Window of Engagement . Each weapon category was assigned 
an engagement window as illustrated in Figure 3"6. Onl> those LVA located 
within the range window could be fired upon by the shore defensive forces. 
The windows are designated by the following input parameters: 



TANK 

ATGM 

Maximum Engagement Range 

TENGMX 

SENGMX 

Minimum Engagement Range 

TENGMN 

SENGMN 


(2) Engagement Rules. Additional defensive tactical criteria 
are implemented into the model logic according to the following rules of 
engagement: 

--A defensive weapon may only engage the two closest incom¬ 
ing waves if more than two waves of LVA are at any time located within the 
weapon's engagement window. 

--If only one wave of LVA is present in a weapon's engage¬ 
ment window, defensive fires of that particular weapon type will be dis¬ 
tributed uniformly against the surviving LVA in that wave. 

--If two waves of LVA are both contained within the engage¬ 
ment window, defensive fires of that particular weapon type will be 
distributed according to a tactical allocation submodel. A weighting 
factor, (DEFWT), input by the user is utilized in establishing the 
proportion of the total weapon strength to be allocated against the 
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surviving LVA's in each of the two waves. Specifically, if DEFWT(l) = 2, 
and DEFWT(2) =1, then each surviving LVA in the closer of the two incom¬ 
ing waves would be allocated twice as much fire as surviving LVA in the 
seaward wave. For example, if waves three and four were both located 
within the ATGM engagement window, then the proportion of OS's fire 
allocated to surviving LVA in wave three would be: 


DEFWT(l) X WV(3) 

--X DS 

DEFWT(l) X WV(3) + DEFWT(2) x WV(4) 

where: WV{3) is the state variable for the 
current number of survivors in wave 3 

c. Attrition-Rate Coefficient Computation 

The classical Lanchester hypothesis for aimed-fire attrition 
is that the casualty rate of a unit is proportional to the size of the 
opposing force. If a Unit "A" is being engaged by a Unit "D", this 
action may be expressed by the differential equation: 

dA 

- = -Betapj X D 

where: Beta^^ is called the Lanchester attrition-rate 
coefficient 

It is assumed that this functional relationship holds for 
each pairing (firing unit, target unit) over the small time interval 
dt. The credibility of the model relates the performance characteristic 
data together with the tactical and physical configurations for each of 
the combat units to derive the attrition-rate coefficients. 

It was decided to express the attrition-rate coefficients as 
the product of the rate of fire (ROF) and the single shot kill probability 
(P(k)) as follows: 
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Betaj^ = P(l<)0A X ROF„^ 

where: DA represents a Unit "D" firing on a Unit "A" 

More complicated models exist [Refs. 9 and 10], however, for the purposes 
of the modeling of the ship-to-shore LVA and defender attrition, this 
method was deemed sufficient. 

Attrition-rate coefficients as described above were derived 
for each pairing (defensive weapon, target) yielding the ten variables: 




I = 1-5 


and 


A switch mechanism is incorporated into the rate of fire (ROF) factor by 
implementing the functional relationship: 


ROF 


D -WV(I) =< 


if WV(I) is located outside 
the engagement window 


if WV(I) is located within 


the engagement window 


where: TBF (Time Between Firings) is evaluated by 

TBF = Aim-Reload Time p.ojecIilUelocity - 

Target Speed 


It should be noted that the relatively slow projectile velocities repre¬ 
sentative of anticipated ATOM assets in the future does cause such 
velocities to become significant in this computation. 

The second factor used in determining each attrition-rate 
coefficient is the single-shot kill probability (P(k)). It is assumed 
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that a hit by a large caliber projectile would constitute a "kill" in 
that it most likely would inflict damage serious enough either to sink 
the LVA, or render it immobile, thus eliminating it from contributing 
to the buildup ashore. A second assumption is that both defensive 
weapon systems addressed would exhibit normally distributed, uncorrelated 
horizontal and vertical errors. Typical dispersion data, both mean and 
standard deviation, for the Tank and ATGM is required as input-data for 
the hit probability computations. 

The suppressive effects of incoming fire upon each of the 
defensive units was considered a significant factor with respect to its 
effect upon the survivability of the incoming assault waves of LVA. It 
was assumed that the suppressive effect would significantly reduce a 
unit's rate of fire, and also increase the error standard deviation. The 
modeling of these suppression effects was accompl.ished by assigning a 
relative suppression factor (SUPFAC) in the interval 1, 2, to both the 
Tank and ATGM units. This factor was determined according to the 
following guidelines; 

SUPFAC = 1 No incoming fires (i.e., the defensive 
unit casualty rate is zero) 

SUPFAC - 2 Maximum incoming fires (i.e., the defen¬ 
sive unit casualty rate is comparable to 
that realized upon full allocation of the 
ATF fire support assets) 

It was assumed that the aim-reload time (^RTM) would be 
increased by approximately 50 percent under the conditions represented 
by a SUPFAC = 2.0. Within the ROF submodel, this is expressed by the 
1 inear relationship: 

ARTM 3 ^jp = ARTM|^Q,^ 3 yp x (0.5 + SUPFAC / 2.0) 
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Additionally, it was assumed that up to 100 percent increase 
in the error standard deviation could be expected under a maximum suppres¬ 
sion environment, hence: 

ERROR STD 5 LIP = error x SUP I-AC 

d. Defensive Breakpoint 

It was assumed that if during the course of the ship-to-shore 
movement phase the defensive forces suffered a cumulative loss in excess 
of 70 percent of their initial force strength, the remaining shore defenses 
would withdraw and commencement of the land combat phase of the battle 
would take place. 

5. LVA Assault Wave Conceptualizatio n 

The model is programmed to handle up to "ive incoming waves of LVA. 
The initial composition of these waves is input by the user by means of 
the variable WVINT. There are no limitations on the number of LVA's 
capable of being in each wave. However, the use'’ is advised that the 

model was intended to model small-unit amphibious operations only. 

« 

a. Wave Posture 

Model functions RNG, HT, and SPD are called upon within the 
model logic to generate the range, height, and speed, respectively, for 
each assault wave as time is incremented throughout the course of the 
ship-to-shore movement phase. The input of tactical employment para¬ 
meters TBW and RD in conjunction with the physical design parameters 
SPDMAX, SPDMIN, HTMAX, and HTMIN for the LVA being evaluated uniquely 
determines the exact range offshore and vehicle configuration (planning/ 
displacement) for each of the five waves. This information then is imple¬ 
mented in the rate of fire and hit probability calculations. 
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b. Ground Forces Ashore 

As each assault wave arrives at the beach, the surviving 
strength of that wave is transferred to the variable TLF (Total Landing 
Force Ashore). TLF represents a ground combat force equal to that trans¬ 
ported by the number of LVA survivors having arrived ashore. Once 
established, the TLF engages the two defensive units allocating its 
fires between the two defensive weapon categories in the same proportion 
as the number of surviving Tanks and ATGM's--that is: 




DT 

X TLF 

DT 

+ DS 



DS 

X TLF 

DT 

+ DS 


The casualty rates applied against the DT and DS state that 

survivor variables are determined by means of the Lanchester aimed-fire 

attrition-rate coefficients WBETAt, r- nr and WBETAt, r nc by 'che 

TLF - DT TLF - DS 

equations: 


dPT 

dt 

dPS 

dt 


I^BETAyj_p _ Qj X TLFp.|. 


-WBETA^LF - DS ^ 


DS 


The computation of these WBETA coefficients is not performed 
within the model utilizing the detailed rate of fire and P(hit) arguments 
described previously, since in the original LVA assault model developed 
by Chadwick, these parameters were considered to be insignificant in 
relation to the overall model. Chadwick assumed his assault model would 
be used as an auxiliary model to a higher-level model, and would receive 
values for these coefficients from that model. 
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6 . ATF Fire Support Conceptualization 

The impact of the amphibious task force's fire support assets 
contribute significantly to the combat effectiveness of the shore 
defense units. Characterizing each of the two defensive force units 
by a simple "located" or "not located" attribute, the attrition rates 
realized by these force units can be simplified substantially by the 
following approach. 

a. "Not Located" Shore Defenses 

At the commencement of the model it is assumed that the 
defensive units DT and DS are emplaced on shore at locations unknown to 
the ATF. The units initially are engaged as "not located" targets by 
area fire for which the following Lanchester area-fire equations are 
applicable: 


dDT 


dt 

= (ALPHAp.|. X 

ATFFS) 

dDS 


= (ALPHAp^ X 

ATFFS) 


The terms in parentheses on the right hand side of the equa¬ 
tions are to be considered a generalized input parameter. The combat 
effectiveness of the ATF fire support assets is also to be considered 
relatively constant during this segment of combat time, and thus it is 
possible to synthesize these input factors by examining the attrition 
losses due to area realized in a previous full-scale model calibration run. 
b. "Located" Shore Defenses 

Once a particular defensive unit has initiated its engagement 
of incoming waves of LVA it is considered located. At this point it is 

43 






dssumed thdt the ATF fire support orgdnizdtion will engage that defensive 
unit through the use of aimed fire. Again it is assumed that the loss 
rate will be in accordance with the Lanchester hypothesis for aimed fire, 
that is: 

dDT 

- = Beta^T X ATFFS 

dt 

dDS 

- = Beta-,- X ATFFS 

dt 

It should be noted that the right hand side of both of the 
equations is to be regarded as synthesized factors to be calibrated 
from a previous high-resolution application. 

C. LAND COMBAT PHASE 
1. Overview 

The land combat phase, like the ship-to-shore phase, has been 
modeled after broad assumptions have been made concerning the type of 
forces modeled and force attrition. These assumptions are quite similar 
in nature to those assumptions made in the ship-to-shore phase of the 
model which would be expected of similar Lanchester-type combat models. 

The first assumption is that of homogeneous forces, which was 
made as a matter of convenience. The defensive forces will be modeled 
as a Tube-Launched, Optical-Guided, Wire Controlled missile (TOW) company 
made up of three TOW platoons located in three separate and fixed defen¬ 
sive positions. Each TOW platoon will be comprised of three TOW sections, 
and have the capability of withdrawing to an alternate position provided 
as input by the user. The aggressor force ashore will consist of the 


44 






consolidated surviving landing force ashore which has been redistributed 
into three offensive units. 

The second assumption is that the aggressor force units will fol¬ 
low three user-defined routes as they advance toward the three defensive 
force positions. Routes can be supplied by the user, or defaulted to 
preassigned routes dictated by the program. However, only three routes 
can be utilized. 

A line of sight (LOS) model written by Professor James Hartman, 
Naval Postgraduate School [Ref. 11], is used, adding great flexibility 
to the modeling of the terrain in the basic scenario. This has a direct 
impact on the probability of detection during any one time period {t,t+dt) 
which is shown as: 


P{Unit i does not detect Unit 
j in a time period t+dt) 

P{Unit i dees not detect 

Unit j in a time per- = x 

iod t+dt) 

P(Unit i does not detect Unit 
j in a time period t,t+dt) 


The first two assumptions provide a basis for applying the third 
assumption which is that attrition of opposing forces will be defined by 
direct-fire Lanchester differential equations similar to the ship-to- 
shore phase, only modeled in more detail. Figure 3-7 provides the scheme 
for the sequence and general flow of events in the model. 

2. LVA Movement Conceptualization 
a. General 

In the original model, aggressor force and defensive force 
locations were modeled in two different ways. Defensive force locations 
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Figure 3-7. Generalized Flowchart for Land Combat 
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were left as user inputs, whereas aggressor force locations had been pre¬ 
determined by the model builder, and could not be altered by the user. 

This allowed a flexibility of modeling defensive positions, tut flex¬ 
ibility was limited because of the method of determining movement routes 
for the aggressors. Glenn Mills provided a user option to the model which 
permitted the user to model a variety of aggressor force movement scenarios. 
This option allows for the choice of attack routes and vehicle speed. In 
addition, the option is highly useful to the unfamiliar user of the model 
since unit locations and attack routes can be initially set to the model's 
default values. Different user selected parameters can be input as the 
user acquires a better understanding of the model's algorithm. 

b. Model 

Three predetermined routes are provided for aggressor force 
movements. Each route is subdivided into 40-meter length intervals, 
since a nonfiring aggressor unit is assumed to move one such interval 
during a time-step of 10 seconds (i.e., average speed of 9 mph). A firing 
aggressor unit is delayed a specified number of time-steps before moving 
to the next interval by the state variable NOD. Each interval in each 
route is represented by its center point coordinates, and by its direction. 
If an aggressor unit enters an interval along its associated route, then 
it is considered to be positioned in the center of the interval, generating 
a possible location error of t 40 meters, since this is the distance 
between two consecutive intervals. 

c. User-Defined Routes 

The user is required to input the original location of each 
aggressor unit, and the locations of each of ten nodes he desires the 
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aggressor unit to move through as it advances on the defensive unit's 
position. This information, along with vehicle speed, is used to calcu¬ 
late route intervals that move the attacking unit through each of the 
designated nodes. A complete route would look like that depicted in 
Figure 3-8. The method used to complete the routes is as follows. 

The straight-line ground distance between the first tv/o adja¬ 
cent nodes (DIST) is calculated as shown in Figure 3-8. The angle between 
the desired direction of movement and a straight west-to-east movement (a) 
is then calculated. Utilizing these quantities and the distance desired 
to be moved during each time-step (DST), the distance to be moved in the 
X and Y direction (XLN and YLN) is now computed as shown in Figure 3-9. 
These distances are added to the coordinates of the previous interval 
endpoint, point C in Figure 3-9, to determine the coordinates of the next 
interval endpoint, point D. This same distance is again added to com.pute 
the coordinates of the next endpoint. Point E. This process is continued 
until the distance from the last endpoint computed to the next node is 
less than DST. This general process is repeated for each pair of nodes 
until the entire route is completed, or the unit's force level is reduced 
to zero or the battle terminated, whichever comes first. 

To insure that all intervals are of equal length, the computa¬ 
tion of the first interval between any two nodes must be considered separ¬ 
ately by taking into account the distance left over from the last 
computation between the previous two nodes. To accomplish this, the 
first interval takes the remaining distance (e) and adds it to an 
interval length (DST - e) for the first interval between any two nodes. 
This insures that each interval along the route is of length DST, which 
is the required length. 
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Figure 3-8. User Determined Routes 
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Figure 3-9. Route Computation 
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3. LOS, Detection, and Fire Allocation 


a, LOS 

The existence of a 1ine-of-s1ght between any two opposing 
units is determined utilizing a 1ine-of-sight model written and program¬ 
med by Professor James K. Hartman, Naval Postgraduate School [Ref. 12], 
and is listed as Subroutine LOS in the land combat phase of the model. 
Professor Hartman's model utilizes a parametric terrain model proposed 
by Needles [Ref. 13], which represented terrain by modeling individual 
hill masses. Each hill is described by a bivariate normal density function, 
and fitted together to form a section of terrain utilizing the following 
information illustrated in Figure 3-10: 


1) (Xc,Yc) - Coordinates of each hill's centerpoint 

2) PEAK - Peak height of each hill 

3) a - Standard deviation corresponding to 

^ the X-axis 


4) 


y 


- standard deviation corresponding to 
the Y-axis 


5) (p) - Rotation factor 


Once the terrain has been "mapped", the existence of a 1ine- 
of-sight can be determined for each pair of opposing units. The informa¬ 
tion required to accomplish this is the location and elevation of each 
unit, as well as the height of the vehicle each unit uses. Professor 
Hartman's model yields the fraction of aggressing Unit A as seen by 
defending Unit B, and the fraction of defending Unit B as seen by 
aggressing Unit A. Figure 3-11 is used to illustrate this. 
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Figure 3-10. Terrain Conceptualization 
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Figure 3-11. Partial LOS Conceptualization 
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b. Acquisition 


The acquisition process was well-modeled in the original land 
combat model devised by Joseph Smoler. The model employs the concept of 
parallel acquisition, whereby the weapon system continuously searches for 
targets, even while engaging other targets. When such a weapon system 
kills its presently engaged target, it immediately can shift its fire to 
a new target, provided that such a target has been acquired either during 
the engagement of the previous target just killed, or earlier [Ref. 14]. 

A general description of the manner in which Smoler modeled target-acquisi¬ 
tion is provided here. However, a more detailed description is provided 
in his thesis. 

The probability that a Unit j is detected by a Unit i at 
time t+dt was modeled for four different combat situations in which the 
opposing forces might find themselves. These situations can be summarized 
as follows: 


Observer 

Not firing (t,t+dt) 
Not firing (t,t+dt) 
Firing (t,t+dt) 
Firing (t,t+dt) 


Target 

Not firing (t,t+dt) 
Firing (t,t+dt) 
Not firing (t,t+dt) 
Firing (t,t+dt) 


The formulas derived to compute the probability of detection 
for each of these situations have a number of common variables, therefore 
their definitions are provided beforehand for clarity; 

P..(t=dt) = The probability that a typical firer in Unit 
"•J i has acquired one or more targets of type 

j by time t+dt 
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Qij(t+dt) = [1 - P.j.(t+dt)] 


QV.j(t+dt) 





The probability that target j is not 
visually detected by Unit i during 
(t,t+dt) provided Unit j does not 
fire during this time interval 

where: 

Sj(t) = 

the number of survivors 
in Unit j at time t 

and: 


the nonfiring detection rate 
of one target in Unit j by 
one observer in Unit j 


QP.j.(t+dt) = (1 - 

The probability that target j is not 
detected by a launch signature during 
(t,t+dt) provided that target j fires 
during this time interval 

where: P. = The probability that one 

observer in Unit i is look¬ 
ing in a direction which 
enables him to detect target j 

and: FR. = The firing rate of one 

^ firer in Unit j 


The first situation occurs when neither the observer nor the 
target is firing during the interval (t,t+dt). This situation allows the 
observer to conduct search operations only, thereby maximizing the prob¬ 
ability of detecting a target in his sector of responsibility, and has 
the target maximizing his probability of not being detected by exposure 
to an observer by a launch signature. Thus, the probability of not 
detecting in time interval (t,t+dt) is 

Q^.j(t,t+dt) = Q.jj(t) X QV^.j(t,dt) 
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The second situation occurs when the target is firing during 
the search interval (t,t+dt) while the observer is conducting only search 
operations. This provides the observer with additional information to 
assist in detection of the target. The observer will detect the target 
by the target's launch signature. Thus, the probability of not detect¬ 
ing in time interval (t,t+dt) is 

Q^j(t,t+dt) = Q^j(t) X (QV,.j(t,ct) + QP^j{t,dt) - 
QV.j.(t,dt) X QP.j(t,dt)) 

The third situation occurs when the observer is firing during 
the search interval (t,t+dt) while the target is maximizing the probabil¬ 
ity of not being detected by not firing during the interval. The 
observer has lowered detection probability by diverting a portion of 
his force to firing on a known targfet. A new factor is introduced 
which will alter the probability of detection, namely the event: 


A = The situation in which Unit j is within the 
field of view of Unit i, with at least one 
of the targets at which Unit i is firing 


This states that Unit j, which is not currently firing, happens to expose 
itself to firing Unit i when firing Unit i is looking and firing on at 
least one other target in j's principal direction. Thus, the probability 
of not detecting in time interval (t,t+dt) is 


Q^.j(t,dt) 


0 


if event A occurs 


g(n) if j is an_aggressor unit 
= <! and event A occurs 


Q..(t) if j is a defending unit 
and event A occurs 


56 






























where: g(n) is an increasing function of n, 

where n is the number of time intervals 
elapsed since time t. 

and: g(0) = QTj(t) 


± g(n) ^ 1.0 for all n 

* J 


The fourth situation occurs when both the observer and target 
are firing during the interval (t,t+dt). In this situation the observer 
has minimized his searching capability, and the target has maximized its 
probability of being detected. Thus, the probability of not detecting 
in time interval (t,t+dt) is 


Q.j(t,t+dt) =< 


Q.jj(t) X QV^.j(t,dt) 

if event A occurs 

g(n) 

if j is an aggressor 
unit and event A occurs 

Q.,.(t) 

if j is a defender unit 

T J 

and event A occurs 


where: 


QV*j(t,dt) 




and: XT- = -• x RF 

I vJ O 

RF = Reduction Factor (the detection rate 
of Unit i has to be reduced since this 
unit fires during (t,t+dt) and the 
search for targets is not as effective 
as for a nonfiring unit) 

syt) = Sj(t) X PTT,,) 


PTT-^ = proportion of Unit i allocated to 
Unit K 


k = (Unit K is engaged by Unit i and Unit j 
is within the field of view of Unit i 
while observing Unit K) 
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If a 1ine-of-sight does not exist between observer i and 
target j, then no accumulation of detection probability will take place 

during the current time interval (i.e., P..(t) will remain the same). 

J 

However, if a 1 ine-of-sight does not exist throughout more than three 
consecutive time intervals, then the P. . is set to zero (i.e., P..(t) = 
0) and the accumulation process will start again from zero if a 1ine-of- 
sight is acquired at a later point in the battle. The motivation for 
this decision rule is seen by the observation that even if observer i 
loses a 1ine-of-sight with target j for a short period of time, he still 
probably has some idea of where to expect the target to reappear, 
c. Non-Firing Detection Rate 

The situations that occurred when the target was in a non¬ 
firing status had detection probability functions that had as a para¬ 
meter X.j, a non-firing detection rate. The manner in which the model 
derives this rate is quite detailed, and deserves attention. 

To begin, each firer in an observing unit is assigned a 
search section (or sector of responsibility) v/hich is characterized by 
two parameters (see Figure 3-12). These parameters are the section 
width (ISECWD), and the primary direction of search (IPRDIR). Further¬ 
more, it is assumed that the search direction within a section of search 
has the following probability density function known as the LIMICON 
Function; 


f(9) = A + B cos 9 -D 4 9 i D 


where: D = 


ISECWD/2 































Figure 3-12. Search Direction 


Y 



Figure 3-13. Observer-Target Scheme 

(A = Observer, B = Target) 
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9=0 corresponds to the observer's 
primary direction of search 

note: A and B are chosen such that 

/ f(0) d9 = 1 

"D 


To determine the probability that observer A is looking in 
a direction which enables him to detect target B, is the value of the 
LIMICON function integrated from an angular value up to 15° on either side 
of the primary direction of fire, specifically: 


,AN6LFT 

P|, = / f(9)d9 = shaded area. Figure 3-13 

ANGRT 


where: 


ANGLFT = < 


ANGLE + 15° 
D 


if ANGLE + 15° < D 
if ANGLE + 15° > D 


and: ANGLE = the absolute value of the angle between the 

the primary direction (IPRDIR) and the ob¬ 
server-target direction (OTANG) 

ANGRT = ANGLE - 15° 


Now, given that observer A is looking in a direction a such 
that ANGRT “ = ANGLFT, the conditional detection rate ( I ANGRT ^ a 
=ANGLFT) is determined by a regression curve [Ref. 15] which is a 
function of the terrain, target horizontal velocity, and the equivalent 
range for a full height target. This detection rate of one observer 
detecting one target becomes 


^ab " ^ ab ' ANGRT 4 “ 4 ANGLFT) x 
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d. Fire Allocation 

Three conditions are necessary for Unit j to be classified 
as a target for Unit i. First, a 1 ine-of-sight must exist between Unit 
i and Unit j. Second, the range between the units must be within the 
maximum range of Unit i's weapon system. Lastly, the probability that 
a detection of Unit j is made by an observer in Unit i in the time period 
t+dt must be greater than O.CO. 

Once these conditions are satisfied, the manner in which fire 
is allocated to a target depends upon how many targets are to share in 
the firepower of Unit i, and what distance exists between i and the 
new target in relation to the other targets under fire. The priority of 
fire naturally will go to the closest target since it is of a greater 
threat to Unit i than the more distant targets. The amount of firepower 
available from Unit i is naturally a function of the percentage of 
surviving force available to fire in Unit i. 

4. Attrition 

Attrition of forces is assessed based upon variable coefficient 
Lanchester equations of modern warfare [Ref. 16]. This method of attri¬ 
tion assessment was used by David Chadwick in the ship-to-shore phase 
of the model, however, in less detail than was modeled by Joseph Smoler 
in the land combat phase of the model. The "extra" detail provided by 
Joseph Smoler is the generation of the conditional probability of a kill 
given a hit. This probability was stated by Chadwick as a user-supplied 
input parameter. 

The restriction of the model to aimed-fire weapons systems and 
homogeneous forces allows for the attrition of forces to be assessed 
using variable coefficient Lanchester equations of modern warfare. 
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The attrition for a defending Unit j is described by the following 
differential equation: 


dSj.(t) 

dt 


(A^-j X PROP..) X S.(t) 


where: S|^(t) 


The number of survivors in Unit k at time t 


Ajj = The rate at which one firer of Unit i kills 
Unit j targets (attrition rate of Unit j by 
one firer of Unit i) 

PR0P|^.| = Proportion of Unit k allocated to fire 

against a Unit 1 


These basic differential equations of force-on-force attrition were 
approximated by the following Euler-Cauchy difference equations: 


S.(t+dt) = Max(0,S.(t) - ZA..(S.(t) x PR0P..)dt 

* W • nJ U * 

for each defending Unit i 

and: 

S.(t+dt) = Max(0,s.(t) - zA..(S.(t} x PROP.Jdt 

0 0 ■ %J ^ * 0 

for each aggressor Unit j 


The manner in which the attrition-rate coefficient A. . is 
derived stochastically already has been discussed in the model 
enhancement chapter, therefore, only a description of how the deter¬ 
ministic attrition-rate coefficient is derived will be mentioned here. 

The attrition-rate coefficient. A.., for each equation is 

• 0 

computed according to the equation: 
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where: T^.. = the time for one firer of Unit i to kill 
one target of Unit j under the conditions 
in the present time interval 


is computed using the Bonder-Farrel 1 formula [Ref. 17]: 




t +1, + ^ - p(Hh) 

^ P(k|h) P(klh) 


+ P(hlh) ~ P(h) 


where: t 

d 


t 


1 


p(h) 

p(hih) 


P(h|m) 

P{klh) 


Time to acquire a target 

Time to fire first round after a 
target is acquired 

Time to fire a round following a hit 

Time to fire a round following a miss 

Projectile's time of flight 

Probability of a hit on first round 

Probability of a hit on a round given 
that the prior round fired was a hit 

Probability of a hit on a round given 
that the prior round fired was a miss 

Probability of a hit on a round given 
that the round fired was a hit 


There are two assumptions of the Bonder-Farrel1 formula that are 
implied by the model. The first assumption is that fire is Markov- 
Dependent in that the probability of a hit of any round depends only 
upon the result of the previous round. The second assumption is that a 
Geometric Distribution describes the parameter P{k[h) in that accumulated 
damage is considered to be negligible. 

The expected value of T.., E[T..], may now be expressed for 

* kJ * w 

each weapon system in the model. It is assumed that for the TOW 
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weapon system P(k h) = 1.0, and P(h m) = P(h h) = P(h), which results 
in the reduced formula: 


E[Ti .] = t^ + t^ + t^ + 


(tm + t^) X (1 - P(h)) 
P(h) 


If the firing weapon system is a tank, then it is assumed that P(klh) = 1.0 
(due to a lack of information), and that t^ = 0.0 (due to the velocity 
of the projectile). Thus, in this case the formula becomes: 

t 

E[Tj^] = +1, +_X (1 - p(h)) 

P(h|m) 


It should be noted that all targets were considered to be stationary 
throughout the attrition process. This is obvious in the case of the 
stationary defending forces, and was assumed to be the case for the 
aggressor forces due to the fact that the hit probability of a TOW 
against a moving target is almost the same as for a stationary target. 

5. Battle Termination 

Two criteria were used as rules governing battle termination. 
The first criterion was the annihilation (zero force level) of one of 
the two forces. The second criterion was that the distance between 
defender and aggressor forces is too small. 

The first criterion is an intuitively obvious reason for ter¬ 
minating the battle, and thus easy to model. However, although the 
reasons for the second criterion might be as obvious, the modeling of 
this is not simple. The manner in which Glenn Mills modeled it was to 
compute the distance between each attacking sub-unit on which casual¬ 
ties were being assessed (i.e., still alive), and each defending 
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sub-unit that was still alive. If any one of these distances between 
active sub-units was too close, the battle was considered to have 
reduced to close-in, hand-to-hand combat. The outcome of this type 
of combat is not currently provided for in the model, and for this reason 
the battle is simply terminated at this point. However, to insure that 
the aggressor units do not pass by the flanks of the remaining defending 
forces ard remain outside termination distance, a check is made of the 
location coordinates of each sub-unit. If any aggressing sub-unit's X 
coordinate places the unit beyond the location of the most forward 
defending sub-unit still in the battle, the battle also is terminated. 

The specification of the distance between forces for battle 
termination is left as a user-input, which provides added flexibility 
of breakpoint distance analysis. In particular, it lends itself to the 
study of optimum breakpoint distances for various weapons systems on 
the battlefield. 
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IV. FUTURE ENHANCEMENTS 


A small-unit amphibious operation combat model has been presented 
in this thesis which emphasizes the simplistic and avoids the abstract 
to provide an understandable and, more importantly, a useable combat 
model for students of combat modeling. However, the combat model pre¬ 
sented has the potential of being developed into a much more refined 
model which could be studied and utilized by more experienced combat 
modelers. Therefore, several enhancements which might be of some benefit 
to the more experienced modeler are mentioned here as possible approaches 
that could be taken in refining the present model. 

A. HETEROGENEOUS FORCES IN THE LAND COMBAT PHASE 

The current land combat phase of the model involves combat between 
homogeneous forces only—that is, each force is comprised of only one 
weapon system type. This type* of force structure was intentionally 
modeled to maintain a relatively simple model to understand and work 
with. However, added flexibility could be attained by modeling multiple 
weapons system types for each of the opposing forces. This would allow 
the user to analyze the effect that different force mixes would have on 
battle outcome. 

The addition of different weapons system types within a single unit 
would reguire extensive restructuring of the attrition process currently 
used in this model. Although Lanchester equations still could be uti¬ 
lized in computing direct-fire weapon system attrition, separate Lan¬ 
chester equations would have to be provided for each weapon system. 
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Furthermore, with the addition of indirect-fire weapon systems (i.e., 
artillery, naval gunfire, and close air support) Lanchester equations 
for area-fire would have to be implemented for each area-fire weapon 
system type. The total attrition of any particular unit then would 
be the summation of the damage assessed by each weapons system type on 
the target being attrited. 

An enhancement of this type would result in more realism at the 
cost of longer execution time, and a more complicated attrition process. 
Since the original intent of the thesis was to provide a simple model 
to understand, it would be advisable to retain a copy of the original 
model prior to adding this enhancement. Then a simple model would 
still be available to the less experienced combat modeling students, 
while a more detailed model would provide the realism that more exper¬ 
ienced modelers would demand. 

B. LOGISTICAL SUPPORT 

Logistical support is one of the most overlooked factors of combat 
in the development of combat models. The influence that the resupply 
of fuel and ammunition alone have on the outcome of a battle is obvious 
and deserves attention. 

Ammunition and fuel consumption could be modeled along the same 
lines as attrition (i.e., through the use of expected values of con¬ 
sumption). When ammunition or fuel on hand reaches a specified critical 
level, a unit could be restricted in movement, or experience a reduced 
level of fighting effectiveness and maneuverability (based on a shortage 
of ammunition and fuel) until resupply of the critical resource could be 
obtained. 
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The amount expended of these resources would necessarily be a function 
of the number of surviving firers in the unit, the number of vehicles 
available to transport the unit, and the number of targets engagec by 
the unit at any one time interval. The expected values of these items 
then could be used in computing the expected rate of consumption of 
ammunition and fuel. Therefore, the overall process could be modeled by 
initially allocating specific levels of these resources (i.e., ammuni¬ 
tion and fuel) to each unit at the commencement of the battle, and 
subtracting the expected expenditure of the ammunition and fuel of a 
particular unit based upon the expected number of survivors firing on 
engaged targets, and the distance traveled by the expected number of 
surviving vehicles of the unit. 

C. GRAPHICAL BATTLE SUMMARY 

A graphical display of what is taking place on the battlefield can 
be worth a thousand words to the user of a combat model. Plotting 
unit locations and force levels on a display of the actual terrairi 
fought upon would eliminate time-consuming interpretation of these 
results from a printed battle summary report. An enhancement of this 
sort would serve both the experienced and inexperienced users of the 
model. The inexperienced user would have results displayed in a much 
more understandable format, while the experienced user would be aole 
to study variant combat scenarios with much less effort and time 
expended. 
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V. FINAL REMARKS 


The purpose of the model that has been developed is to illustrate 
a number of underlying concepts of combat modeling which have been 
addressed in this study. Therefore, it seems appropriate to readdress 
these concepts to allow the reader to reflect upon them in light of 
what has just been presented. 

A. INTEGRATING INITIALLY INDEPENDENT COMBAT MODELS 

The model developed here was made up of two sub-models: ship-to- 
shore and land combat models. These sub-models, as discussed earlier, 
utilized similar combat modeling methodology (i.e., Lanchester equa¬ 
tions) in computing force level attrition. However, each sub-model 
was developed by different individuals, which created several problems 
when the two separate sub-models were integrated into a singular con¬ 
tinuous flow algorithm. In particular, individualized FORTRAN coding 
techniques and documentation of state variables within the program 
structure required the restructuring of major portions of FORTRAN 
code to make the overall combat model tractible and understandable. 
This serves to illustrate the need for a standardized programming 
technique to be applied to programming of combat models, and high¬ 
lights the need for proper planning and coordination in development 
of large scale combat models by teams of combat modelers. 

B. THE USER-ORIENTED APPROACH TO COMBAT MODELING 

This thesis illustrates the desirability of a user-friendly 
approach to combat modeling. It was a major contention of the thesis 
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that this approach to combat modeling has not been closely addressed 
by combat modelers providing combat models for the United States mili¬ 
tary. Furthermore, it was mentioned that the lack of concern given 
to this approach of combat modeling might help to explain the lack of 
enthusiasm exhibited by the United States military in utilizing combat 
models for the training of field commanders and staffs. The thesis had 
as one of its purposes, the presentation of a combat model designed to 
be easily understood and utilized by intended users, combat modeling 
students. Combat models should be designed and documented with the 
user's capabilities and needs in mind, as opposed to those of the 
programmer. 

C. A COMBAT MODEL FOR STUDENT USE 

The small-unit amphibious operation combat model presented here is 
a basic Lanchester-type combat model which has been designed with a low 
level of complexity in order that it might be understood more easily, 
and studied by students of combat modeling. It has been recognized 
that combat modeling students may have little or no experience of the 
governing theory, and therefore would comprehend the theory of combat 
modeling more easily by utilizing and understanding its basic application. 
For this reason, enhancements that would increase the complexity of the 
model are discouraged, and enhancements that would make the model more 
understandable (i.e., graphical battle summary reports) are strongly 
encouraged. 
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I. INTRODUCTION 


The purpose of this manual is to familiarize the user with the model, 
and to provide administrative information describing how the potential 
user would access and run the model. 

The small-unit amphibious operation combat model is a two-phased 
combat model which conducts both ship-to-snore and land combat. The 
model uses both aimed and area-fire Lanchester-type equations for 
casualty assessment. The battle is initiated by an amphibious task 
force positioned 25 miles offshore from an opposing defensive force 
which is illustrated in Figure A-1. If an amphibious landing is success¬ 
ful, land combat will be conducted inland over a 10 x 10 km piece of 
terrain representing an area east of Fulda, West Germany, known as the 
Fulda Gap, which is depicted in Figure A-2. 
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Figure A-1. LVA Approach Conceptualization 
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Figure A-2. Lend Combat Terrain Model 
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II. AVAILABLE OPTIONS 


The model has been developed with a number of options available to 
the user to provide more model flexibility for the more experienced 
user. Each of these options, including user responsibilities, is 
discussed here with the input requirements for each being outlined in 
the next section. 

A. STOCHASTIC VS. DETERMINISTIC ATTRITION 

The user has the option of using stochastic or deterministic attri¬ 
tion computation. Both methods utilize Lanchester aimed-fire equations; 
the difference between the two is the method of calculating the attrition- 
rate coefficients used in the Lanchester equations. 

Deterministic attrition can be thought of as the expected value of 
attrition, and is implemented by using the Bonder-Farrell method of 
calculating the attrition-rate coefficient. A... The stochastic method 

* vJ 

can be thought of as the randomization of attrition, and is implemented 
by using random deviates from a Beta Distribution in conjunction with 
the range of a target to generate individual attrition-rate coeffi¬ 
cients for each unit at each time-step. 

B. VARIANT ATTACK ROUTES 

The user has the option of providing variant aggressor force attack 
routes. The user can utilize the program’s straight west-to-east routes, 
or can input desired altered routes for aggressor force units to follow. 

To select new routes, a user must input the number of nodes desired on 


75 











f 




T 



each of three routes, and the coordinates of ^ach of these nodes. The 
program then will compute routes through each node. The nodes must be 
inputed in order from west to east, and should not create an angle 
between the west-to-east axis and the route direction that exceeds 45°. 

C. ALTERNATE DEFENSIVE POSITIONS 

The user has the option of implementing alternative defensive unit 
locations. This option permits the user to add more realism to the 
model by allowing the defending units to withdraw to alternate positions 
when their primary positions become untenable (i.e., distance between 
opposing forces is too close). This breakpoint distance is determined 
and inputed by the user, and also is used as the distance for battle 
termination 1n the event that the battle reduces to close-in combat (i.e., 
hand-to-hand). The alternative to moving the defenders is to terminate 
the battle when the breakpoint distance is initially reached. 

D. BATTLE SLMMARY PRINT-OUT 

The user has the option of limiting the printed output of the model. 
The user can receive a battle summary print-out at the completion of 
each 10-second time interval, or this information can be suppressed, 
printing out the results only after each phase of combat. 
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III. REQUIRED INPUT 


The small-unit amphibious operation combat model presented in this 
thesis has been provided with a blank data set (see Appendix D) which 
includes each variable of the model requiring input provided by the 
user, and space available following each variable for the user to place 
the desired variable value. However, the definition of each input 
variable may not be familiar to the first-time user of the program. 
Therefore, the following list of input variables and their definitions 
is provided as a quick reference for the user of the model. 


Input Variable 

IPRINT 


SPDMAX 

SPDMIN 

HTMAX 

HTMIN 

WIDTH 

TENGMX 

SENGMX 

SENGMN 

TARTM 

SARTM 


$h ip-to-Shore Phase 

Definition 

User option for selecting type of battle 
summary report desired: 

0 - Each Time-Step 
1 - End of Battle 

Maximum speed of LVA in the water. 

Minimum speed of LVA in the water. 

Height of LVA above water at maximum speed 
Height of LVA above water at minimum speed 
Width of an LVA. 

Tank maximum engagement range. 

ATGM maximum engagement range. 

ATGM minimum engagement range. 

Tank aim-reload time. 

ATGM aim-reload time. 
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Input Variable Definition 


TVEL 

Tank projectile velocity. 

SVEL 

AT6M projectile velocity. 

TSIGV 

Standard deviation error in the vertical 
axis for Tank fire. 

TSIGH 

Standard deviation error in the horizontal 
axis for Tank fire. 

TMEANH 

Bias error in the horizontal axis for Tank 
fire. 

SSIGV 

Standard deviation error in the vertical 
axis for ATGM fire. 

DEFWTS 

Defensive force tactical allocation weights 

WVINT(i) 

Initial strength of assault wave I. 

DINIT(i) 

Initial strength of defensive Tank (1=1) 
and ATGM (1=2) units. 

A(i) 

Aggressor force attrition coefficients. 

B(i) 

Defensive force attrition coefficients. 

WB(i) 

Aimed-fire attrition-rate coefficients 
for defensive force Tank and ATGM units. 

GAINL 

Defensive force attrition level at which 
remaining defending forces withdraw and 
ground assault comniances. 

GAMMA 

Aim-reload time suppression factor. 

DELTA 

Aiming error caused by the suppression 
factor of ATFFS 


The remaining portion of the input data refers to the terrain model 
developed by Professor James Hartman. It is suggested that this portion 
of the data set not be altered until the user has studied and fully under 
stands the Hartman terrain model. 
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Land Combat Phase 


Variable Definition 


ITRIT* 

Input variable denoting whether attrition 
will be stochastic or deterministic: 

0 - Stochastic 

1 - Deterministic 

DSEED** 

Double precision seed used in the Beta 
Distribution Random Deviate Generator. 

PP - QQ 

PD - QD 

Input parameters for the Beta Distribu¬ 
tion Random Deviate Generator: 

PP-QQ Aggressor force 

PD-QD Defensive force 

NBU 

Number of defensive units. 

NRU 

Number of aggressor units. 

RMINTK 

Minimum effective range of an LVA weapon 
system. 

RMAXTK 

Maximum effective range of an LVA v^eapon 
system. 

RMINTW 

Minimum effective range of a defensive 

TOW weapon system. 

RMAXTW 

Maximum effective range of a defensive 

TOW weapon system. 

IRTE 

User option for selecting type of aggressor 
force attack routes: 

0 - Program determined 

1 - User determined 

ISPD 

Speed of aggressor force units: 

1 - 9 MPH 

2-12 " 

3- 15 ” 

4- 18 " 


■f* 

XIC(i,j), YIC(i,i) Coordinates of the interval endpoint 


N(i) 

of the route for Unit i. 

Number of nodes for aggressor route i. 


Note: *There are two ITRIT variables in the data set. The first ITRIT 
refers to the aggressor forces. 


**There are two DSEED variables in the data set. 
refers to the aggressor forces. 
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Variable 


Definition 


XL0C(i,j),YL0C(i,j) 

Coordinates of node i for aggressor route i 

X(i),Y(5) 

Location of defensive Unit i. 

FL(5) 

Force level of a defensive Unit i. 

IPRDIR(i) 

Principal direction of fire of defensive 
Unit i. 

lALT 

User option for selecting alternate 
defensive positions: 

0 - Yes 

1 - No 

BREAK 

Breakpoint distance between aggressor 
units and defensive units. 

ITEM 

Input variable denoting number of time- 
steps allowed for aggressor unit moves. 

XA(i),YA(i) 

Coordinates of alternate position for 
defensive Unit i. 

P(i J) 

Probability of first round hit by Unit i 
in range band j. 

PHH(i,j) 

Probability of a hit following a hit by 

Unit i in range band j. 

PHM(i,j) 

Probability of a hit following a miss by 
Unit i in range band j. 

PKH(i,j) 

Probability of a kill given a hit by 

Unit i in range band j. 
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IV. EXP&CTED OUTPUT 


The small-unit amphibious operation combat model's output is 
designed to be self-explanatory. Each phase of the amphibious opera¬ 
tion is reported in the output of the model. The output format for 
each phase will include an initial information section to provide the 
user with feedback concerning the operation of the model as read-in 
by the model from the user-supplied input data. This serves as a check 
and a record for the user to insure that the model was run according to 
the intended design of the user. Secondly, battle summary reports are 
provided at specific points of the battle depending upon the desires 
of the user as input by the user option variable IPRINT. An example 
of the model's output is displayed in Appendix F. 
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V. ACCESSING AND EXECUTING THE MODEL 


The prospective user who wishes to study the small-unit amphibious 
operation combat model must first contact Professor James Taylor of the 
Operations Research Department and obtain the user identification 
number and password for the disk space containing the model and its 
support programs. 

A. ACCESSING THE MODEL 

Once the required infcrmation is obtained, the user should proceed 
to LOG ON to his OWN disk space entering the CMS mode of operation. 

Upon entering CMS, the following commands should be executed: 

LINK TO (USER ID*] 191 AS 192 RR 

PASSWORD 

ACCESS 192 B/A 

COPYFILE AMPHIB FORTRAN B = = A 
COPYFILE SEA DATA B = = A 
COPYFILE LAND DATA B = = A 
COPYFILE 3SEA DATA B = = A 
COPYFILE BLAND DATA B = = A 
COPYFILE WAR EXEC B = = A 
RELEASE 192 (DET 


*Note: USER ID refers to the user id provided by Professor Taylor. 
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What is r6C6iv6d on the user's disk is a copy of the following files; 

1. The Small-Linit Amphibious Operation Combat Model (APPENDIX B). 

2. A complete data set; SEA and LAND (APPENDIX C). 

3. A blank data set; BSEA and BLAND (APPENDIX D). 

4. The model's executive program: WAR (APPENDIX E). 

B. EXECUTING THE MODEL 

To execute the model utilizing the data se" provided, the user must 
first compile the FORTRAN program, AMPHI6, by entering the following 
commands: 

DEF STOR IM 

I CMS 

FORTGI AMPHIB 

Once the program is compiled, the user enters the name of the 
executive file WAR, which then executes the program and displays the 
listing file of output from the model, (i.e., AMPHIBl LISTING (APPENDIX 
F)) in the BROWSE mode of XEDIT. 

C. ALTERING THE DATA SET 

The user may desire to invoke one of the available options provided, 
or alter specific elements of the existing data set to "play out" various 
combat scenarios. To alter the existing data set, the user first decides 
whether to alter the ship-to-shore phase of combat, or the land combat 
phase. Once this has been established, the user can simply XEDIT the 
appropriate data file, replacing the old input data with the new input 
data. 

To construct an entirely new data set, the user should make use of 
the blank formatted data set provided. The user simply XEDIT's the 
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BSEA or BLAND data files, inputting new data by typing over the spaces 
provided. The variable names are listed in both of the data sets, as 
well as in Chapter III of this user's manual. The space provided in 
the blank data sets is designed to be compatible with the READ font'at 
statements of the program. 

D. EXECUTING THE MODEL AFTER ALTERING DATA 

If the user has just altered specific elements of the data set pro¬ 
vided without altering file names, the user will once again enter the 
name of the executive file WAR, and enter the new data set file names 
where appropriate. Once this editing of the executive file has been 
accomplished, the user simply enters the executive file name WAR to 
execute the model again. 
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VI. PROGRAM STRUCTURE 


The small-unit amphibious operation combat model is a computerized 
model written in FORTRAN. It consists of a main program and 19 sub¬ 
routines. To assist the user in understanding the operation of the 
model, a brief description of the function of each subroutine, as well 
as the functioning of the main program, is provided. 

A. MAIN PROGRAM 

The main program serves as a director program for the model. It call 
for the initialization of data for the ship-to-shore phase of combat, and 
then commences the execution of that phase of combat. The results of 
the ship-to-shore phase of combat as provided by subroutine SEA are 
then reviewed to determine if the land combat phase of combat should 
begin, or if the battle should be terminated. If the results warrant 
a continuation of the battle, the reason for continuation is printed 
and land combat is initiated. 

B. SUBROUTINES 

There are 19 subroutines in the model. The function of each has 
been provided at the beginning of each subroutine in the coded program, 
and also is presented here for clarity. 

1. Subroutine SEA 

This subroutine is the main driver program for the ship-to-shore 
phase of the amphibious operation. Its main purpose is to initialize 
key parameters, and to direct program flow in the ship-to-shore phase 


of combat. 
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2. Subroutine RKINT 

This subroutine provides the interface between the EULER numerical 
integration routine (RKLDE6) and the subroutine ATTR which determines 
each unit's status as time progresses throughout the amphibious operation. 

3. Subroutine ATTR 

This subroutine determines the attrition rates and updates the 
status of each unit with respect to shore movement based upon the given 
state variable strengths, and implements this information into the attri¬ 
tion loss-rate computation. 

4. Subroutine DTGTS 

This subroutine determines the wave numbers that are to be engaged 
by the defensive Tank and ATGM units, based upon the engagement window 
criteria and LVA wave survivor force levels. 

5. Subroutine DATAIN 

This subroutine reads in all user-supplied information required 
by the ship-to-shore phase of the model. 

6. Subroutine OUTPUT 

This subroutine provides an input summary printout based upon 
the data received by subroutine DATAIN. A printout of dispersion data 
generated as a result of data supplied also is provided. 

7. Subroutine PHIT 

This subroutine computes the probability of a hit based upon 
the range, width, and height of a given target. The type of weapon 
being employed against the target then is taken into consideration for 
computing the specific probability of a hit. 
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8. Subroutine INTRP 


This subroutine is a check to insure that the range of a target 
and the dispersion data are compatible for comauting the probability of 
a hit in subroutine PHIT. 

9. Subroutine RATE 

Given the range and speed of a target, along with the type of 
weapon being used to fire upon the target, and the suppression factor 
the firer is being subjected to, subroutine RATE computes the rate of 
fire used against a particular target. 

10. This is the primary subroutine of the land combat phase of the 
amphibious operation. Information required for the operation of the 
land combat phase is read-in and printed in a summary table for user 
review. The information provided by all other subroutines used in the 
land combat phase is used in this subroutine as input to the basic land 
combat algorithm. 

11. Subroutine SETUP 

This subroutine is used to read-in the terrain data and create 
parametric terrain. The terrain data v^ill be used when computing 
line-of-sight between targets and observers, as well as providing a 
grid system for unit locations and movement. 

12. Subroutine ROUTE 

This subroutine computes the route of each aggressor unit when 
the user has selected the option of inputting aggressor routes. It 
calculates the coordinates of each interval endpoint along the route, 
making each interval length (distance moved during a ten-second time- 
step) the same. The interval length is determined by the speed the user 
has selected and inputted for the current battle. 


87 





13. Subroutine LAMBDA 


This subroutine used in conjunction with the lirie-of-sight 
routine computes the detection rate (DETRAC) of target j by the observer 
i, given the percent of target visible (PCTVIS) to the observe>‘, 

14. Subroutine ELEV 

This subroutine determines the terrain elevation for a given 
set of X, Y coordinates. This function is used in conjunction .vith the 
1ine-of-sight subroutine in computing a 1ine-of-sight between observer 
and target. 

15. Subroutine STOCH 

This subroutine determines the attrition coefficients when a 
user has selected a stochastic attrition option. The calculation is 
a function of the original stochastically determined attrition coeffi¬ 
cient, as well as a function of range. 

16. Subroutine ETK 

This subroutine computes the expected time for a given firer 
to kill a given target. The calculation is a function of range, time 
of flight for a round, and hit and kill probabilities for the firing 
weapon system. It is a number that is used in computation of fhe 
deterministic attrition coefficients. 

17. Subroutine SORT 

This subroutine is used to sort targets in ascending range 
order. This is used to determine the priority of a target for fire 
allocation. 

18. Subroutine KOVER 

This subroutine determines what portion of a particular target 
is covered by the terrain between the target and observer. 
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This number is used both in the detection of the target, and in the 
attrition computation. 

19. Subroutine LOS 

This subroutine was written by Professor James Hartman, Naval 
Postgraduate School. It computes a percent of a target visible to a 
particular observer, given the location coordinates of both. 
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VII. DEFINITION OF VARIABLES IN COMPUTER PROGRAM 


A. VARIABLES USED IN THE SHIP-TO-SHORE PHASE 


A(i) 

Aggressor force attrition coefficients. 

ATFFS 

Amphibious Task Force Fire Support. 

B(i) 

Defensive force attrition coefficients. 

CDSURV(i) 

Current strength of defensive Unit i: 

1 - Tank 

2 - ATGM 

CSURVE(i) 

Current strength of assault wave i. 

DA(i) 

Attrition rate for defensive Unit i due to 
the effects of ATFFS/TLF. 

DEFWTS 

Defensive Force Tactical Allocation Weights. 

DELTA 

Aiming error caused by the suppression factor 
of ATFFS. 

DINIT(i) 

Initial strength of defensive Unit i. 

DSl 

That portion of the defensive force ATGM 
unit assigned to engaging the closer of two 
multiple waves in the ATGM engagement window. 

DS2 

That portion of the defensive force ATGM 
unit assigned to engaging the farther of 
two multiple waves in the ATGM engagement 
window. 

DTI 

That portion of the defensive force Tank 
unit assigned to engaging the closer of two 
multiple waves in the Tank engagement window. 

DT2 

That portion of the defensive force Tank 
unit assigned to engaging the farther of 
two multiple waves in the Tank engagement 
window. 

DTI PH 

Hit probability of rounds fired by DTI against 
the assault wave in its engagement window. 
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DTIROF 

Rate of fire utilized by DTI against the 
assault wave in its engagement window. 

GAINL 

Defender attrition level at which remaining 
defending forces withdraw and land combat 
commences. 

GALF 

Denotes whether the landing force buildup 
is sufficient for land combat: 

0 - Insufficient 

1 - Sufficient 

GAMMA 

Aim-reload time suppression factor. 

GATK 

Denotes whether the landing force has 
initiated the land combat: 

0 - Not started yet, 

1 - Started already. 

GATM 

Time at which land combat commenced. 

IL(i) 

Denotes if wave i has reached the shore: 

0 - Wave i not ashore, 

1 - Wave i ashore. 

IPRIMT 

Denotes whether the user desires battle 
summary at each time-step, or just a 
final summary: 

0 - Battle summary printed after each 
time-step, 

1 - Final battle summary only. 

rwp?: 

Weapon-type code: Tank = 1, ATGM = 2. 

IWSTAT(i) 

Current status of assault wave i: 

0 - Not engaged, 

1 - Landed, 

2 - Under fire by ATGM, 

3 - Under fire by Tank, 

4 - Under fire by both ATGM and Tank. 

RD 

Distance offshore at which waves initiate 
their transition. 

RKSURV(i) 

Concatenation of CSURV and CDSURV. 

SA(i) 

Attrition rate for wave i due to ATGM. 

SARTM 

ATGM aim-reload time. 

SENG(i) 

Wave number of the closer of two assault 
waves in the ATGM engagement window. 


91 




























SENGMN 

ATGM mini mum engagement range. 

SENGMX 

ATGM maximum engagement range. 

SRNG(i) 

Firing range to wave SENG(i). 

SSIGH 

The standard deviation error in the 
horizontal axis for ATGM fire. 

SSIGV 

The standar'd deviation error in the 
vertical axis for ATGM fire. 

SVEL 

ATGM projectile velocity. 

SWTS(i) 

The proportion of the total defensive 
force ATGM strength to be allowed to 
engage wave SENG(i). 

TA 

Time first assault wave initiates its 
transition. 

TA(i) 

Attrition rate for assault wave i due 
to Tank fire. 

TARTM 

Tank aim-reload time. 

TB 

Time first assault wave completes its 
transition. 

TBW 

The interarrival time between waves 
arriving at the beach. 

TFF 

Time first assault wave reaches the beach 

TENG(i) 

Wave number of the closer of two assault 
waves in the Tank engagement window. 

TENGMX 

Tank maximum engagement range. 

TMEANH 

The bias error in the horizontal axis 
for Tank fire. 

TMEANV 

The bias error in the vertical axis 
for Tank fire. 

TRNG(i) 

The firing range to assault wave TENG(i). 

TSIGH 

The standard deviation error in the 
horizontal axis for Tank fire. 

TSIGV 

The standard deviation error in the 
vertical axis for Tank fire. 
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TSURV 

Total number of surviving LVA ashore at 
the current time. 

TVEL 

Tank projectile velocity. 

TWTS(i) 

The proportion of the total defensive 
force Tank strength to be allowed to 
engage wave TENG{i). 

WB(i) 

Aimed-fire attrition-rate coefficients 
for defensive force Tank and ATGM assets. 

WID 

Width of LVA. 

WVINT(i) 

Initial strength of assault wave i. 

WVRNG 

Firing range to assault wave i. 


B. VARIABLE USED IN THE LAND COMBAT PHASE 


ALPHA(i) 

Initial attrition-rate coefficient for 
stochastic attrition option. 

■ ANGH(i) 

Orientation angle of the hill ellipse 
measured in degrees counter-clockwise 
from East to the major axis. 

AP0A(i,j) 

The average proportion of the j^^ attacker 
of Unit i allocated to fire on Unit i. 

AVD 

Average distance. 

AVSP 

Average speed of moving aggressor units. 

BASE 

Overall terrain elevation above sea level. 

BREAK 

Breakpoint distance between aggressor units 
and defensive units. 

DISMAX 

Maximum distance allowed between aggressor 
units before the leading unit is delayed. 

DIST 

The straight-line distance between two 
movement nodes input by the user. 

DST 

The distance in meters to be moved each 
time-step by aggressor units. 

ECC(i) 

The eccentricity defined as the ratio of 
major axis length to minor axis length. 
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FL(i) 

Force level of Unit i. 

FO(i) 

Initial force level of Unit i. 

lALT 

Denotes whether the user desires alternate 
defensive positions or not: 

0 - Yes, 

1 - No. 

IC 

Counts number of time units a defender 
has been moving. 

IDIR 

t h t * 

Direction of interval in i^*^ route. 

II(i) 

Interval index for Unit i. 

IITIME 

Current time. 

I MAX 

Maximum number of time intervals allowed. 

IMOVE 

Number of time units a defender is allowed 
for moving to an alternate position. 

IPRDIR(i) 

Primary direction of fire for defensive 

Unit i. 

IRAN 

Range. 

IRTE 

Denotes whether user wants to input routes 
or not: 

0 - Program determined routes, 

1 - User determined routes. 

ISE 

A switch variable set to 1 when the defen¬ 
sive force AT6M unit initiates its fire. 

ISECWD(i) 

Width of search sector for defensive Unii i 

ISPD 

Input variable to denote user's desired 
speed for aggressor force movement; 

1 - 9 MPH, 

2- 12 MPH, 

3- 15 MPH, 

4- 18 MPH. 

IT 

Current time period. 

ITE 

A switch variable set to 1 when the defen¬ 
sive force Tank unit initiates its fire. 

ITEM 

Input variable denoting number of time- 
steps allowed for aggressor unit move. 
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ITIME 

Current time, in seconds, of battle. 

ITRIT 

Input variable denoting whether attrition 
will be stochastic or deterministic: 

0 - Stochastic, 

1 - Deterministic. 

lUSTAT(i) 

Current status of Unit i: 

0 - Alive, not firing, 

1 - Alive and firing, 

2 - Killed, 

3 - Moving. 

LVAFR 

Firing rate for LVA weapon system. 

LATOB 

Indicator variable for one- or two-way 

LOS calls: 

0 - Do not compute LOS from Unit A 
to Unit B, 

1 - Compute LOS from Unit A to Unit B. 

LISTH(i) 

List of hill numbers for each grid square. 

LOA(i,j) 

"h h 

The number of the target of Unit i. 

LOST(iJ) 

Denotes whether 1ine-of-sight exists 
between Unit i and Unit j. 

LOT(i,j) 

h 

The number of the target of Unit i. 

LST 

Index number for the first hill listed 
for grid square (i,j) in LISTH(i). 

MVTDIR(i) 

Movement direction of Unit i. 

N(i) 

Number of nodes inputted by user for 
route i. 

NA(i) 

Number of aggressors of Unit i. 

NBU 

Number of defensive force units. 

NOVELS 

Number of forest ellipses in terrain. 

NF(i) 

Number of time units a Unit i is allowed 
to fire at the same location. 

NHL(i,j) 

Number of hills in each grid square (i,j). 

NHILLS 

Number of different hills to be modeled. 

NHTOT 

Total number of hills modeled on battlefield 
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NLOSC(i,j) 

Number of continuous time-steps that LOS 
does not exist between Unit i and Unit j. 

NOD 

Number of time intervals Unit i delayed 
in movement. 

NOI(i) 

j. u 

Number of intervals in the i^^ route. 

NRU 

Number of aggressor force units. 

NT(i) 

Number of targets of Unit i. 

CFL(i) 

Force level of Unit i during previous 
time-step. 

F(i ,j) 

Probability of first round hit by Unit i 
in range band j. 

FHH(i,j) 

Probability of a hit following a hit by 

Unit i in range band j. 

FHM(i,j) 

Probability of a hit following a miss by 
Unit i in range band j. 

PKH(i,j) 

Probability of a kill given a hit by 

Unit i in range band j. 

FM 

The proportion of time a moving unit is 
searching for targets. 

FOA(i ,j) 

t h 

The proportion of the attacker of 

Unit i allocated to fire on Unit i. 

FOL(i) 

Percent of Unit i lost during the current 
time-step. 

FTT(i) 

Proportion of surviving firepower allocated 
to the i^^ target if there are j targets 
available to be engaged. 

RANGE 

Current minimum distance between aggressors 
and defenders. 

RKATTR 

Vector containing the current attrition 
loss rates to be applied within the 

Euler integration routine. 

RF 

Detection rate reduction factor for a 
firing unit (in comparison to a non¬ 
firing unit). 

RMINTK 

Minimum effective range for an LVA 
mounted weapon system. 
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RMINTW 

Minimum effective range for a TOW weapon system. 

RMXTK 

Maximum effective range for an LVA mounted 
weapon system. 

RMXTW 

.Maximum effective range for a TOW weapon system. 

ROF 

Rate of fire. 

ROT(i,j) 

Range of the target of Unit i. 

SIZETK 

Size of LVA weapon system. 

SIZETW 

Size of TOW weapon system. 

SPRD(i) 

Measure of hill size which is defined to 
be the distance in meters measured along 
the major axis from hill center to the 
contour line which is 50 meters down from 
the peak. 

SUM BO 

TotaT defensive force level. 

SUMRO 

Total aggressor force level. 

SUPFAC 

Suppression factor. 

TA(k) 

j. i* 

Time tc acquire a target for k^^ weapon 
system type (k = 1, 2). 

TFl(k} 

Time of flight to lOOOm for k^^ weapon 
system type (k = 1, 2). 

TF2{k) 

lime of flight to 2000m for k^*^ weapon 
system type (k = 1, 2). 

TF3(k) 

Time of flight to 3000m for k^^ weapon 
system type (k = 1, 2). 

TH(k) 

Time tc fire a round following a hit 
for weapon system type (k = 1, 2). 

TI(k) 

Time to fire first round after target has 
been acquired for weapon system type (k = 1,2). 

TM(k) 

Time to fire round following a miss for 
weapon system type (k = 1, 2). 

TMACIJMACJ 

Elevation of Unit i and Unit j in LOS model. 

TOWFR 

Firing rate for TOW weapon system. 
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TPOL(i) 

Total percentage of lost since battle began 
for Unit i. 

VISFR(i,j) 

The fraction of Unit i as seen by Unit j. 

VISFRA 

Fraction of Unit A as seen by Unit B. 

VISFRB 

Fraction of Unit 8 as seen by Unit A. 

X(i),Y(i) 

Coordinates of Unit i. 

XA(i),YA(i) 

Coordinates of alternate position for 
defensive Unit i. 

XC(i).YC(i) 

Coordinates of center of hill i. 

-< X 

1—1 1—1 
o o 

Coordinates of the j ' interval endpoint 
of the route for Unit i. 

XL,YL 

Distance added to previous interval endpoint 
for vehicle to move DST during a time-step. 

XL0C(i,j) 
YL0C(i J) 

Coordinates of the node inputted by 

the user for the route of Unit i. 


98 








APPENDIX B 


COMPUTER PROGRAM 
for the 

SMALL-UNIT AMPHIBIOUS OPERATION COMBAT MODEL 


The small-urn't amphibious operation combat model is a computerized 
model written in FORTRAN. It consists of a main program and 19 sub¬ 
routines. It was designed to serve as a reference to itself in order 
that the reader would not be forced to refer to various manuals outside 
of the program each time an explanation of the functioning of a parti¬ 
cular aspect of the program was desired. A listing of the computer 
program follows. 
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ooo or\o r^ooo 


This PRCGPAM IS A SMALL-UMT APPHIBICUS CPERATICN CCMBAT MODEL 
UTILIZING LANChESTER-TYPE EQUATIONS TO COMPUTE ATTRITION. 

IT CONSISTS OF TWC BASIC PHASES. THE FIRST BEING THE SHIP-TO-SHOR 
CCPBAT PHASE, ANC THE SECOND BEING THE LAND COMBAT PHASE. 

SHIP-TO-SHCRE PHASE CCMMON BLOCK VARIABLES ***** 


COMMON /AMPH/IL(5 ) ,HB(2> , / (2 1 ,B{2),ITE,ISE,RO,WVINT(SI,WID, 

♦ TBW,DINIT(2),GAINL,H«STAT(5) 

COMMON /ENGR/ SPCM AX,SPOM IN.HTMAX,HTMIN,TTS,TAA,TB,TFF 
COMMON /OISPER/T?IGV(6,2) ,TSIGH(6,2),TKEAMH(6,2), 

*SSIGV(7,2i,SSIGH(7,2) 

COMMON /CcF/TENGMX,S£NGMX .SEN6MN,TARTM,SARTM,TVEL, 

#SVeL,OEFWTS( 2) 

COMMON /SUPEFT/GAMMA,DELTA 
COMMON /lOUT/TSUFV,ICRINT 

***** LAND CCMBAT PHASE COMMON BLOCK VARIABLES ***** 

COMMON /GRPl/ IPP0IR(6),I IECWC(6) ,MVTDIR(6) ,X(6) ,Y<e ),SP0(6I 
COMMON /GRP2/ TA(2 ) ,T1(2) ,TH( 2 l,TM(2 I ,TF1(21,TF2{2) , TF2(2) , 
*P(2,6),PHH(2,e),FHM(2,6.) ,fKH(2,6) .Tr(21 
COMMON /GRP3/ NBU,NRU,FL< (),FC(6) ,NOI(3) ,XIC(3 ,2001 ,YIC(3,200), 
*I0IR(3,2C0I,AVSF,ISPC 

*, IUSTAT(n,II(6),L0ST(6,6 ) ,VISFRA.VISFPB,SIZETK, 

*S IZE'^W, NT ( 6) ,NF( 6 I ,SRF,OI SMAX , 

*NLCSC(6,6I,VISFR(6,6),PMIHTK,PMXTK,RM INTW,RMX~W,OP,tqWFR,LVAFR, 
*PTT(3,3) ,RF,PCA(e,6) ,AFOA 16,61 ,L0A(6 ,6),NA(6) ,OFL (6) , POL(6) 
COMMON /GRP4/ TFCL(6l,CLDC(6,cI,Q(6,6) 

COMMON /COPS/ LOT (6 ,6) ,R0T(6,6 1 

COMMON /HILLS/ XC( ICC) ,YC (luC ) ,PEAK( 100) ,ANGH(1001 ,SPRO(IOO) 
CQMM.CM /HILLS/ SC C < 100 I, P >X ( 1 CO I , PY Y ( 1 CO ) , PXY ( 100 ) , BASE 
COMMON /HILLS/ NHILLS . 

COMMON /CCV5R/ CXC(1501,CTC(1 :C),CPEKk( 1501,CPXX(i50) ,CPYY(150) 

CQMMCN /CCVER/ CPXY(150 I,HCVELS 

CDMMCN /CCLNTR/ KH .KhV, ,KV ,KN , HGRS ,KE LL ,K I NT 

COMMON /CR ID/ LS7(5,4) ,NHL(5,A),L1STH(1501,KHREP(150) ,KTREP 
COMMON /GRID/ LSTC(5,4),NC(5,4),L1STC(400 I,KCREP(150 I 
COMMON /GRP6/ ALPHA)6) 

COMMON /CRP7/ XA(6 ), YA(6 ) ,IMCVE(61 

main DRIVER PROGRAM ***** 


GATM=0. 

GATK=0. 

c *** initialize data fcp ship-to-shcre phase 

CALL DATAIN 
C 

CALL SETLP 

C *** CCNCUCT SHIP-TC-SHCRE COMEAT PHASE 
CALL SEA (GATM.GATk ) 

IF (GATK.NE.O.) GC TC 5 
HRITE(6,600) 

STCP 

5 IF<GATK-2.C) 1C,2C,3C 
10 WRITE(6,610I 

GC TC 40 

20 WPITE(e,62C) 

GC TC 40 

3C WRITE(6 ,63CI 

40 WRIT£{6,640) GATM 

C *** CCNCUCT LAND CCMEAT PHASE 

CALL GRCLNDC GATM ,TSURV ,IPRINT,TTS) 

C 

C 

60C FORMAT) IX,'TOTAL LANDED LANDING FCPCE STRENGTH IS INSUFFICIENT', 
* ' FOR GROUND ATTACK'1 

610 FCR)-'AT) 1 X ,'LAND CCMBAT STARTS WHILE SHCRE COMBAT IS GOING ON') 
620 FORMAT)IX, 'LAND CCMBAT STARTS AFTER DEFENDER BREAKS CONTACT') 

630 FORMAT i lx ,'L ANC CC’^BAT STARTS A=TER ALL WAVES LANDED'I 
640 FORMAT)/IX,'LAND CCMBAT ATTACK TIME =',F6.1,' SECONDS') 

STCP 


TOO 
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ENC 

SIBROUT I.NE SEA (GATN* ,GATK J 

SLEROUTINE IS ThE ^^A IN DRIVER PROGRAM FOR THE SHI P-TC-SHORE 
Ph^SE OF THE AMPFIBICUS OFERATICN. ITS MAIN PURPOSE IS TO 
INITIALIZE KEY PARAMETERS AND TO DIRECT PROGRAM FLCW FOR THE 
SUF-TC-SHCRE PHASE CF COMBAT 

CCMMCN /AMPH/IL{5),wa(2) , A{2) ,B(2J»ITE,ISE»RDtWVINT(5) ,WID, 
*TEVtDlNn(2) ,GAIM,IV.STAT (5) 

CCMMGN /ENGR/ SPCMAX,SPDM IN,HTMAX»HTMIN,TTS,TAA*T8,TFF 

CAIL OUTFIT 
IRC=500 
ITEW=12C 
RC^l.O^IFD 
TeW=1.0*ITBW 
TINT =0.0 


♦ CCMPUTATICN OF FIRST WAVE TIME PARAMETERS 
TA-TIME FIRST WAVE INITIATES TRANSITQN 
TB-TIME FIRST ^AVE COMPLETES TRANSITION 
TFF-TIME FIRST UAVE REACHES THE BEACH 

TAA=(5C0C.-RD)/SPCMAX 
TE=TAA+TTS 

TFF=TB+ ( FC-( 0. 5^ ( S F DMAX-S FDMI N )=»TTS ) -150. )/SPOMIN 
DEL=10. 

WFITP{5,6C0) RCtTEW 

600 FCFMATI/,lXt • ITERATION IN IT I A TED. ..R D= •»FIO.3,1X,» TBW = 
, F10.3) 

CALL RKINTIOEL ,TINT»N»GAT M,GATK ) 

RETURN 

ENC 


SUBROUTINE RKINT{H,TI»NfGATM,GATK) 

SUBROUTINE RKINT FPCVICES THE INTERFACE BETWEEN 
THE EULEF NUMERICAL TMTEGFATICN RCUTI\E(RKLDEQ) 

ANC THE SLBRQUTINE ATTR WHICH CETERMINES EACH 
UNIT'S STATUS AS TIME PRQCRESSES THROUGHOUT THE 
AMFHIBICLS OPERATION 

COMMON /AMPH/IL{5) ,WE(2) , M2 ) r B(2), ITE,ISE,RD,WV INT(5 ) ,WID, 
^TBHtOINITi2J ,GAINL,IWSTAT (5) 

CCMMCN / ICUT/TSURV t IPR!MT 

OIMEMSICN CSUPV(5) fCCSURV (2) ,TA{5),S A{5)tDA(2) 

DIMENSION RKSURV{7) ,RKAT1R{7 ) ,TA^TR (2C0,12) ,TI ME(200) 

VARIABLE CEFINITICNS 

IMAX - MAXIMUM ALLCWAdEL NUMBER OF TIME INTERVALS 
IL(I) - A SWITCH VARIABLE wHCSE ELEMENT I IS S 5"^ TO 1 WHEN 
WAVE I ARRIVES AT THE BEACH 
ISE - A SWITCH VARIABLES SET TC 1 WHEN THE DEF.ATGM 
UNIT INITIATES ITS FIRE 
IT - CURFENT TIME PEPICO 

ITE - A SWITCH variables $ET TC 1 WHEN THE DEF.TANK 
UNIT INITIATES ITS FIRE 
T - CURRENT TINE 

TSLRV - TOTAL NUMBER OF SIRVIVING LVA AT THE CURRENT TIME 

STATE VARIABLE DEFINITIONS 
CDSURV(I) - CURRENT STRENG‘"H CF DEFENSIVE FORCE I 
I = 1 TANK 

I = 2 A1GM 

CSURV(I) - CURRENT STRENGTH GF ASSAULT WAVE I 
DIMT(I) - INITIAL STRENGTH QP DEFENSIVE FORCE I 
RKSURVd) - CONCATENATION CF CSURV ANC CDSURV 
WVINT(I) - INITIAL STRENGTH CF WAVE I 
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XO 

15 

2C 

25 

30 

120 


40 

45 


46 

50 


110 

130 

65 


I> = 2850 
G^LF^O. 

IV^X=1S9 

nE==o 
ISE = 0 
TSIRV-0* 

7 i^E( n-^c. 

T = TI 

DC 10 1=1 ,5 

CSURVd ) = WVINT( n 
TSUR\ = TSUPV+C StFV (I ) 

U( I ) 

jWSTA7(I)=C 
CCMINU6 
DC 15 1=1.2 

CDSURV ( I) =C INIT( I ) 

CCMINUF 
DC 20 J=l,12 
T^TTRx i »'J )=0« 

17 = 1 

CC 25 I=lt5 
PK5UPV (I )*-=CSljRV^ I) 

RKSLRVC 6 )=CSLR VI1) 

RKSURV(7 )=CSURV(2) 

CC 30 1=1,7 
RKiiTTP ( I )=C. 

M=0 

C^IL A T7F(T,CSLRV tCCSURV,7A,SA,CA,GA LF,GATK,GATM,IX) 

S'^^TE VARIABLE DEFINITIONS * = 

0A(1) - 4TTRITirK RATE FCP DEFENSIVE UNIT I DUE 
THE EFFECTS CF ATFFS/7LF 
S4(I) - 41TRITICN PATE FG F WAVE I DUE TO ATGM 
TA(I) - 4TTRITI:N rate FCF wave I DUE TO TANKS 


TO 


RK4TTR(I) IS A /ECTCF CCN7AIMNG THE CURRENT ATTRITION 
LOSS RATES TO Be APPLIED W^hIN THE EULER 
INTEGRATION ROUTINE TO THE STATE VARIABLES. 
1=1,5 LV4 WAVES 1-5 

1 = 6 OT 

1=7 OS 


I> = IX 4 7 

IF(IL(1).LC.99) GC TC 130 
00 4C 1=1,5 

FKSURV n ) = CSLRV( I ) 

RKATTF (1) = (TA ( n+SA( I) )*(-!.0) 

CC 45 1=] 2 

RKSUPV I 1+5)=CCSURV( I) 

RK dTRI I 45 l = -l .O^CA ( I ) 

S=PKLOEG / T ,RKSLkV ,RKATTR , T,H ,M) 

LS INK = 1 

DC 50 1=1 ,5 

IF (LSINK.EC.CJ GC TO ^6 

CSLRV(I) = CSUPV(I) + RKATTR(I) 
GC TC 50 

eSURV (II = PKSLRV(I) 

CCiNTiNUE 
DC 55 1=1,2 

COSUPVxII=RKSURV( 1 + 5) 

CC^TINU£ 

IF(S-i. I no , 12C, 12C 
WRITE (6 ,600 ) 

STOP 

CCNTINUE 

IT=IT+1 

TSIRV=0. 

CC 65 L=l,5 
T51RV=TSIPV+CSLPV (L) 

IF(TSURV.LE.C.I TSLRV=0. 

TI.ME( IT ) = T 

IFdPRINT.EO.l ) GC TC 999 
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C PPINT RE5LLT CF 5h IP-TC-SFOR E i^OVEMENT AFTER 

VjPITE(6,dlC5 I 
^FITE(6 ftZC) 

CC 70 1=],^ 

PL0ST = 1.-CSLRV( I ) /WVINK I ; 

WRI TE (6 f6 30 ) I ,CSUPV( n. IWiTAT (I) , OLOST 
7C CCMINUE 

FLCST = i .-CSURV (5 )/\^VINT( 5 ) 

WRITE(6,640J C SUFV(5},lUSTAT(5}.PL0ST,TSURV 
FLCST=1*-CCSURV( : I/OINHl 1) 

WFITF(6,6ECJ CCSORV(1J,PLCST 
FLCST= l.-CCSUR V( 2 )/CINr\ 

T ASURV = CCSURV ( 1)-►CCSUR V( 2 ) 

WRITE (6,660 CCSLRV(2 ) , ? L CST « T ASU RV 
SSS CONTINUE 


:ACH TIME STEP 


DETERMINE R: THE FIRING R4NGE TO THE LAST INCOMING ASSAULT WAVE. 


R = PNG(T-4.*T BW ) 

DETERMINE IP ALL WAVES LA^DcC ANO LANj COMBAT STARTED 
NOTE: THE MODEL IS TERMINATED IF: 

1. THE PIPING RA^GE TC THE LAST ASSAULT WAVE IS LESS 
THAN 75 METER I. 

2. THE DEFENSIVE BREAKPOINT HAS BEEN REACHED 

2. THE MAXIMUM NLMBEF CF ITERATIONS HAS BEEN EXCEEDED 


IF(R.LT.75.) GC TC 200 

IF(IT .GT. IMAX J GG TC 2C0 

IF ( ILiII.EC.991 GO TO 2CC 
GC TO 120 

200 N=IT 

*** PPINT RESUL^ OF THE SHIP-TO-SHORE COMSAT PHASE 
WFITE<6,6iOJ T 
WPITE(6 , 62CJ 
CC 90 1=1,4 

FLGST = l.-CSLRV(n/WVlMT( n 
WRITE (6,6 20 J i ,CSLR V( I J , I WSTAT (I I ,FLCST 
90 CONTINUE 

FLCS'r = l . -CSURV (5 )/WVlNT( 5 I 

WP ITE( 6 , 6401 C EuF V( 5 I , IWS TAT( 5) , => LCS T , TSUR V 
PLCST=l.-CC3UPV(i )/ClNlT( 11 
WPITE(6,65CJ CCSLFV(1J ,=LCST 

FLCS'r=i.-CDSURV(2)/0INIT( 2 ) 

TASUPV=CCSuRV(1)^CCSlRV( 2 I 
WFITE(6,6cCJ ClSLRV( 2J ,PL CST , TASURV 
WRIT = (6,67CJ T SURV 
IF(GATK.Gc.l. ) GC TC 9999 

IF( TSOP V. LT.9. I GC TO =999 
GATK=3. 

GATM=T 

9999 CONTINUE 

600 FORMAT(IX,» ERPCR..S .N 8.1. CR • 2 » J 

61 C FORMAT (///IX, »TIME = ’ ,F6. 1 .IX, 'SECCMDS*// ) 

620 format ( 1>, »v/AVE* ,2X, ‘ccrc E LE VEL • ,2X , ' ST A TUS * , 2X , • LOST-PCT» , 
*^2X,»TGTAL SURVIVING* I 
620 FORMAT(3X,I 1,3>,F1C.A,EX,I 1 ,5X ,F8.3J 
64 0 FOFMAT(2X.*5 * ,3>,F10.4,5X ,11,fX,-8.3 , 7X,F5.2) 

650 FOFMAT (1X,*TANK’ ,2X,=10«*+ ,11X,FS.3J 

66 C FCFMAT (lX,*ATGM*,2X,FlC.4aiX,r8.2,7X,F5.2l 

670 FCFMAT ( IX, *F INAL LVA SURVIVORS AShORE= ’,F10.3) 

RETURN 

END 


FUNCTICN PKLOEC(^ ,T, P,X, H ,NT ) 
OIMcNSICN T( 1 ) ,F (1 ) ,C(23 ) 
NT=NT+1 

GC TO (1 ,2,3 ,4 ),i\T 
1 hl = H 

AA = H1/4.C 
CC 11 J=1,N 
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11 0(J)=0. 

X=X+AA 
GC TO 5 

2 ;<s>'«'AA 
GC TO 5 

3 X=X+AA 
GC TO 5 

A CC 6 L=1,N 

6 Y { L)=Y( D + ^iA’^F (L) 
N7 = C 

X=X+AA 
PKLDE0=2 . 

GC TO 8 
5 CC 7 1=1,N 

7 >r{ I)=Y( I ) + AA«F (1) 
RKlOEQ=l.C 

8 RETURN 
ENC 


SUEROUTINE A TTR ( T t CSUR V, 0 EUR V , TA , SA , CA,GALF ,GA^;< ,GA1 , IX ) 

GIVEN THE CURRENT TIME ANC STATE VARIABLE STREf^GTHS, 

SUEROUTINE ATTR CETERMINE5 THE ATTRITION RATES AND UPDATES 
THE STATLS OP EACF UMT WITH RESPECT TC SHORE MOVEMENT 
ANC IMPLEMENTS THIS INFCR^ATIC^ INTO THE ATTRITION LCSS RATE 
CCMPUTAT ICN. 

♦ STATE VARIABLE CEFIMTICNS 

CA(I) - CURRENT ATTRITION LOSS RATE FOR OEF. FORCE I DUE TC 

/TFFS(AMFHielOUS TASK FORCE FIRE SUPP0RT)/TLF EFFECTS 
IL(I) - ^•HEN ECLAL TC 99 INDICATES THE DEFENSIVE BREAKPOINT 
FAS BEEN PEACHED 

SA(I) - CLRRENT ATTRITION LCSS RATE FCR WAVE I DUE TC ATGM FIRE 
TAd) - CURRENT ATTRITION LCSS RATE FCR WAVE I DUE "“C TANK FIRE 


CCMMGN /AMPH/IL(5),W8(2) ,4(2) ,6(2),! TE,ISE,RD,WVI NT(5) ,WI0, 

♦ T EV,DIN IT(2) , G4 INL,IWS7AT (5) 

CCHMGN /CEF/TEFGMX,SENGMX ,SENGMN,TARTM,SARTN,TVEL, 
♦SVEL,DEFWTS( 2 ) 

INTEGER TENG(2) ,SENG(2) 

CIMENSICN -^RNG (2) ,TWTS (2) ,SRNG(2) ,DSLRV(2) ,SWTS(2) , 

♦ CSLRV(5 ) ,TA( 5 ) ,S4 (5 ) ,CA(2 J, ASX20 ) 


LSINK = 1 


CC 10 1=1 ,5 
TA( I )=C, 

SA(I)=C. 

1C CCNTINUE 

VARI/BLE CEFINiTIQNS 

DTI - THAT PCRTICN GF THE OT UMT ASSIGNED TO ENGAGING THE CLOSER 
OF TWO MULTIPLE WAVES IN THE TANK ENGAGE^^ENT WINDOW 

DT2 - TH4T PCRTICN CF THE DT LMT ASSIGNED TO ENGAGING THE FARTHER 
OF TWO MULTIPLE WAVES IN THE TANK ENGAGEMENT WINDOW 

OSl - THAT PCRTICN OF THE OS LMT ASSIGNED TO ENGAGING THE CLOSER 
OF TWO MULTIPLE WAVES IN THE ATGM ENGAGEMENT WINDOW 

CS2 - THAT PCRTICN OF THE DS LMT ASSIGNED TO ENGAGING THE FARTHER 
OF TWO MULTIPLE WAVES IN THE ATGM ENGAGEMENT WINDOW 


CS1=0^ 

CS2=0. 

DT1=0. 

CT2=0. 

FAC=1. 


0 
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CE7ERMINE IF P4R7 CF LANDING FORCE ADVANCE TO ATTACK INLAND 
KEY TERR/ IN 

IF{GATK.EC.1.0) GC TC 15 

IF(GALF.EC.1.C. /NC.(CSLRV(1)+DSURV(2) ) .LE.GAINL* (DINIT{1) 
»+CINIT(2)}) GATM=T 

IF (GALF.EC.l.O.ANC. (DSURV(1)^0SURV{2)) . LE.GAINL*(DINIT{1) 
*+CINIT{2))) GATK=i.O 

DETERMINE IF DEF. BREAKPOINT FAS BEEN REACHED 

15 IF((OSURV(l)+DStPV(2))*LT.O,2*(DINIT{l)+DINIT(2))) GC TO 20 

CETERINE ATTRITICN RATE C^ DEFENSIVE FCRCES 3Y ATFFS 
BASED UPCN AREA CR AIMED FIRE STATUS 

IpVPNG = FIRING RA^G£ TC AN ASSAULT WAVE 

c/(i)=e(!) 

DA(2) = e(2) 

IF( ITE.EC.C) DA( i ) =:A( D'^O 5URV {1) 

IF( ISE.EQ.O) D A (2 )=A{ 2 )*DSLRV{ 2) 

GC TC AO 
20 CSLRVC 1 

CSLRV(2)=C. 

DA(1)=C. 

CA(2)=0. 

IF(GATK.EC.l.) GC TO 35 
GAT=T 


TIME=* fFT.lf* SECONDS *) 


DETERMINE !<= C E F. ER E AK PC I NT HAS BEEN REACHED BEFORE SUFFICIENT 
LANDING FCPCE IS BLILT UP ON THE SHOPE FCR INLAND ATTACK 

25 DC 30 I=1 f5 1 

WVRNG = RNG (GAT-Tew*{ I«1 )) 

IF{WVFNG.LT.75. ) IL (I )=l 

IF ( IL{ I ) .EC.: ) TLF=TLF + CSURV(I ) 

30 CONTINUE 
GAT=GAT+1C. 

IF (TLF.LT.9.0.ANC.IL(5). EC.l ) RETURN 

IF(TLF.LT .9.0.AND. IL( 5 J.NE. 1) GO TC 25 
GA 1K=2 . 

GAlF^l. 

GATH=GAT 

^RITE{6,61C) GATV 

610 FORMAT (/ax ,'LAND C CMBAT INITIATED 

35 IL(1)=99 

WPnE(6t62C) T 

620 FORMAT ( IxaePEAKPOI NT REACHED AT TIME =• fF7.lt* SECONDS') 

RETURN 

SLEROUTINE DTGTS DETERMINES THE FIRING STATUS FOR THE 
TWC DEFENSIVE LMTS. 

AO CALL DTGTS{T fTENGfTRNGfTWTSfSENGf SRNGfSWTSfCSURV) 

state variable definitions <*<•*>'* 

TEN'G(l) - THE ^^AVE NUMBER OF ThE CLOSER 0^ TWO 
V^AVES IN THE TA^K ENGAGEMENT WINDOW 
TRNG(l) - the FIRING RANGE TC WAVE TENG(i) 

TWTSd) - the FFCrCRTICN CF THE TOTAL DT STRENGTH 
TO BE ALLCCATED TO ENGAGING TENG(l) 

TENG(2) - THE WAVE NUMBER OF T^E FARTHER OF TWO 
WAVES IN THE TAfK ENGAGEMENT WINDOW 
TRNG(2} - similar IN*^ERP R FTAT I CN AS TRNG(l) 

TWTS(2) - SIMILAR INT E RP R £T A"^ ICN AS TWTS(l) 

SENG(l) - THE WAVE NLMBER CF THE CLOSE? OF TWO 
WAVES IN the ATGM Ei\GAGE*^ENT WINDOW 
SPNG(l) - FIRING range tq WAVE SE.NG( 1 ) 

SWTSd) - THE FRCPCRTICN CF THE TOTAL CS STRENGTH 
TO BE ALLOCATED TO ENGAGING SENG(l) 

SENG(2) - THE WAVE NUMBER CF THE FARTHER OF TWO 
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WAVES IN THE ATGM ENGAGEMENT v^INOOW 
SRNG(2) - SIMILAI= IN TERP R ETATI CN AS SRNG(l) 

SWTS(2) - SIMILAP interpretation as SWTS(1) 

CETERMINE THE CUMULATIVE hUMBER OF SURVIVING LVA’S 
TF^T have REAChEC THE EEACH - TLF 

TLF=0* 

CC 45 J = 3 ,5 

IF( IL (J). EQ.ll TLF=TLF4CSURV( J) 

45 CONTINUE 

CETERMINE IF TLF BUILT UP IS SLFFICENT FOR LAND CCMBAT 

IF(TLF*GE*9. ) GALF=1. 

IF(GATK.EO.1.) TLF=TLF-9. 

♦ ALLOCATE THE FCRCE STRENGTH OF TLF BETWEEN THE ^wO 
CEFENSIVE FORCE UMTS 

CSUM=OSUFV(I )4.CSLRV( 2) 

TIF1={0SLPV{1)/CSUMI^TLF 
TLF2={DSLRV( 2) /D5LM)-^TLF 

ACC TO CAl ANC CA2 THE ATTRITION LOSS PATE DUS 
TC THE EFFECTS OF TLFl ANC TLF2 

OA(l)=OA(l)-**TLFl=^W°(U 
CA(2)=0A(2)+TLF2^W0( 2 ) 

IF(DSURV(l).LE*O.OI CA(1)=0. 

IFCOSLPVI2) .LS.C.C) DA<2)=0.0 
ICA = 1 
N = 20 
MUL = 1 
ISCRT = C 

call LRNC( I)( »ASX,Nf NULt ISCRT) 

CETERMINE IF TKEPE EXISTS AN INCOMING WAVE IN THE 
TANK ENGAGEMENT VINCCW (I.E. TENG(1 ) .NE.0) 

IFCTENG(l).EC.C.) GO TO ICO 
ITE = 1 


CETERMINE THE TI^'E SINCE WAVE TENGCl) CROSSED THE 
5CCC. MEIER OFFSHORE MARK -T1 

Tl=T-TeW4 (TENG (1 )-l ) 

CT3=TWTS{1)^0SLRV(1) 

FAC=1. 


CETERMINE THE SUPPRESSION EFFECT ^*0 EE IMPOSED 
CN THE CT UNIT EASED CN THE ATTRITION LOSS RATE 
CLPRENTLY IN EFFECT 

SUPFAC = ATFFS SUFRESSICN FACTOR 

SLFFAC = OA( 1) 

CALL RATE (TRNG (1 ), SPC(T1 ) ,1,SUFFAC,DTIROF) 

OTIROF - RATE CF FIRE LTILIZED BY OTl AGAINST WAVE 
CALL PHlT{TRMG(i),WlC,bT( 71), ItSUPFACfOTiPHl 

DTlPh - HIT FFCeAeiLlT^ OF RCUNOS FIRtC BY DTI 
AGAINST WAVS TENG(1) 

DETERMINE THE ATTRITION LCSS RATE FOR WAVE TENGCl) 

CLE TO CTl FIRES 

TA(TENG( 1 ) )=0TlPH-'^DTlPCF4tCTl 
IF (LS INK.EQ.O ) GO TC 55 

IF (AS;«(ICA).GT.TA(TENG(1 ) ) ) GOTO 50 
TA(TENG( : ) ) = 1. 0 
CC TC 55 

50 TA(TENGCl)) = C.O 


TENGd ) 
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55 IC4 = IC^ + 1 

DETERMINE IF TFEPE IS 4 SECOND INCOMING WAVE THAT 
IS IN the tank engagement WINCCW, if THERE IS THE 
ATTRITION PATE CC PUT AT ! 0 NS ARE SIMILAR IN FORM 
TC THOSE FPEVICUSLV rcRFCFMEC FOR THE CLOSER WAVE 

IF(TENG(2).EC.G) GC TO IOC 
T2=T-iev»*^( TENG(2)-1 ) 

0T2 = TV^TS( 2 ) ^DSlRVI 1) 

CALL PATE(TFNG(2) »SPO(72) ♦ 1, SUPFACtDT2ROF) 

CALL PHIT(TRNG(2) ,WI0,FT(T2),1 ,SUPFAC»DT2PHJ 
TA(T ENG(2) )=CT2FHiOT2RT2 
IF(LSINK.EC.O) GC TC 6f 

IF(ASX(ICA).GT.TA(TENG(2))) GO TC 60 
1A(TENG(2)) = 1*0 
GC TO 65 

60 TA (TENG( 2) ) = 0.0 

65 ICA = ICA 1 

DETERINE IF THERE EXISTS 4N INCOMING WAVE IN THE ATGM 
ENGAGEMENT WINDOW, IF THEPE IS, DETEFMIfJE THE ATTRITION 
EFFECTS AGAINST THAT WAVE OUE TO ATOM THE ATTRITION 
RATE CCMFUTATICN ARE SIMILAR IN FORM TO THOSE FOR THE 
EFFECTS CLE THE TANK FIRE. 


ICC 


IF(SENG(1).EC.C) GC TO 20C 
ISE=1 

SI=T-TeW*(SENG(ll-l) 

DS1=SVTS(1)»DSLRV{2) 

SUPFAC=CA(2 ) 

CALL PATE(SRNG(1),SP0( SI) , 2 ,SUPFAC,DS IROF) 
CALL FHIT(SFNG(1) ,W ID , FT(T1 ),2,SUPFAC,OS1PH) 
SA(SENG(iJ )=0S IPH^OSIRCF^DSI 
IF(LSINK.EC.U) GO TC 7 5 

IF(ASX(ICA).GT.SA(S5NG(1))) GO TC 70 
SA(SENG( 1 ) ) = 1. C 
GC TC 75 

SA ( SENG( 1) ) = C. C 
ICA = ICA + 1 

If (SEMG(2).EC.C) GC TC 20C 
S2=T-ie^^( SENG(2)-L) 

DS2=SWTS(2)^05UPV(2) 

CALL PATE (SPNG(2I ,SPC( S2 ) , 1 , SU PF A C , CS 2R0 F ) 
CALL FHIT(SRNG(2J , W10,FT(S2 ) ,2,SUPFAC,0S2PH) 
SA(SENG(2 ) )=CS2Ph-=CS2RCF*CS2 
IF(LSINK. EQ.G) GO TC 35 

IF(ASX(ICA).GT.SA(SENG(2))) GO TC 80 
SA(SENG(2) ) = 1.0 
GO TC 85 

80 SA (SENG( 2 ) I = 0. 0 

E5 ICA = ICA + 1 


7C 

75 


2CC RETURN 
END 


SIBROUTINE DTGT 5 (T,TENG,T FNG,TWTS,SENG,SRNG,SWTS,C$URV ) 

GIVEN TFE CURRENT TIMS AND LVA WAVE SURVIVOR POPULATIONS, 
SLERCUTINE OTGTS DETERMINES THE WAVE NU 1 M 6 ERS THAT ARE 
TC BE ENGAGED BY DEFENSIVE TANK ANO ATGM UNITS BASED 
CN THE E^GEMENT WINDCW CRITERIA 

COMMON /AHPH/IL(5).WS(2) ,A(2) ,6(2) ,ITE,ISE,RD,WVINT(5) ,WID, 
♦ TeV^,DINIT( 2 ),GAlNL, I/<STAT (5 ) 

COMMON /CEF/TENGiO ,SENGMX ,SENGMN,TARTM ,SARTM,TVEL, 
♦SV£L,DEFWT5(2) 

INTEGER lENG(2 ) ,SENG (2 ) 

DIMENSION TP.NG(2) ,5RNG(2) ,TWTS(2) ,SWTS (2) ,CSURV(5) ,DEMC(5 ) 
00 iO I=l»2 
TENG( I )=0 
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TV»TS( I )=0. 

TRNG( I )=0. 

SPNG(I)=0. 

SENG( n=o 
SWTS( I )=0. 

10 CCMINUE 
J1=0 
JS^O 

SSIM=0. 

CC 100 1=1,5 

NVRNG = rNG(T-TEy**( I-l) ) 

IF(WVPNG.LT.75.) !L(!)=1 

IF (VVRKG.LT.7 5. ) U-'STAT f I)=l 

♦** IF THE PIPING PANGE TO A VAVE IS LESS THAN 75 METEPS, 

THE WAVE IS CCNSICEREG TC HAVe REACHEC A COVERED AND 
CC^CEALEC POSITION CN THE BEACH 

IF( (WVRNG.GT.TENCMX) .. C P. (CS UR V ( I ) .LT . 0.05 ) .CP . 
♦(WVPNG.LT.75.).CP.(JT.GE.2)) GC TO 5C 

JT=JT4l 

TENG(J1 )=I 

TWTS( JT) =OEFWTS( J" )*CSURV( I) 

TSLH = TSUH' + TWTS( JT ) 

TRNGIJ1)VRN G 

50 IF (( WVPNG.GT.SENGNXI .i:r. ( CSURVCI ) .LT.O. 05) .CP. 

♦ (KVRNG. n.SENG^N ) .CP. ( JS. CE.2 I ) GC TO 100 
JS^JS+1 
SENG(JS ) = I 
SPNG(JS)=WVP 

SWTS( JS) =OEFhTS( JS)*CSUP V(I) 

SSUM=SSUH + SWTS(JS ) 

IOC CONTINUE 

CETEPMINE WAVE STATUS 

CC 20 1=1,2 
CC 25 J = 1 ,5 

IF( IWSTAK jl.NE.l.A.’^O.SENGd). EC.J) IWSTAT(J)=2 
25 CONTINUE 

2C CCMINUE 
DC 20 1=1,2 
CO 35 J=l,5 

IFdWSTAK J).EC.l) GO TC 35 

IF( IWSTAT( J ). £0. 2.ANC.TENG( D.EO. J) IWSTAT( j)=4 
IF dWSTAT( J ) .NE.2 .AND .IE N'G( I ) . EG. J ) !WSTAT( J)=3 

25 CONTINUE 

30 CONTINUE 

IF (TENG(1).EO.C) GC TC 500 
CO 2CC 1=1,2 

TV IS (I ) =TWTS(I)/TSU M 
CONTINUE 

IF (SENG( 1 ).£0.0) RETURN 
00 6C0 1=1,2 

SWlSd ) = 5WTSd )/SSUM 
CONT INUE 
RETURN 
ENC 


200 

500 


600 


SUEPOUTINE OATAIN 

SUeROUTINE OATAIN RE^DS AIL USER SUPPLIED INFCPMATICN 
PECUIRED EY THE SHIP-TC-SFOR £ PHASE OF THE MODEL 

CCMMON /AMPH/IL(5),WB(2) ,A(2),B(2),ITE,ISE,RD,WVINT(5),WIO, 
*TeV,0INn(2),GAINL, IWSTAT (5) 

CCHMON /ENGR/ SPCMAX,SPCM IN,HTMAX,HTM IN, TTS,TAA,TB,TFF 
CCMMCN /CISPER/TSIGV(6,2) ,TSIGh(6,2) ,THEANH(6,2) , 
*SSIGV(7,2) ,SSIGH(7,2) 


108 










CCVMON /CEP/TE^Gyx »SHNGMX tSENG'^N* TAR TMtS ARTw, VVEL , 
♦ SVELtOEFVTS( 2) 

CCNMON /SLPEFT/GAMMAtOELTA 
CC^^'CN / ICUT/TSUFV 11 FR INT 
REAO< 5» 5(C) IPRINT 

SPCMAXtSPOMIN ,HTMAX,HTMIN,V<ID 
ITS 

TENG^Xt SENGMX ,S£NGNN 
TAFT-''»SA-TM,T'vEL tSVEL 
( (TSIGV( ItJ) 1 1 = 1 ♦£ ) tJ= 1t 2) 

((TS IGH( ItJ> . 1=Ij 6 ) ,J = 1, 2 ) 

( (TM6/!NH { I t J J ,1 = 1 »6 ) t J=1 t2 ) 

( ( 55IGV(I.JJ , 1 = 1T 7 ) tJ=l,2) 

( ( SS IGH( I* J ) t I=i, 7 ) ,J = 1, 2) 

(DfcF\^TS( I) Tl = lr2) 

(UVINK I ) t 1= 1 ,5) 

<CINIT( I >, I=i ,2) 

(Ml) ,I=lt2) 

(6(1; ,I=it2J 
(hE(I ),I=lt2 ) 

GA IM 

___ GAPMAtCELTA 

5CC FCRMAT(IFX ,I1 ) 

510 format( 10,F6.2t9XtFA. 1,8>,F5. lt7X,F 5.1t 
520 FCFMAT( lEXtrS. ItlEXt ^a.l, 16X,F£.1) 

530 FCRMAK 7> ,7F 10.3) 

540 FCFMAT(3:>tFd. ItlAXtFS.l) 

550 FCFMAT(27Xt5FlC.5) 

560 FCFMAT(40tF<:.5tICXtP9.r) 

570 FCFMAT(17X,F7.2, i2X,F7.1 ) 

RETURN 
ENC 


Re4C(5,510 
PEA0(5 ,5 10 
RcAC(5,520 
READ(5,510 
BEAO( 5, 530 
PEAC(5, 530 
RE/C(5,530) 
REA0(5, 530 
RE/C(5,530) 
PE;0(5,54C) 
PEAC(5,550) 
FE/C(5,540 ) 
PEA0(5,5£C) 
PCA0C5,5(CI 
RE40(3,560) 
REAO( 5,570 
RE£D(5,57C) 


7X,F6.3) 


SUERCUTINE OUTPUT 

*** SUeRQUTUE OLTFIT FRCVIOES AN INPUT SUMMARY INTOU'* BASED UPON 

THE DATA RECEIVEC BY SL3RCUTINE DATA IN. A PRINTCUT CF DISPERSION 
CATA GENERATED AS A RESULT OF DATA SUPPLIED IS ALSO PPCVIDEO 

COMMON /AMPH/IL( 5 ) ,Wa( 2) , A ( 2 ) , 3 ( 2 ) , I TE ,I SE , W VI NT ! 5 ) ,WI 0 , 

♦ Tew,DlNIT(2) tC-AUL, InSTAT (5) 

CCMMCN /C I SP ER/T SIC-V(c,2 ) ,T SIGH (6 ,2) ,TMEANH(6f 2) , 
★SSIGV(7,2),SSIGH(7,2) 

CCFMCN /ENGk/ SPCMAX,SP0M iNtHTMAXtHTMIN, TTStTAA,TB,TFF 
CCMMOM /CEF/TENGMXtSENGMX ,SENGMN, TAR TM rSAP. TM ,T VEL , 

*SVELtO£F^TS(2) 

CCMMON /5UPEFT/GAMMA,DELTA 


INPUT SUMMARY PRINTOUT 




UPITE( 6,E 1C) 
UFITE(6,E12) 
URITE(6,E13) 

UPITE(6,E14) 

UPITE(6 tEil) 
URITE( 6, E15) 
URITE(6,Ei6) 
URITE(6,E17) 
UR IT6(6,E 16) 
UFITE(6tcl9) 
URITE(6,E2C) 
UPITE(6,E21) 
WRITE(6,622) 
URITE(6,E23) 
UFITE(6 ,E24) 
UPITE(6,E25) 
URITE(6,E26) 
UPITE(6,627) 
U P ITE ( 6 , E 2 E) 
hPIT5(6,E29) 
UPITE(6 ,630) 


(UV INT( I),1= 3,5) 
(DINIT(I) ,I=1,2) 


SFOMAXtSPCMI N,HTMAX,HTM IN,WIC 


TENGMXtTARTM ,TyEL 
SENGMX tSENGMN,S ARTM, SVEL 
l)EFUTS( 1) ,DEfUTS (2) 

A(i),6( i) 

A(2),E (2 ) 

We( 1 ) ,W8(2) 
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c 

c 

c 


WRITE(6,631) G4INL 
WR ITE(6,632) G^NM;i,CEL7A 


0I3PERSICN CATA FRINTOtT 


C 


IDISP=1 

IF (IDISF.EQ.C) RETURN 
WRITE(6,633 ) 

VkRITE (6 ,624 ) 

CQ 55 1=1,6 

WRITE (6,63 5 ) TS IGV(1,1 ),TS I GV( 1,2 ),TSIGH(1,1) ,TSIGH(I,2) , 
*TMEAMH(1,1),TWEANH(I,2) 

55 COMINUE 
WRITE(6,626) 

CO 56 1=1,7 

WRITE(6,637 ) SSIGV( 1, 1),SSIGV(1,2 ) , SSIGH(I,1) ,$SIGH(I,2) 

56 COMINUE 
WRITE(6,63E) 

W R i T t: t 6,6 3 9) 


610 FOFMAK » 1«,1X,’=4'^ INITIAL SHI F-TO-SH CRE PHASE INFCRMATION >*») 

612 FOFf^AT (/6X,* INITIAL FORCE STRENGTH*) 

613 FOFMAT (1> ,»WAV£»,:X,*1»,5>,*2*,5X,*3»,5X,*4»,5X,*5») 

614 FOFMAT( 1), 'LVA'Tf (2X,F4. 1 n 

611 FORMAT(/IX,*CEF, TANK ASSETS = ’ ,IX,F^. 1, 5X, 

*»DEF. ATGM ASSETS =»,1X,F^.1) 

615 F0FMAT( / nx, *LVA ENGR SPEGS*) 

616 FORMAT(1> , *S PCMAX * ,2X,* SP CM IN * ,5X,* HT^AX*,2X, » HTMIN *,3X,’WID» ) 

617 F0FMAT( 2>,F5*2 ,3X ,F3*2,3X iF5*2,3X,F4«2,2X,F4«2) 

618 F0RMAT(/1X,*CEFENSIVc TACTICAL PARAMETERS*) 

619 FORMAT(ICX,* RANG8 * ,^X,* A IF-RELCAO *,3X, *PPOJECTIL=* ) 

6 2C FORMAT ( c > , 'M A X * , 3 X , * v I N * , ^X , * ^^ I ME * , 7 X , * VELCCI TY * ) 

621 FORMAT ( IX , *J ANK * , IX , F6.1, «X, F C . 2, 7X , F6.2 ) 

622 F0FMAT(1>,*ATGM* , IX , F6.1 , IX,Fc . i,2X, F5-2,7X,F6.2 ) 

623 FORMAT(/,iX, *CEFENSIVE ^ACTICAL ALLOCATION WEIGHTS;*, 

♦ /,1X,*WA\;E 1 = * , F5^2,1X , *«AVE 2 = *,F5.2) 

624 FORMAK/lx,* DEFENS! VE FCRGE ATTRITION COE FFI CI ENT S » ) 

62 5 FORMAT (9X, *ALPFA'*^a * , lOX , * SET A ^ A * j 

62 6 F0FMAT(1X,*0T* ,6X,R7.5,9X ,F7.5 ) 

627 FORMAT ( 1 X, ’D S • ,6X , F 7. 5 ,9 X ,F 7, 5 ) 

628 FORMAT(/,* AIMEC FIRE ATTPITICN RATE GCEFFICIENTS FOR*, 

*/,* DEFENSIVE TANK /NO ATGM ASSETS*) 

629 format(/ lX,*W8£TA(iJ = *,F7.5,ix,*ABETA(2)=* ,F7.5) 

630 FORMAT (/IX,* BREAKPCINT ASSUMPTICN: 0.3-MTGTAL DEF FORCE)*) 

631 FORMAT(/,IX,*OEFENCER ATTRITIGN LEVEL ALLCWiNG FOR LANC COMEAT* 
’^/,1X,F5.2, •'«(TCTAL CEFENDER FCRCE)*) 

632 FORMAT (/IX,* ARTN SUP '^AC T CR= * , F5. i, 2 X, * ERROR SUP F AC TOR = * , F 5.1 ) 

633 FCRMAK /1X,*DI SPERSICN DATA*/) 

634 FORMAT (2X , *RANGE* , 2X , *TSICV* ,2X,*RAN GE*,2X,* TSIGH * , 

^2X,'RANGE *,2X,*TNEANh» ) 

635 format ( lX,F7oi ,2X,F5ol,lX ,F7. 1 ,1X,F5.1 ,1X,F7. 1 ,1X,F5.1) 

636 format(/3X,* range *,2X,* SS IGV’ ,2X,*RANGE*,2X, *SSIGH* ) 

63 7 FORMAT ( 1 X , F7.1 ,2X , F 5.1 ,1 X ,F7.1 , IX ,F 5 . 1) 

628 FOPMAT(/,* CURRENT STATUS OF WAVE I VARIABLE DEFINITIONS*,//. 

^* C - NOT ENGAGING *,/,* 1 - LANDED*,/,* 2 - UNDER FIRE BY ATGM* 
3 - UNDER FIRE 2V TANK*,/* h - UrjDER FIRE BY BOTH ATGM G TANK 
639 FCFMAT ( * 1 * ,*'^^^^’^ THE SHI F-TO-SHORE PHASE BECI.NS v:«;:>vv*///) 
RETURN 
END 


SLBROUTINE PHIT(PANGE,W,h , IWPN,SUPFAG,PRHIT) 

C GIVEN Tl-E RANGE, WIDTH AND HEIGHT OF A TARGET, AS WELL AS, THE 

C TYRE OF WEAPON EEING USED, SUBROUTINE PH IT COMPUTES THE 

C PPCBABILITY OF A HIT 

C 


CCNMON /AMRH/IL(5) ,W0(2) , 4(2) , E(2 ), ITE,ISE,RD,WVINT(5),W10, 
*TeW,OiN IT(2),GAINL,IWSTAT(5) 

CCMMCN /CISPER/TSIGV(5,2 ) ,TSIGH(6,2) ,TMEANH(6,2), 

^£SIGV(7,2),SSIGH(7,2) 

CCMMON /SLPEFT/GAMMA,DELTA 


C 

C IWRN CODE: TANK 


ATGM =2 


no 
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*** 

50 

100 


60 

70 




c 


10 


20 

30 
C 6CC 
610 
62C 


V^RUSLE CEFI.MTICNS 

ISIGH - THE STC CEV ERROR IN THE HORIZCNTAL FOR TANK 
TSIGV - THE STC CEV ERROR IN THE VERTICAL FOR t^i^k 
TVEANH - THE EIAS ERROR IN THE HORIZONTAL FOR TANK 
THEANV - THE EIAS ERROR THE VERTICAL FOR TANK 
SSIGV/SSIGH - Sli^ILAR INTERPRETATIONS FOR THE ATGH 

PI = ARCGS (-1. C) 

IF(RANGE.LT.25.) STOP 

IF(IWPN.EC.l) GO TO 50 
ATGM FIRING DATA CCMPUTATIGNS 
V<HEANH=0 .0 
HMEANV=C.C 

CALL INTFPISSIGVtRANGEiWSIGVt?) 

CALL INT FF(SSIGV,PANGE*WS IGHt7 ) 

TANK FIRRG OATA CC^PUTATICNS 
kHEANV=C.C 

CAll INT FP (i iVEAKrfRANGEtWr'EANHf6) 

CALL INTFP(TSIGV,PANGE ,i*.S IGV t6 ; 

CALL INTRP(TSIGH,RANGEtWSIGH,6) 

CONVERSICN TC MLS 
Z = ARSIN(H/RANGE) 

WSIGV=wS IGV^( l,+CELTfl'*SUP FAC) 

V»SIGH=WS IGH^ ( 1 . + CELTA^SUP FAC) 

TGTH=( Z^ ^^CO.C )/< 2.C*PI) 

TGTlrt=(APSIN(W/RANGc) ) ^(64 CO.0/(2.0*P I ) ) 


INSTITUTE NORMALITV ASSUMCTICNS TC COMPUTE HORIZONTAL 
AND VERTICAL HIT FRO EA EI L IT I E S 

C = -1.0'*$QRT( 1./2. ) 

HCR1=((lGTW/2.)-H^EA^H)/WSIGH 

H0R2 = ( ( (-I.O^TGTVn )/2.0)-^ ^EA^H)/W SIGH 

FHITX=1.C 

IF(ABS(HCRl).GT.E.) GO TC 60 

PHITX = C.5>( ERFC(C'^H0R1 )-E R F C ( C’‘HO R 2 ) ) 

VER1 = ( ( 7GTH/2. )-v>'.mEANV )/WSIGV 

VER2 = { ( (-1.0^TGTH)/2.)-wHEANV)/'«SIGV 

PHITY=1.C 

I F(ABS( VtPl) .GT.8.) GG TC 70 

Phi TY = C.5--( ERFC ( G-VER 1 )-£ R FC {C^’^VE R2 )) 
PPHT=PHITX^'PH ITY 
RETURN 
END 


SLBPOUTI^E INTPP(X ,ARGtVAL,N ) 

SLERCUTUE INTPP IS A CHECK SUEROUTINE TO INSURE THAT RANGE OF 
TARGET ^NC OISPERSICN OATA ARE CO^PATABLE FOR PRCBABILITY CF 
HIT GGMFUTATION IN SLEROUTINE PHIT 

CINENSIGN X(N,2) 

WFITE(6,6C0) AFC 
IF(ARG.LI.X( 1,1)) GO TO 30 
CO 10 1=1,N 

IF(ARo.GT.X(I+l,l)) GO TC 10 
D1FF=X(I + 1, !)-;{( I, 1 ) 

OELTA= ARG-X{1,1) 

VAL = X( I,2) + (CELTA/CIFF)=»(X( I + l ,2)-X(I,2) ) 

RETURN 
CCNT UUE 

IF(ARG.G1.X( N, 1) ) GO TC 2C 
VAL=X (N ,2 ) 

PETUP^ 

VsRITE(6,aC) 

STCP 

ViP ITE{ 6,620 

FQFHAT ( 1 ) , »ARG»‘)'*^^ = ' ,F10.3) 

FGFMAK* ERROR IN INTRP A PG . GT . X( N, 2 ) * ) 

FCFMAT{* ERROR IN INTPP AFG.LT.X(1, 1 )' ) 
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or>r>or>«-no M'tr'.ri /%r^o 


SICP 

ENC 


SL'EPGUTINE R AT E ( P ^ ^G E, SP E ED, iWPNt SU P ^AC, RO^ ) 

GIVEN THE RANGE ANC SPEED OF A TARGET ALONG WITH THE TYPE OF 
WEAPON BEING USEC TC FIRE UPON THE TARGET AND THE SUPPRESSION 
FACTOR THE FIRER IS BEING SUBJECTED TO, SUBROUTINE RATE COMPUTES 
THE PATE OF FIRE USEC AGAINST A PARTICULAR TARGET, 

COMMON /CEF/TENGMXtSENGMX ,SENGMN,TARTM,SAPTM,TVEL,SVEL 

COMMON /5UPEFT/GAMMA,DELTA 

RCF=0.0 

IF(PANGE.LT.25.) RETURN 
IF( IWFN.E0.2) GC TO 10 

IF(FANGE.GT.TENGMX) RETURN 

1PTM=TAPTM^ (1.0-^GAM^ A*SU?FAG) 

CT=TPTM4PANGE/(TVEL^SPEED) 

PCF=1.0/CT 

RETURN 

10 IF(RAMGE.GT.SENGMX) RETURN 

IF(RANGE.LT.S5NGMN') RETURN 

SRTM=SAFTMfc (1.0+ CAMM A^*'SUPFA C ) 

DT=SRTM + RANGE/(S VEL+SPEED) 

RCF=1.0/CT 

RETURN 

EMC 


IN THE FUNCTIONS HT,SPD, /ND PNG, THE ARGUMENT T 
IS the TIME SINCE THE WAVE BEING ADDRESSED 
CPCSSED THE 50C0 METER OFFSHORE MARK 

FUNCTION SPO(T) 

CCMMCN /ENGR/ SPCMAX,SPDM IN,HTMAX,HT MIN, TTS,TAA,TB,TFF 
If(T.GT.lAA) GO TG 5C 
SFO = S PCMAX 
RETURN 

5C IF(T.GT.Te) GC TC ICC 

SPO=S FCMIN+ ( (T6-T)/TTS )^(SFCMAX-S PCMlN) 

RETURN 

IOC' SFC = SP0M1N 
RETURN 
ENC 

FUNCTION HT( T) 

CCMMCr4 /ENGR/ SPCM AX , SPDM IN, HTMAX ,HT MI N, TTS , TA A , T3 ,TFF 
IF(T.GT.TAA) GC TC 50 
HT=HTNAX 
RETURN 

50 IF (T.GT.TE) GC TC 100 

HT=HTMIN+( (TB-D/TTS)* (HT M AX-HTMI N ) 

RETURN 

ICC H1=HTMIN 
RETURN 
ENC 

FUNCTION PNG(T) 

COMMON /AMPH/IL(5),we(2) ,A( 2) ,B (2) ,ITE,ISE,RD,WVINT( 5) ,WID, 
♦ Tev,DINn(2) ,GAINU,IWSTAT (5) 

COMMON /ENGR/ SPCMAX,SPDMIN,HTMAX,HTMIN,TTS,TAA , T8 , TFF 
IF (T.GT .TAA) GC TC 5C 
RMG = 5CC0. C-(SRCMAX'5'T) 

RETURN 

50 IF(T.GT.TE) GO TO 100 

PMG=PC-C. 5^(T-TAA)=^(SPCMAX + SPD(T) ) 

RETURN 

100 FNG=PD - ( ( (TB-TAA )/2.0)*(S FOM IN +SPDM A X ) )-( (T-TB )^SPOM IN ) 
IF(RNG.U7.75.) RNG=0.0 
RETURN 
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} 


ENC 

C 

C 

C 

c 

SieROUTUE GRCUNC(GATM,TSIRV »IPRI NT , ITS) 

C This IS THE PRIMARY SUBRCITINE CF THE LAND COMBAT PHASE OF 

C THE AMPHIBIOUS CiERATION. INFCPMATICN RECUIRED FOR THE OPERATION 

C CF THE LANC CCHBAT PHASE IS RE40 IN AND PRINTED IfJ A SUMMARY 

C TABLE FCF USER P-EVIE«, ThE INFCRMATICN PROVIDED BY ALL OTHER 

C SIBRCUTINES USED IN ’ME LAND CCMBAT PHASE ARE USED IN THIS 

C SLBRCUTINC AS INPUT TO THE BASIC LAND COMBAT ALGCP.ITHYM 

REAL*8 C5EED 

PEAL TSLrV,TTS fP (5 ) 

COMMON /GRPl/ IPRCIR(6 ) t I 'ECUC(6) tMVTDIR(6) ,X(6) ,Y(6),SP0(6 ) 
COMMCN /GRP2/ TA ( 2 ) t T :. ( 2 ) ,TH ( , TM( 2 ) , TF 1 ( 2) » TP2 ( 2 ) , TF3 ( 2 ) , 

♦ P (2 *6) ,PHH(2 t6)tFHM(2.6) f FKH(2»6),TF (2) 

CONMCN /GRP3/ NBUtNRUo FLI «) ,FC (o) ,NG I i 3i »XIC(5 t200J .YIC(3,200) , 

*ICIR(3,2CC),AVSr,ISPC 

*♦ IUSTAT(E),II(6),LaST(6,6),VISFRA,VlSFF8,SIZETK, 

♦ SIZETW,M(6) ,NF( 6),SRF,DMMAX , 

♦NLCSC(6,6) ,VISFR(6,6),RM LNTK,RMXTK, RMINTW,RMXTW,OP,TOWFR,LVAFR, 

*PTT(3,3) ,RF,PCA(6,6) ,AFOA f6,6 ) ,LOA(o ,6),NA(6),OFL(6),POL(6 ) 

COMMCN /GRP4/ TPCL(6) ,CLDC(6,6),Q(6, 6) 

COMMCN /GPF5/ LOT (6 ,6 > ,R G T ( 6,6 ) 

COMMON /HILLS/ >^C(1C0) , YC (10 C ) . PE AK ( 100) , SX (10 0 ) , 3Y ( 100 ) , .RHO ( 100 ) 
COMMCN /HILLS/ S C AL E (100 ) ,TV<C P HC ( 100 ), TWO SCL ( 1 00 ), B A SE 
COMMON /HILLS/ NFiLLS 

COMMCN /CCVE.R/ C / C ( 15 3 ) , C K ( 1 5 C) , CPE AK (1 5 0) , CP XX (150 ) , CPY Y ( 150 ) 

COMMCN /COVER/ C-XY(130),NCVELS 

COMMCN /CQLNTR/ K H , KH ^ , KV ,KN , HGRS , KE LL ,K I NT 

COMMON /GRID/ LSTdOtlO) ,^HL( iC,10) , LISTH<45 0) ,KHREP(100) ,KTREP 
COMMON /GRID/ LSTC(10,10 ) ,NC(10,10), L 1STC(400),KCREP(150) 

COMMCN /GRP6/ ALFHA(d 

COMMCN /CRP7/ XA(6 ),VA(6) ,lMCVEf6) 

C=**’*^ IMTIALI/ATI cn* 

C *** VARIABLE CEFINdICNE 

C BL - INITIAL CEJ'ENSIVE pCFCE LEVEL 

C RL - INITIAL AGGFESSCP FCFCE LEVEL 

C ZL * ZERO LEVEL 

ei=o,o 

PL=C.C 
MP = 0 

FAI=3.1415G 
ZL=.00CC1 


C 

c 

c 

c 

c 

C 

c 


c 

c 

c 


READ terrain DA*A FCR LIN I OF SIGHT 
CHECK FCR STOCHASTIC OR DETERMINISTIC ATTRITION 
ITRIT-A7TRITICN MODE 1=DETERMIN ISTIC 

0=S7CCHASITIC 
CSEEC-DOUSLE PRECISICN SEED NUMBER 

PP AND CC ARE The BETA D 1STR IBU T IG.N PARAMETERS FOR DEF UNITS 
PO AND CD APE THE BE7A DISTRIBUTICN PARAMETERS FOP ATK UNITS 

REAC(9,5CC) ITKIT ,CSEEC 
REA0(9,5C1) PPfCC 
CC 5 1 = 1,3 

CALL GG3TR ( CSEED ,PP, Q 0,1, R ) 

ALPHA (I) = P(1) 

CCMINUE 

READ(9,500) ITRIT,CSEEC 
READ(9,501} PC,CC 
DO e I =A,6 

CALL GGBTR ( CSEED,PD,QD,1,R» 

ALPHA(I)=P(1) 

CONTINUE 

READ IN NUMBER CF ATTACK AND DEFENSE UMTS 
NBU - NUMEER CF CEFE^SIVE UNITS 
NRU - NUMBER CF AGGRESSCR UNITS 
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c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


RE40(9,5C2) NBLtNnU 


INITIALI2E WEAPON SIZES 

SI2ETK - SIZE CP IVA WE4FCN SYSTEM 
SIZETW - SIZE CF TOW WEAPON SYSTEN 


SI2ETK=2.5 

SI2ETW=2*5 


RE^O IN EFFECTIVE WEAPON FANGES 

pyiNTK ANC RMXTK ARE MAX 4ND MIN RANGES OF LVA MOUNTED WEAPON 
RMNTW ANC RMXTW FOR MAX /NO MIN RANGES OF TOW DEFENSIVE WEAOON 


REAC(9»5C2J RM INTKtRMXTK,FMINTW,RMXT W 


INITIALIZE PM,RF ,TCWFR,LV/FR ANC NOD 

PM - FFCPORTICN CF TIME A MOVING UNIT IS SEARCHING pQR TARGETS 

RF - OETECTICN PATE RECLCTICN FACTOR FOR A FIRING UNIT 

(IN CCMFAPISCN to / NCNFIRING UNIT) 

TGWFR - FIRING R/TE CEFENDING TCW WEAPON SYSTEM 

IVAFR - FIRING RATE ATTACKING LVA WEAPON SYSTEM 

NOO - NUMBER OF TIME INTERVALS U^a T I DELAYED IK MOVEMENT 
(ICO FAR IN FRONT OF OTHER JMTS) 


PM=.352 
RF=.5 
T0WFR=.C3 
LVAFR=. I 
NOD = 2 

DO 10 I-ItNRU 
NOI ( I ) = 125 
10 CONTINUE 
K=NRU+I 
L=NRU+NEL 
00 15 I=3tL 
IK I ) = C 
15 CONTINUE 
C 

C*** READ IN FCRCE LEVELS CF EACH AGGRESSOR UNIT 
C 

ISURV= INT(TSURV/NRU1 
DC 20 1 = 1,NRU 

FL( I )=FLOAT(ISUPVi 
20 CCNTIMUE 
C 

CHECK FOR TYPE CF RCUTE OE TERM I NIT I ON 
C 

REA0(9,504) IRTE.ISPC 

VARIABLE DEFIMTICNS 

IRTE - DENOTES WHETHER USER WAN'^'S TO INPUT ROUTES OR NOT 
0 - FFOGRAM DETERMINED ROUTES 
I - LSER DETERMINED ROUTES 

INPUT VARIAELE TC DENOTE USER'S DESIRED SPEED FOR 
1 - 9 MPH 

2- 32 MPH 

3- 15 MPH 

4- 36 MPH 

AGGRESSOR FORCE NOVEMENTS 

AVERAGE SPEED CF AGGR5SSCP FORCE MOVEMENTS 
DISTANCE IN METERS TO 86 MCVED EACH TIME STEP BY 
AN AGGRESSCR UNIT 


ISPO - 


AVSC 

DST 


IF(ISPD.cG.l ) AVSP=S.O 
IFdSPD.EC.l ) CST = ^0.232 
IF(ISPO.EC.2) AV5P=12.0 
IF( ISPD.EC.2 ) CST = 53.6A3 
IF(ISPO.EC.3) AVSF=15.0 
IF(ISP0.EC.3) CST=67.053 
IF(ISF0.EC.4 ) AVSP=13.0 
IF(ISPO.EC.4) CST=80.463 


READ IN INITIAL ACGPESSOR LNIT'S LOCATICNS 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


30 

250 

255 


CO 25 1 = 1,NRU 

PEAC(c,506) XIC(I,1),YIC( 1,1) 
COMINUH 

IFdPTE.EG.l) GO TC 250 


DO 


30 

DC 


1=1 ,NPU 
3C J=2,125 

YIC<I,J)=YIC( I,J-1)^CST*(J- 1) 

XlCd ,J)=XIC< I,J-1)+CST*(J-1) 

IDI R(I,J)=0 
CONT INUE 
GO TO 255 
CALL ROUTE 
SL^PO=O.C 

STATE VARIABLE CEFINITICNS 
FLd) • FORCE LEVEL CF UNIT I 
SUMRC - TOTAL KGFESSOP FCRCE LEVEL 
MVTDIP(I) - ^^CVE^'E^T DIRECTIC^ CF UN I 
IDIR(I.J) - OIPECTICN CF THE JTH INTERVAL 
lUSTAT(I) - CUPPENT STATUS UF UNIT i 

0 - UNIT ALIVE A^C NOT FIRING 

1 - UNIT ALIVE ANC FIRING 

2 - UNIT DESJFGYEC 

3 - UNIT MCVUG 

NUMBER CF TIME INTERVALS UNIT I IS ALLOWED TC FIRE 
AT SAME LOCATION 
INTERVAL INDEX FOR UNIT I 


IN THE ITH ROUTE 


NF(I) - 
Ild) - 


35 


CC 25 1 = 1,NRU 
F0( I ) = FL( I ) 

SUMRC = SUMPG+FOd ) 

X( I) = )IC( I ,1) 

Y( I )=Y IC( 1,1) 
MVTOIR( I) =ICIRd,l) 
SPO( I )=AVSP 
lUSTAT (I )=0 
IPRDIRd) =ICIR<I ,1) 
ISECUD( I) = 120 
NF( I ) = 1 
IK I ) = 1 
CONTINUE 


READ 
C 

c 
c 
c 
c 


IN CEFENSIVE LNIT'S LCCATICN3 

STA^E VARIABLE CEFINITICNS '**'^^* 
IPPCIF(I) - PRINCIPLE CIRECTION OF FIRE FOR 
ISECWCd) - UICTH OF SEARCH SECTOR FCR UNIT 
SUM8C - TOTAL DEFENSIVE FCRCE LEVEL 


UNIT 

I 


AO 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SLMBa=0.0 
CC 40 I=K,L 

READ ( c ,50 7) X( 1 ) , Yd ) , f L ( I ) , IP PO I R d) dS ECW 0 ( I) 

F0( I ) = FL( I ) 

SUMBC=SUMeC+FO(I) 

MVTOIF(I)=C 
SPO( I ) = C.C 
lUSTAT(I)=0 
IMOVE d )=C 
CONTINUE 

CHECK FCR ALTERNATE DEFENSIVE POSITIONS AND READ IN IF WANTED 
STATE VARIABLE DEFINITIONS 

lALT - I^CICATES IF ALTEPrATE CEFENSIVE POSITIONS DESIRED 
0 - YES 
1 - NO 

BREAK - CLOSEST DISTANCE ALLOWED 8ETUEEN OPPOSING FORCES 
BEFORE breakpoint REACHED FOR PROGRAM 
ITEM - MMBER CF TIME INTERVALS ALLOWED FOR DEFENDER'S MOVE 
TC THE ALTERNATE DEFENSIVE POSITIONS 
TA(K) - TIME TO ACQUIRE A TGT FCR KTH WEAPON SYSTEM TYPE (K = 1,2) 

TdK) - TIME TO FIRE 1ST FOUND AFTER TGT ACQUIRED BY KTH WEAR SYS 

TH(K) - TIME TO FIRE A RCLND FOLLOWING A HIT FOR KTH WEAR SYS 

TN(K) - TIME TO FIRE A RQLND FCLLQWING A MISS FOR KTH WEAR SYS 
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TFl(K) - TIME OF- FLIGHT FCR KTH WEAR SYS PROJECTILE TO 1000 METERS 

TF2(K) TIME OF FLIGHT FCR KTH WEAP SYS PROJECTILE TO 2000 METERS 

TF2(K) - TIME OF FLIGHT FCR KTH WEAR SYS PROJECTILE TO 3000 METERS 

RE^O(9t5C£l lALT ,EREAKtlTEM 
IFdALT.EC.l ) GO TC 260 
00 45 I=K»L 

REAC(9,506) XAdltY^d) 

45 CONTINUE 
26C CELT=^10. 

JA (l)=2C, 

Tl (i ) = 3. 

TH(1)=8* 

TM(1 )=10 . 

TFKl) =1 . 

TF2<1) = 1. 

TF3(1)=!• 

TA(2)=20. 

Tl(2) = i3. 

TH(2) = 8, 

TM(2)=15. 

T Fl(2) = 10. 

TF2(2)=12. 

TF3(2) = lf* 

REiSC IN HIT AND KILL FFCBA EILIT lES 

STATE VARIABLE DEFINITIONS 

P(I.J) - FROB 1ST RC'JNC HIT BY UNIT I IN RANGE BAND J 
FHH(I.J) - PRCe OF HIT FCLLCaED BY A HIT 

PFMdtJ) - PRCe OF HIT FOlLQWcD BY A HISS 

PKHdfJ) - PRCE OF A KILL GIVEit A HIT 

PTTdfJ) - PRCPOPTICN SURVIVING FIRE PCW5R ALLOCATED TO 

ITH TARGET IF J TARGETS ARE AVAILABLE 
NLCSCdfJ) - NUMBER GF CCMINUCUS TIM? INTERVALS THAT A LINE OF 
SIGHT(LGS) DC5S NOT EXIST BETWEEN UNIT I AND UNIT J 
OdtJ) - FPOEAEILITY UNIT J NCT DETECTED BY UNIT I AT CURRENT TIME 
VISFRdtJ) - FPACTICN CF FEIGHT OF T GT J VISIBILE TQ rlRER I 
IRAN - RANGE 

DO 55 1=1,2 
DC 5C J=l,6 

RE AC (9,5 09 ) P ( I, J ) , FHH ( I, J ) , PHM( I, J ), PKH (I , J ) 

50 CONTINUE 
55 CONTINUE 

PTT(1,1 ) =1.0 
PTT(1,2)=C.B 
PTT(2,2)=C.2 
PTTd,3) =C.8 
PTT(2.3)=C.15 
FTT(3,3)=0.05 
DO 6C I = 1,NRL 
DO 6C J=K,L 

NLCSC( I , J)=0 
NLCSC(J,I)=C 
0( I, J)=1.0 
Q(v,I)=1.0 
VdFR( I , J)=C.O 
V ISFR(J, I)=C.O 
6C CONTINUE 
IC = 1 

C*** PRINT INITIAL BATTLE INFORNATICN 
C 

WRITE(6,6C0) 

WRITE(6,6(jl) 

CO 65 1=1,L 

WRI TE (C,603 ) I ,Xd ) ,Y ( 1) , FL (I ) 

65 CONTINUE 

IFdTRIT.EO.l) GC TO 265 
WRITE(6,6G4) 

GO TC 270 
265 WFITE(6.605) 
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27C IF (IRTP.EC.O ) GG TC 275 
WRIT£<d*6C6) 

275 V^PITE(6,6C7) ^iVSP 
VsPITE( 6,6C€) BFE^iK 
IF( lALT.EC.O) GC TO 28C 
WRITE (6t609 ) 

GO TO 2F5 
280 WPITE(6,(eiC) 

CC 70 I = t-,L 

WRITE <6 ,6 11 \ I ,XA(n r V^( n 
70 CONTINUE 
265 IPAN=50C 

WPITE(6,612) 

CO 75 1=1,6 

WRITE <6,6 13) IRAN ,P (1 , 1) , PHH(1 ,1) ,FHM(l,I),PKH(l,n 
IRAM=IPAN+500 
75 CONTINUE 
IRAN=5C0 
WF ITE(6,61'») 

00 80 1=1,6 

WRITE <6 ,612 ) IRAN ,P(2 , 1) ,FHH(2,I ) , FHM <2,I),PKH(2 , n 
IRAN= IRilN + 5C0 
8C CONTINUE 

WRITE< 6,61 = ) 

WFITE<6,617) 

WFITE< 6,61E) 

C'M* UPDATE LCCATICN CF REC UNITS. 

NF - NUf^BER OF TIVE INTERVALS U^a T I COES NOT FIRE 
NOD - NUMBER OF TIME INTERVALS UNIT I OELAYEC IN MOVEMENT 

OI5MAX=5CCC.0 
2S0 00 SO 1=1 ,NRU 

IF( lUSTAT< I ).EQ.2 ) GO TO SO 

IFdUSTATd ).E0.0) GC TC 295 
NF( I )=NFd) + l 

IF <NFd ) .LT.NCO) GC TO 90 
NF d ) =1 

295 00 85 J = 1, NPU 

IF <J .EO. I) GO TO 85 

IF (IUSTAT<J) .EO. 2) GO TO 85 
GIST = Xd ) - X< J ) 

IF (CIST .GT. OISMAX) GO TO 90 

85 CONTINUE 

IK I ) = I I( I) 4 1 
K7=IIdJ 
X<I ) = >IC( I ,K7) 

Y( I )=>101 I,K7 ) 

MVTO IF(I) =IClF d,K7) 

IPRDIF(I)=ICIR(I,K7) 

SO CONTINUE 

line—OF-SIGHT CHECK BETWEEN UNI^^S AND TARGETS SELECTION 
STATE VARIABL E OEF INITIONS 
NTd) - NUMBER OF TARGETS DETECTED BY UNIT I 
XXI,YYl - COORCKATES CF UNIT I LCCATICN 
yX2,YY2 - CJCRCINATES CF UNIT J LCCATICN 
TMACI,TMACJ - ELEVATICN OF UNIT I ANC UNIT J 
0,0 - INCICATES NO UNITS UNDER GROUND 

SI2ETK, SIZETW - SIZE GF L\^A VEHICLE ANC SIZE OF TOW VEHICLE 
LATCB - INDICATOR VARIABLE FCR CME OR TWO WAY LOS CALLS 
0 - CC NOT COMPUTE LOS FROM UNIT A TQ UNIT B 
1 - COMPUTE LCS FFCM UNIT A TC UNIT B 
(YIELDING VISFPBJ 

VISFRA - FPACTICN CF HEIGHT CF TGT AS SEEN BY UNIT A 
LC£T(I,J) - INDICATES IF LCS EXISTS eETWcEN UNIT I AND UNIT J 
0 - NO LCS EXISTS 1 - LG S EXISTS 

CO 95 J=K,L 
NT(J )=C 
S5 CONTINUE 

CO 105 1=1,NRU 
NT(I )=C 
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n;ST/iT (I ).5C.2 ) GO TO 105 
OC 100 J=K,1. 

IF( ]UST4T(J).EQ. 2.0F.IUSTAT(J). EQ.3) GO TO 100 
>X1=X (I ) 

>Y 1 = Y < I) 

CALL ELEV(XX1 ,YY1,TMACI) 

XX2=X(J) 

YY2=Y( J) 

CALL ELEV(XX2 iYY2fTMACJ) 

LATOe=l 
L B7CA = 1 

CALL LOS (XXI •'lYl t THAC 1,0 ,0 ,S I 26TK , X X2 , Y Y 2, TMAC J , 0 . C , 
LA7CB,L0TOA , V ISFRA ,VISFF3» 

VISFR(I,J)=VISFRA 
\ISFR(J,I)=VISFRe 
IF ( VISFRA^GT. 2L) GOTO 30 0 
LOST( I, J) = 0 
L05T(J,I)=0 

NLOSC (I, j) =i\'L0 5C4 I ,j; ^'l 
NLOSC (J , I) =NLOSC( I ,J ) 

GO TO ICO 

300 10ST(I,JJ=1 

l.OSK J,I ) = 1 ‘ 

NLOSCd ,J) =0 
MLOSC(J, I)=0 

RANGE = $aRT((X (I )-X(J) )^^2 + (Y ( I )-Y( J) ) 

)F(RANGE.LT.RMINTK.GR*RANGE.GT.RMXTK) GC to 305 
IF (0( I, J) . cQ.l .0 ) GO TO 305 
lUSTAK I) =1 
NT( I )=N T( I J + 1 
H=NT( I J 
LOT(I,M ) = J 
RCT( I,M )=RANGE 

I F( M.EG.l ) GC TO 20 5 
CALL £GRT(I,M) 

305 ..F (RANGE.LT.RMirJTW.OR.RANGE.GT.RMXTW) GO TO 100 

IF(0( vM). EG.l.C) GO TO 100 
IUSTAT(J)=1 
NT(J) = N1(J)41 
M = NT(J ) 

LCT( J,^^ )=I 
RCT( J,M )=RANGE 
IF(M.EC.l I GO TC 1 CO 
CALL SORT(J,M) 


ICC CCNTI^UE 
K5 CCNTI\'LE 


CO 110 


a »NRU 


IF( lUdATI I ).EC.2) GO TO 110 
IF(NT( .1 ).NE.C) GC TO IIC 
IUSTAT(I)=0 
NF( [ ) =0 

no CCNTHiUE 

CO 115 J = K,L 

IF( IU<TAT{J).cC.2.0R.USTAT(J).EQ.3) GC TO 115 
IF(NT(J ) .EC.O) lUSTAT( J)=0 
115 CCNTINIJE 

UPCATE CF TFE ACCUMULATED CETECTION PROBABILITIES. 

C 

lAA^'l 
IEE^NRU 
ICC = K 
ICC=L 
FF=:TGWFF 
CF = PM 

CO 120 1=1,6 

12C CONTINUE 
3C7 DO 135 I=IAA,ieB 

IF( lUSTAT ( I ).EC.2.QR. ILSTAKD.EQ .3) GO TO 135 
DO 120 J = ICC,ICD 
FRO P = C. C 

IF( lijST/iT(J).E0.2.0R.IUSTAT (J). £0.3) GO TO 130 
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21C 

315 


3 50 
315 

»0ELT«=L(J 
TO 31G 
V 

=^DEIT*FL( 

G?-CV*C?) 


320 

325 

330 
•: -a c 


12f 


34C 

345 

35C 


130 CCNTINU3 
135 CCMINJE 

iF( IX^.EC.O GC TC 3: 
rR=L / 

l4i: = K 

IB3=L 
ICC=1 
IDD = Nt:U 
CP=1. c 
GC TG 3C7 
C 

ALLCC^TICN. 

C 


CLCC{I,J}-C(I.j) 
lF(LC51(i«J)*tG«C) GC TO 
IF(M (I J.GT.O ) GO TG 
FCTvI5 = MSFr(J» I) 

CALL L A'-OA { I ,J.PCT 
CV=EX?{-(CETFAT^OP 
IF (N'T( J /. GT .C ) GO 
Qd , V)=C ( 

GC TG 13C 

CP = ( 1. C-PSL5K)-^*{F 
C{ I, J) = C( 1 + 

GO TO 130 
N5=NT(I) 

CC 12 5 

K1=LGT( !,II > 

ANGi=AT;N2(^(<I)-T 
A^GZ = AT^^d( > {J)-Y{ 

IFU ANGI^ANC-ZJ .Gz. 

I F{ ANG2, L’.O.O) 
a>'g=2’‘pa: +A.\ 

GO TG 325 
ANG= :*A,\’G2- 
iF(A^G.GT,PAI} 

GO TQ 33 5 
AN’G^^ABS (A^^Gl-ANGl ) 
AA=15.0 «pa:/130.0 
!F( ANG. CT. AA ) GO 
5RC? -PRO - -"^TTC J 
CCNTI VJ3 

I^■{FKL^#cC a^« ) \S\J • C’ 

IF{NT(J :.GT .0 )GCTC 
CALL LA.-CiC I .J, 

5v=E>Px'-{.-AC?^ 

0(1, J) = :< :,j) = 

GC TC 13C 

c (I, J) = 0 .: 

GC TC 13C 

IF( IAA.3C.: J GC 130 
0(1.Ji=1.C 
GC TO 130 

IF(\LCSC( I,J ) GC TO 130 

ClI Vw) '1«0 


,PSUEK) 
} } } 


J) ) 


( ! ) , X (<1 ) - 
I ) ,X( J)- <( 
J.O) GC TO 
GC -C 320 
G1-ANG2 

ANSI 

A\C*=2*=>AI^ 


G 12 5 
1,N5 ) 

345 

3^*0 

fCTVlStDE" 

i^FAT^^OELT 

V 


X ( I) ) 
I ) ) 
330 


A.NG 


PATtPSU 
^-L(J)) 


state VAP:A3Lz OE^IM’ICN 
APCA(I,o? - AVEPaGE -'CPC=’:CN C?= Trt 

allccate: TC fire :n’ uni 


5 *MMXX 

c jT H AGGP: 
T I 


355 CO 1^0 I=ItL 


140 


145 


>A{ I )=C 


DO 


1^5 

IF( 


1=1 fL 

:letat(: j. ec. 2 .of. iusta- 


IF(NT(I ) .5C.0 ; GC TO 155 
CC 145 J = I ,5 

APCAC, .)=C.C 
CCMINJE 

IF(NT(I ). EC.l ) GO TO 
IF(NT( I ), 53.2) GO 
NXT=3 

v>'i = icT(: ,1) 

MY2= LOT{I ,2) 


(IJ* cC.3 J GC TG 1 


370 

TC 365 


EK) 

) 


r UNIT : 


55 
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36C 


150 

365 


37C 


MM3=L0T(1,3) 

PF0E-{1.C-C( IfiM.vi ) )^C(I,MM2 I, MM3) 

APOA (I , 1 ) =APaA( I ,1) i-PTT( 1 ,1 )*ocr0 
PROe-Q( I ,MM1)4=( 1.0-C( I,MM2) )-^Q( i,MM3) 

APOA (I ,2 ) =APCA( I,2)-^PTT(1 ,1 I^PPCE 
PROS =Q(I )KQ(I,MM2)*(1.0-C(I,MM2) ) 

AFC A (1,3 )-APCA( 1,3 ) ^PTT (1, i ) ^pRC^B 
PROB -( 1. C-C( I, ) » -« (1.0-C< )^Q{ I ,MM3 ) 

AFOAd, l)rAPCA( I, 1) -►PTK 1,2 )>?Ff.B 
AFOA (I ,2) = APCA( 1,2) +PTT(2,2 )^PRr:3 
PROB = ( 1 . C-Q( I ,M M ) ) I ,MM2 ) * (1.0-Q( ) , MM3 ) ) 
APOA (I , 1 ) = APOA( 1,1) +PTT( 1,2 ) '^PRQB 
APOA (I ,3 ) =A?CA( I ,3) -^P'T (2 ,2 ? -C2 

PROB-Q( I ,MM1 )*( 1. C-C( I ,MM2) )^ ( 1.0-0 ( ) ,MM3 ) ) 
AFOA (I,2)=APCA( 1,2 )+PTT f 1, 2 » OB 
APOA (I ,3 ) =A?CA{ I,3 )-k?TT( 2 ,2 l-^PRCB 
FFCB= (1.0-C( I , MMl) )i^( l.C-0( I ,MM2 ) )-t( i.c-0{ ) ,MM3) ) 
APOA (I , 1 )=APCA( 1,1 ) 4PTT(1 ,3 )^P^^C3 
APOA (1,2) =APCA( ! ,2) +PTT( 2 ,3 1 
APOA (1,3 ) = APGA( 1,3)+?TT{3,3)^^RQB 
00 150 J=1,NCT 
— KK=LOT( I , J) 

NA(K X) = NA(KK)+l 
IN = NA(KK ) 

LCACHK, IAJ)=I 
PCA(t‘K,IN) = APOA{I,J ) 

CONTINUE 
GO TO 155 
NCT = 2 

MM1 = L0T{I , ]) 

MM2=L0T(1,2) 

PFCB= (1 .0-C( I ,N^'l) )*Q (I ,MM2) 

APCA( I, 1) = /PO A (I, 1 ) -HP TT( 1, D^PR OB 
PFCB=C(I ,MM)^( 1.0-Q( I,MM2) ) 

AFCA(I ,2) =/PCA(1,2)+PTT(1,1) *PRC3 
PFCE= (1.0-C( I ,MM1) )A( l.C-0( I ,xm:) ) 

AFCA ( I , 1 ) = / PCA ( 1,1 )-hPTT (1 ,2 )-PR 03 
APCA( 1,2) = APOA ( I, 2 )«*-P TT(2 ,2)^PR03 
GC TC 360 
NC1 = 1 

MM 1 = L0T(I, 1) 

PRCE=1.0-G( I, ^M1 ) 

APCAd ,1)=AP0A( I ,1)+PTT( 1 ,l)^PROe 
GO TC 360 


155 CCNTINUE 
C 

ATTFITICM CCMFUTATICN 

STATE VARIABLE OEFIMTIGNS 

RANGE - CURRENT MINIMUM DISTANCE BETVsEEN AGGP-ESSIR /NO DEFENDER 
FOA - PPCFCRTICN CP THE JTH ATTACKER OF UNIT I ALLOCATED TO 
FIFE ON LNn I 

TFQL - TOTAL PERCENTAGE LOST SINCE START CF BATTLE ^CR UNIT I 
AVD - AVERAGE CISTANCE 


375 

3dO 

335 


SUMR=0.0 
5UMB=0,C 
CO 165 1=1,L 

IF( IU5TAT(I).EC.2.OR, ILSTAT (1) .EQ .3 ) GC TO 165 
M6 = NA( I) 

SUM=C .C 

IF(M6. EC.O) GC TO 3^5 
CO 16C J = 1 ,M6 
M7=LQA(I,J) 

IF(M7.LT.K) GC TC 375 
nTPE = 2 
GC TO 380 
nYPE = l 

RANGE = SQR'^( ( X (I ) - > (M7 ) ) * *2+( Y ( I )-Y ( M7 ) )**2) 
IFdTRIT.EO.l ) GC TO 385 

CALL STCCH (ITYFc,RANGE,AJI ) 

GO TO 350 

CALL ETKdTYP E,RANGE,T ) 
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AJI = UO/T 

SUM=SUM+AJlAFl(M7)’fPGA( I , J)*OELT 
CCNT INUE 
QFL(I) =Fl(n 
FL(I J=FL(I J-SUM 
IF (FL( I )*GT.2L ) GO *10 AOO 
FL(I)=C.O 
IUSTAT( I ) = 2 
IF( I.LT.K) go to 405 
SLMB=SUMB + FL(I) 

TPQL( I )= (FO (I )-FL( I) ) /FO(I ) 

GO TO 165 
SL^P=SUMR+FL(I ) 

TPCLd )=(FO(I)-FL( I ))/FO(I) 

165 GQNTINUE 


290 

160 
-C c 


ACC 


405 


PRINT AND CHECK FCP BATTLE TERMINATION. 

C 

c 

ITIME=ICMNT (TTS) 

DC 175 I=KtL 

1F( IU5TAT(I ).EC.2} GO IQ 175 
DO 170 J=ltNPU 

IF( lUSTAK J). EG. 2) GO TO 170 
CHECK = X( I )-X( J) 

AVC = SGRT((X(I )-X(J))*^2+(Y(I )-Y(J) ) 

IF(AVD.LT.6RE4K.OR.CHECK.lt. 50.) GO TC AlO 
170 CCKTINUE 

175 CONTINUE 
GO TO Ai5 
C 

C*^* COMPLETE 4GGRE35CR UNITES ^QVE 
C 

AlO CO 180 I=K,L 

IF( IAIT.EC.1.0R.I'-CVE( I). EC. ITEM) GO TO AAO 
IF( IUST4T(I ).EC.O) lUSTATd )-3 
IMOVEd ) = IMCV£( I ) + l 
IF( IMOVE(I).LT.I lEM) GO TO 180 
X( I ) = X4( I) 

Y( I ) =Y4 d ) 

IF ( IUSTAT( I ) . FQ. 2 ) lUSTA Td) =0 
180 CONTINUE 

415 IITIME = ITIME+IFIX(GATM) 

IF( IPR INT.EC. 1) GO TO 430 
WR ITE(6,622 ) IITI 
^RITE(6,622) 

WRITE(6,624) 

DO 185 1=1,NRU 
M6=NT(I) 

IF(N6 .NE.O ) GO TO 420 

WRITE (6 ,62 5 ) I, X( I ), Y( I ) , FL ( I ) , lUSTAT (I } , TPQL ( i) 
GO TO U5 

420 VRITE(6,6 26)I,X(I),Yd),FL(I), IUSTAT( I ),TPOL( I ) , 

*(LOT(I,J ),J=1,N6) 

185 CCNTINUE 

V<PITE(6 ,627) 

WRITE(6,628) 

NNN=NRU+1 
DO ISC I=NNN,L 
N6 = NT( I) 

IF(N6.NE.C) GO TO 42 5 

WRITE! 6 ,62 9 ) I, Xd ) , Y( I ) , FL( I ) , lUSTAK ! ) ,TPOL( II 
GO TO 190 

42 5 WRITE(6,6 30 ) I ,X(I ), Y (I ), FLd ), IUSTAT( I ) , TPOL( I ) t 

♦ (LOT(I ,J ) ,J= 1,N6) 

190 CCNTINUE 

43C CONTINUE 
C 

C*** CHECK FCR EATTLE TERMINATICN. 

C 

ICT=0 

C*** CHECK IF 4N AGGRESSOR FORCE UNIT IS STILL ALIVE 


121 





00 195 1=1 ,NPU 

1F(FL( I ).EO.C.C) GO TO 195 
I0T*1 
195 CONTINUE 

IF( IOT.50.1) GO TO 935 

C*** NO AGGRESSOR FORCE UNIT ALIVE SC DEFENDER WINS 
WRITE(6,621) 

GO TO 950 
935 ICT=0 

C*** CHECK IF A DEFENSIVE FORCE UNIT IS STILL ALIVE- 
CC 200 I=K,L 

IF(FL (I I.EC.O.Oi GO TC 200 
IOt= 1 
200 CONTINUE 

IFdOT.ec.l) GO TO 995 

C*** NO DEFENSIVE FCRCE UNIT ALIVE SO AGGRESSOR WINS 
WRITE (6,632 ) 

GO TO 950 

C*** CIS'^ANCE EETWEEN FORCES TOO CLOSE 
99C WRITE (6,632) 

GO TC 950 
995 IC=IC+1 
GO TO 290 
95C continue 

WRITE(6,622) IITIME 
WR nE(6,623) 

WR ITE( 6,629) 

CO 205 1 = 1,NPU 
N6 = NT (I ) 

IF(N6.NE.C) GO TO 955 

WR ITE(6,625 ) I ,X (I ) ,Y( I ),FL(I) , lUSTAT(I),TPOL( I ) 

GO TC 205 

955 WRITE (6,626 )I,X(I ),Y( I ),FL (I),IUSTAT(I),TPOL(I) , 

♦(LOT(I,J),J=1,N6) 

2C5 continue 

WRITE(6,627) 

WRnE(6,62£) 

NNN=NRU+1 
DC 210 I=NNN,L 
N6=NT (I ) 

]F(N6.Ne.C) GC TC 960 

WRIT5(6,629 ) I,X(I ) ,Y(I ) ,FL (I) , IUSTAT(I),TPOL( I ) 

GO TO 21C 

960 WRITE (6,620 ) I,X( 1 ), Y( I ), FL ( I) , IUSTAT( I ) ,TPOL( I ) , 

*(LOT(I,J) ,J=1,N6) 

210 CONTINUE 

500 FORMAT(18X,I2,12X,F1C.2) 

501 FORMAT(35>,F12.2,7X,F12.2) 

502 FORMAT(16X,I 2 , leX, 12 ) 

503 FORMAT(10;<,F7.1,10X,F7.1 , lOX , F7.1,10 X, F7.1 ) 

5C9 FORMAT!If>,I2, ISX,12) 

506 FORMAT(12X,Fa.1,12X,F3.1) 

5C7 FORMAT(12X ,F8.1 ,12X , F3.1,6X,F9.1,7X, 13,7X,13) 

508 FORMAT(19X,I 1,12X,F8.1,29 X,I 1) 

509 FOFHAT(5X,F5.2,5X,F5.2,5X,F5.2,5X,F5.2) 

600 FORMAT! • 1*,IX,INITIAL LA ^ 0 C0M3AT INFORMATION «* ’, 

*////,' LCCATICN') 

601 FOfMAT!' UNIT',7X ,'X' , £X ,'Y SAX, • FCRCE LEVEL') 

602 FORMAT!1X,13,3X,FT.1,2X,FT. 1,7X,F5. 1 I 
609 FGRMAT(/1X,'ATTRITION IS STflC HAST IC ' / ) 

6C5 FCRMAT(/IX,'ATTRITICN IS ££T£FMN1ST IC'/ ) 

606 FORMAT! 1> ,'ROUTES DE'EPMUEO 6V USER'/) 

607 format(1>,'ATTACK VEHICLE SPEED IS ',F9.1,' M.P.H.',/) 

6Ce FORMAT!IX,'BREAKRCINT DISTANCE IS ',F6.i,' METERS',/) 

609 FORMAT!IX,'DEFENDER WILL FOT MOVE TO ALTERNATE PCSITICNS'/) 

610 format(IX,'DEFENDER WILL HOVE TO ALTERNATE POSITIONS' /IX, 

♦ •ALTERNATE POSITICNS ARE: '/1X , •UN IT• ,5X,'X',3X,'Y ' ) 

611 format(IX,13,3X,F7.1,2X,F T.1) 

612 F0FMAT!/9X,'ATK KILL PROB/61LITIES'/IX,'RANGE•,9X, •P', 

♦ 9X ,'PHH • ,2X, 'PHM ' ,2X,'PKH ') 

613 F0PMAT!2X,I9,9!2X,F9.2)) 

619 FORMAT!/9X,'DEF. KILL FPC EABILITIES•/IX, •RANGE',9X, • P •, 
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I 







, 'PHH • ,2X. 'PHf'* ,3X,*PKH •) 

£15 FOPFAT(//,' CURRENT STATU' OF UNIT I VARIABLE DEFIN ITI CMS• i//i 
*• C - ALIVE MOT FIRING*,/,* 1 - ALIVE AND FIRING*,/, 

♦ * 2 - KILLED*,/,* 3 - FCVING*) 

617 FOfMATI//,' VEFICLE SPEED VARIABLE DEFINITIONS*,//, 

♦ • 1 - 9 f'PH*,/,* 2-12 .‘PH*,/,* 3-15 MPH',/,* 4-18 MPH*) 

618 FORMAT!* 1 *,10X,• ***** THE LAND COMBAT PHASE BEGINS *****•,//) 

622 FORMAT (//IX, •TI.vE=' , 14,IX ,*SECCNOS*/ /) 

£2? FOF.'1AT(/,' AGGPESCR LNIT INF C R NAT ION ' ) 

624 FORMAT!1>,*UMIT',5X,'X*,e>,*Y* ,5X.'FCRCS LEVEL* ,2X,•STATUS* , 

*2LGST-FCT*,2X,'TARGETS•) 

62 5 FORMAT! 2>,Il,i>,F7.1,2X,FI. 1,6X,F5. 1,9X, I1,6X,F5.3I 

626 FORMAT !3>,I1,3X,F7.1,2X,F-;.1,6X,F5. 1,9X, I1,6X,F5.3,3>,3(I1,1X)) 

627 FORMAT!/,* DEFENSIVE UNIT INFCPMAT10N* ) 

628 FOPMAT!1>,*UMIT*,5X,•X *,g>,*Y*,5X,' FCRCE LEVEL* .2X,*STATUS* , 

*2X ,* LOST-FCT *,2X,* TARGETS *) 

629 FOPMAT!3>,II,3X,F7.1,2X,F7.1,6X,r3. 1 ,9X,11,6X,F5.3) 

630 FOPMAT!3>,!1,3X,F7.1,2X, = 7.1,6X,F:. 1,?X,I1,6X,F5.3,3X,3!I1,1X)) 

631 FORMAT! 1 >, *****.i AGGRESCP FORCE IS ELIMINATED. END OF BATTLE.') 

632 FORMAT!!),* **** DEFENSIVE FC^CE IS ELIMINATED. END OF BATTLE.*) 

633 FORMAT ! IX, ****** DISTANCE SETV*EEN FORCES IS TOO CLOSE. *, 

** END OF EATTLE ****** ) 

RETURN 

END 

SUBROUTINE SETLP 

♦*■* SUBROUTINE SETUP IS USED TO READ IN THE TERRAIN DATA AMO 
CREATE PARAMETRIC TERRAIN. THIS TERRAIN DATA ,<ILL Sc USED 
WHEN COMPUTING LINE-CF-SIGHT EETwEEN TARGETS AND OBSERVERS 
AS UELL AS PFCVIDING A GRID SYSTEM FCR UNIT LOCATIONS AND 
ROVtiMENT. 

CCMMCN /HILLS/ XC !100),YC (100 ),PEAK!100),ANGH!100),SPRD!100) 
CONf^ON /HILLS/ EC C ! 1 CO) , P )X! 1 C C) , PYY ! 1 CO) , PX Y! 100 ) , S AS E 
COW'CN /HILLS/ NFILLS 

COMMON /COVER/ CXC!150),CYC!150),CPEAK!150),CP XX!150),CPYY(150) 
COMMON /COVER/ CPXY ! 15C) ,RCV £LS 
C 0.VMCN / COUN TP /KH, K KV., KV , KN , K C-RS, KEL L , KI NT 

COMMON /GRID/ LST!5,4) ,NHL! 5,4),L1STH!150),KHR5P!150),KTREP 
COMMON /GRID/ LSTC!5,4),NC(5,4),LI STC!400),KCREP!150) 

FA I==5. 14 159 
R£AP,!5,5COI NHILLS 
REA(;!5,£1C) BASE 
CO 50 1 = 1 ,NHILLS 

READ! 5,52 0 XC! I ) ,YC! I ) ,PE AK ! I ) , A NGH! I ) ,SPRD! I) , ECO ! I) 

50 CORtINUE 

PEA0(5,£3C)LST 
REA()(5, 53C)NHL 
FEA(:!5,5CC)NHTCT 

READ(5, 530! LI STH! I ) ,I =1 ,MHTCT) 

65 CO .100 1 = 1,NHILLS 

AAMGLE = ANGH!I )*FAI/130 . 

SAMG='IN!AANGLE) 

CANG=COS!AANGLE) 

A=PEAI‘!I) /! FEAK!I )-50, ) 

A=ALCG!A) 
e=A* ECC !I )X‘*2 
SSPD= 5 FRO!I )**2 

PXX! I )=-!A*CAMG*CANG+B*SANG*SAMG) /SSPO 
PYY!I )=-(A*SANG*SANG + B XCANG*CANG) /SSP 0 
PXY! I ) = !2.=f'£ANG*CANG*< E-A ) )/SSPD 
KHREF( I )=-2 147483600 „ 

C ALL VALUES MOV, IN METERS CN Q — 10,000 GRID 
ICC CONTINUE 

REAC!5,5CO) NOVELS 
IF!NCVEL'.EQ.O ) GC TC 200 
00 15C 1=1,NOVELS 

REAC!5,54C)CXC!I),CYC! I ),CPEAK!1),CPXX!I ),CPYY!I),CPXY!I) 
KCFEF!I)=-21474836CC 
150 CONTINUE 

READ !5,530)LSTC 
REA0!5,53ONC 
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READ ( 5,50C?NC TCT 

REAC(5,53C) (LISTCd ),I=I,NCTOT ) 

2CC K iP6P=-2l474S3cC0 
KH=0 
KHW^O 
KV=C 
KN^O 
KGPS=0 
KeiL=0 
K IM = 0 

5CC FCPiMAK 20,14) 

510 FOPf^AT( ie>,r5. 1) 

520 FGP,HAT(12X,F7. i , 2 X ,F7. 1 , 5 >, F6 . 1,5X, F 6 . 2, 5X , F8. 2 , 4X , F4 . 1) 

520 FOPMAT( If >aci 5) 

540 FCRVAT(2 f1C.4,3E12.7 5 
RETURN 
ENC 

SUcRGuTI^E RCUTE 

*** SUEROUTINE RGUTS CCMPUTES THE ROUTE OF EACH AGGRESSOR UNIT 
WHEN THE USER HAS SEIECTEC THE OPTION CF INPUTING AGGRESSOR 
ROUTES. IT CALCULATES THE CQCRCINATES CF EACH INTERVAL EMOPCINT 
ALCNG the ROUTE,faking EACH IMERVAL L EN GTH ( 01 S T AN C E MOVED DURING 
A 1C SECCND TIME THE SAFE. THE INTERVAL LENGTH IS DETERMINED 

eV THE SPEED THE USER HAS SELECTED AND INPUTED FOR THE CURRENT 
BATTLE. 

CCfMCM /GPP3 / NcU .NRU. FL ( S ! , FO (6 ) ,N0 I (3 ) ,X ICO ,200) , YIC( 3, 200) , 

»ICIR(5,2CC),AV SP,ISPC 

ILSHAT(6 ) ,I K c),LCST( c,6 ),VISFRA,VI SFRB,SIZETK, 

♦ SI2ETW,NT(6) .NF(b ; .SRF ,01 S"ax , 

*NLCSC( 6 ,A ) .VI SF- ( 6 r6 ) .RM! .< ,RFXTK,RMI fTT w , RMXT W , C P , T OW FR , L V AFR , 
*PTT(3,3) ,RF,PCA(t ,c) PCA (1,6 ) ,l0A(6,6),NA(6),OFL(6) ,?0L(6) 
CIFENSION XLCCO ,?C5) ,VLCC (2 ,2C),N(3) 

IF( ISPO. EC.4 ) CST = Fr..^t2 
IF(IS?0.£C.3 ) !:ST=67,053 
IFdSPD.EC.a ) CS^ = 53.c43 
IF( ISPD.EC.l ) CST=.4C. 132 
DC 110 1=1,NRU 

REAO(S,50C) N(I) 

NL = N( I )-l 
OC 1C lN=2aL 

R£AD(9,5]C) xlg:s,ylc:cs 

XLCCd ,IN)=XLCCS 
YLCC(I A)=T LCCS 
CON T I ME 

XLCC(I,1)=A1C( 1-1) 

YLQC( I ,1)=VIC d ,1 ) 

IOIR( 1 , 1) =C 
NL=N ( I ) 

NUM=2 

CO ICC J=1,NL 

XL = /LOCd,J-l)-<LaC (:,J) 

YL = YLOCI I , J+ 1) - MCC (1 . J ) 

0 IST=SCRTl XL>-^2-YL^^2) 

Y = A£S( YL ) 

Z=i/XL 

ANCL=ArAN(7 ) 

DEC=ANGL»5T.2958 
IF(J.EQ. 1) GC tq 50 

XLN=(CST-PXTRA)-C0S(ANGL ) 
cist=(oist+e:<tra )-ost 
YLN=(CST-EXTRA)* dN( ANGL ) 

XIC(I,NUF) = XIC(I ,NUM-1)+XLN +XLE 
IF( YL .GT.C.) GO 10 2C 
YLN=-YLN 

YIC( I ,NUM)=YIC(I ,NUM-1)+YLN+YLE 
IF( YL.GT.C.) GO TO 3C 

ID IP(1,NUM) = - ]FIX(CEG) 

GO TO 40 

IOIR(I,NLM)=IFIX(DEG) 


20 

20 
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I 


40 

NUM=NCM+1 


50 

XLf< = CST=»CCS (^NGL ) 
YLN=0S7»«IK(ANGLI 
IF(YL.GT.O. ) GO TO 

^0 


yiN=-YLN 


6G 

IF(OIST.LT.C$TI gc 

TO 9C 


70 

80 


SO 


XIC( I ) = XIC( I tNUM-I) + XLN 
VIC( I tNU M levied iNUf'-D + YLN 
IF( VL.GT.C. ) GO *10 7C 

ICIR( I,NUM)=:~ ]FIX(C£G) 

GO TC 80 

ICIR(ItNUM)=IF!X(D£G) 
CIST=CIST-CST 
NLM=NIM+1 
GO TO 6C 
EX1P/5=CIST 
XLE=£XTR;i*CG$( >iNGL ) 

YLE = EXTP;i*SIN(4NGL) 

if(vl.gt.o.) gc to iOO 

YL£=-YL£ 

CCNT I^U£ 


100 _ 

lie CONTINUE 

500 FOFM^iT (36X,I2 ) 

510 F0PMAT( 12XtF8.ia2>»F8.1 ) 

R ETURN 
ENC 

SLEROUT INE LiiMCA ( I» J ,PCTV IS , C ETRAT, P K ) 

SLERCUTINE L^iMCA U CCNJKTICN WITH THE LCS RnUTINf CCMPUTES 
THE DETECTION P4TE(CETRAT) OF TARGET J BV THE OBSERVER I GIVEN 
THE PERCENT CF TARGET VISIBLE (PCTVIS) TO THf: OBSERVER. 

CChNON /CPFl/ IPFC IR (6 )f HECWC (6 ) tMV TDIR( 6),X( 6) , Y ( 6), SPO( 6) 

TCFACT=1.C 

2EF0L=C.OOCOl 

FAI=3.14159 

7 C=:( ISECWC(I) *PAI /laC.C) / 2.0 

e6E=(1.0/(2.0'*(S IN(0 )-C^CCS( C ) ) ) ) 

IF(AfeS(BES».LT-Z£RCL) eBe=O.G 
AAA=(-Bet I^COS(D) 

IF ( A6S( AAA).LI .ZERCL ) AAA=0.0 
CTANG=ATAN2((Y(J)“Y(I)),(>(J)-X(I))) 

IF(CTANG,LT.-PAI/2,AND.CT^NG.GT.-PAI ) OTANG = I^PAI + CTANG 

FD=iPRCIF (I)^PA!/i80.0 

IF((PO^aiANG).GE.C.O) GCTC 1 

1F<PD.LT ,C.O ) GOTC 9 

ANGLE=2«FAI+OTANG-FC 

GOTO 10 

9 ANGLE=2TFAI+PD-OTANG 
10 IF(ANGLE .GT. PAI) ANGLE=2»PAI-ANGLE 
GOTO 2 

1 ANGL£=AES (FO-CTANG) 

2 IF(ANGLE.GT.D) GC TO 3 
CLF=PD+0 

CLCVi = PO-C 

ANGLFT=ClAN3+( 15.C*PAI/18C. ) 

IF(ANGLFT.GT.CUP) ANGLFT=CUP 
A^GLRT=0lA^G-(l£.^FAI/180•) 

IF(ANGLRT.LT.OLCW) ANGLRT^OLOW 

FK=EBB* AES (ABS (SIN(ANGLFT J)-AeS(SlN( ANGLRT)) ) +AAA^ ( AMGLFT- 
^ANGLRT) 

IFIPK.LT.C.O) GO TC 3 
IF (PK.GT.1.0) GC TC 5 
GO TO 8 

3 PK=0.0 
DETRAT=0.C 
GO TO 6 

5 PK=1.0 

8 RANGE=SQRT (( X( J)-X( I ) )=«‘^2 H(Y ( J)-Y(I ) )^«2 ) 

RR = 0.001 GRANGE/PCTVIS 
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I 


c 


c 

c 

c 


T0ANG=ATfN2((Y(I)-V(J)),(>(I)-X(J))) 

AD = f'VTOIF( J)*P/>I/1S0.0 
HORveL=Ae'(SPD<J)»SI^(TC^^G-^CI) 

HOPVEL = HCPVEL* 16C<?. 3/2600 .0 

DENQ.M=1.4 53+TCFAC7* (0. S'?? £ + 2.128* (RR**2 »-0.503 8 *H0RVEL 1 
IFCDENCM.LE.ZEPCLt DENCM=2ERCL 
CETRAT=0.CC3*1.O 80 /OENCM 
CE1RAT=CETRAT#PK 
6 RETURN 
E^C 

SUERCUTIhE ELcV(X,Y,7HAC) 

*** SUERCUTUE EL5V CETERMINEE THE TERRAIN ELEVATION FCR A GIVEN 
SET OF X. Y CCCRCINflTES. THIS RUMCTIOM IS USED IN CONJUNCTION 
WITH THE LCS SLBPCL'TINE U CONFUTING L IN E-OF-SIGHT BETnEcN 
08SERVER AND TARGET. 

COH^l-'N /HILLS/ XC (IJC) ,YC liOO ) ,P = AK( 100) , ANGH( 100 ),SPRC( 100 ) 

COHMOM /hills/ ECC(1CC),P)X( ICC),PYY(IOO) fPXY(lOO) .base 
CQHNON /HILLS/ NHILLS 

COH.NCfl /GRID/ LST( £ .A) ,NH l( 5 , A ( , L 1ST H (15 0 ) ,KHREP(150 ) ,KTR£P 
CCHPON /GRID/ LSTC(5,A) ,NC( 5, A) ,LISTC( AGO) ,KGREP( 150 
DATA GSIZE/IOOC./ 

** FUNCTION TC CONFUTE TERRAIN ELEVATION FCR GIVEN X, Y CCCRDINATES. 
ZMAX=BASE 
IX = 1-H F!X(X/GSIZ£) 

IY=1+IFIX(Y/GSIZE) 

if(nhl(Ix,:y).ec.o) go to iso 

LS=LST(IX,IY) 

L£NC=LS-»NHL( IX,IY)-1 
CC 100 L=L£,lend 
I=L iSTH(L) 

OX=X-XC( I ) 

CY=Y-YC(I) 

cxso=ox*cx 

OYSC=OY*CY 

CXY=CX»CY 

FACTCR=PXX(I )*CXSC+PYY(I l^QYSC + PXYd )-*CXY 
IF (FACTCF.LT.-3. ) GO TO ICO 
HT = PEAK (I ) = E xp(FACTOR) 

If (FT.l£ .ZYAX) GC TG ICO 
ZHAX=riT 
ICC C0NTP4UE 
150 tyac=ZHAX 
RETURN 
ENC 

SUBROUTINE STCCh (I,RANGE,A) 

*** SLERCUTINE STCCH CSTER.NINfS THE ATTRITION CGERFICIENTS WHEN 

A USER HAS SELECTED A STOCHASTIC ATTRITION OPTION. THE CALCULATION 
IS A FUNCTION Cr THE ORIGINAL STOCHAS^ ICALLY DETERMINED ATTRITION ' 
COEFFICIENT AS WEcL AS A FUNCTION OF RANGE. 

CCMMCN /GPP 6 / ALPHA! 6 ) 

COKMCN /GfiP3/ NcU.NRU,FL(£),FC( 6 ),N0I(3),XIC(3,200),YIC(3,200), 

*ICIR(3,2CC),AV SP,ISPC 

*, IL>STAT( 6 ).II( 6 ),LCST( 6,6 ),VISFRA,VI SFR 8 ,SIZSTK, 

*SI2S*'W,NT(6) ,NF(c ) ,S?f,0!SMAX , 

*NLCSC( 6 , £ ) .VI £FR( 6,6 ) ,RMI ^TK . ? MXT K, R M NT W , RMXT’W ,OP,TCWFR,LVAFR, 
*PT7(2,3 ) .RF, FCA ( 6 , 6 ),A?0A ( 6 , 6 ),LOA( 6 , 6 ),MA(6),0FL(6) ,P0L( 6 ) 

IF (I.EG.2) GC TC 10 
A = ALPHA ( I )*((1.0-RANG£/RH>TW)**2) 

GC TC 20 

1C A = ALPHA ( I )<•( (1 .0-RANGE/P..M>TK)=»x'2) 

20 RETURN 
END 

SUERC'JTINE £TK(I .RANGE,T ) 

SUBROUTINE ETK CCMPUTSS THE EXPECTED time FOR A GIVEN FIRER TO 
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KILL A GIVEN TARGET. THE CALCULATION IS A FUNCTION OF RANGE, 
TIME OF FLIGHT FCP A RCUNC ANC HIT AND KILL PR OB ABI L IT IE S FOR 
THE FIRI^G WEAPON SYSTEM. IT IS A NUMBER THAT IS USEC IN THE 
CCMPUTATICN OF THE DETERMINISTIC ATTRITION COEFFICIENTS. 

CCMMCN /GRP2/ TA(2),T1(2) ,TH<2),TM(2 ),TF1(2),TF2(2) ,TF3(2) , 
*F(2,6),PHH(2.6),FHM(2,EI,FKH(2f6),TF(2) 

IF(I.E0.2) GO TO 5 
TF( I )=TF1( I ) 

GOTO 6 

5 IFCRANGE.GT.iOCC.C) GO TC T 

TF( n = TFl ( I ) - (TFK I )*( lOOG . C-RANG E) /I COO.O) 

GO TO 6 

7 IF(RANGE.GT.2CCO.O) GO TO 8 

TF( I J=TF2(I )-( (TF2(I)-TF1(I))*(20CC.C-RANGE)/I 000.0) 

GO TO 6 

8 TF(I) = TF3(I)-(<TF3(I)-TF 2(1) )*(3000. C-RANGE)/I 000.0) 

6 J= (PANGE-1250.0 )/fC0.C 
IF(J.CT.6) J=6 

T=TA( I)+T 1( I )-TH(I)+( ( 1H( I )^TF(I) )/PKH(I,J) ) + ((TM(I)+TF(I) )/ 
*FHM(I,J) ) + ((l.C-PHH(I,J) )/PKH( I,J ) + PhH(I,JI-P( I,Jl ) 

PEILRN 

ENC 

SUERCUTINE SGR T(I ,M) 

♦ ♦♦ SUEROUTINE SORT IS USED TO SORT TARGETS IN ASCENDING RANGE 
OPCER. THIS IS USED TO DETERMINE THE PRIORITY OF A TARGET 
FOP FIRE ALLCCATICN. 

COMMON /GFP5/ LOT (6 ,6 ) ,RQT<6,6 ) 

DO 1C J= 1 ,M 

IF(ROT( I,M) .GE.RCTdtJ )) GCTO 10 
R=POT(I,J) 

NN::LGT(I ,J) 

RQT(I,J)=POT(I,M) 

LOT(I,J)=LQT(I,M) 

ROT(I,M)=R 
LOT(I,M)=NN 
1C CONTINUE 
RETURN 
END 

SUERCUTINE KOVER(ZC•TMACT tSI2ETtZT,S.HTS,ZS,VISFRTI 
*** SUERCUTINE KCVER CETEPMINES WFAT PORTICN OF A PARTICULAL TARGET 
IS COVERED BY THE TERRAIN BSTV^EEri THE TARGET AND OBSERVER. 

THIS NUMEER IS USEC IN THE DETECTION AND ATTRITION COMPUTATION. 

IF(S.£Q.O.) GO TC 10 

IF(HTS.GE.ZS) GO TO 20 
HEXT=ZO+ (HTS-ZO) /S 
EVIST=AMAX1(HE>T,TM£CT) 

IF(EVIST.GE.ZT) GO TO 20 

IF( EVIST.LE.ZT-SIZET ) RETURN 
VI$=(ZT-EVIST )/SIZET 

IF(VIS.LT. VISFPTI VISFRT=VIS 
RETURN 

10 IF(HTS.LT.ZO) RETURN 
20 VI5FRT=G.O 
RETURN 
ENC 

SUBROUTINE LOS ( XA , YA , TMAC , T M ICA , SI Z EA ,XB , YB , TMAC B ,T MI C3 , S I ZEB, 
♦ LATCB, LBTC;3, ViSFRA.y ISFPB ) 

THIS SLBrCLTINE V^AS VRITTEN BY PRCFESSCR JAMES HARTMAN, NA\FAL 
POSTGRACLATE SCHCCL. IT COMPUTES A PERCENT OF A TARGET VISIBLE 
TO A PARTICULAR CESERVEB, GIVEN THE CCCROINATES OF BOTH 

COMMON /HILLS/ XC(IOC) ,YC (10C) , PEAK( 1CC) ,ANGH( 100 ) ,SPPO(100) 
COMMON /HILLS/ ECC (100),P>X( ICC),PYY (100),PXY(100),BASE 
COMMON /HILLS/ NHILLS 

COMMON /COVE R/ CXC(150),CYC( 15C),CPE AK(15CI,CP XX(15C),CPYY( 15 0) 
COMMON /COVER/ CPX Y (150) , NCV EIS 
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IOC 


lie 

12C 

125 

13C 


lAC 

15C 


UC 


170 


180 

19C 


20C 

C GRID 
C 


COWMCfsi /COt^TR/KH, KHW , KV, KN , KGRS » K6LL, KINT 

COl^MCN /GRID/ LST{ 5,A) tNHL{ 5, A) ,L ISTH{ i50» ,KHREP( 15C ) iKTRSP 

COMMON /GRID/ LSTC(5,4),NC(5,A),LlSTC(400)»KCPEP(i50) 

CIMENSICN IGX( ICCI .1 GY( IOC) , lELdCOI ,C£1 (100) ,CS2 (100 ) 

Otlfi GSI/E/IOOC./ 

SUBPOUTINt TG CCyPlTE FRACTION VISIBLE FOR OBSERVER TARGET PAIRS 
VISFRA = 1<, 

VISFP8 = 1 . 

>£A=XB-)i A 
Y£A=Y8-¥A 

IF ((X3A. tC.C . ) .AM . ( YBA.EO.O . ) ) RETURN 
1F(SI2EA+TMICA.LE.0.) GO TO 510 
IF(S IZES+Tvice.LE.O.) GO TO 510 
IFCrmCA. 1 T.O. ) VISFRA=1.0-*T«iCA/ SIZEA 
IF(TMTC3. LT.C. ) VI SFR3 =1.0 + TM 03/SI ZE3 
ZA^TyACA •* TMICA + JIZEA 
Ze=TFACB + IMIGB + SIZEB 
KTRE?=KTREF-t 1 
ZEA-ZB-ZA 
XeA£C--Xfi A»XE A 
YBASC^YBA^YEA 
XY EA-XEA^YEa 
TW 0SBA=2.*>EA 
TLGt*iA = 2 .*YEA 

CONFUTE GRIC SCUAPES CROSSED BY A TO E LINE 
NGPSO=0 

IFIXBAi 310,95,100 
X3A=0.l 
ISG>=-1 

XINC=GSIZE/XRA 
GC TO 12C 
ISGX=1 

XINC=-GSIZE/XEA 
IF(YBA) 1^3,125,330 
YBA-C.l 
I£GY=-1 

YINC=GSJZE/Yf A 
GO ~0 15C 
ISG¥=1 

Y1NC = -G SIZE/YBA 
IX = 1+IF IA (X5/GSIZE1 
n = l-^iFI>fYc5/GS!ZE) 

XN6XT=C-S IZS» <FLCAY( IX) -t-O. E*( ISGX-1. ) ) 

YNEX’=GSIZE*{FLCA~(IY)+0.5*( ISGY-1.) ) 

XST = r-=( XE-VI = XT) /X9A 
YSTgP=(Ye-YM£xT)/veA 
hGF£a=:JGFSC*-i 
IGXINGR EC ) •'X 
IC-Y(NGkSC )=IY 

IF((XST£C.GT.1.KAND.(YSTEP.GT.1. )) GC TO 200 
IF(XSTEP-fSTEP ) 17C, iSC,1EO 

IX=JX ^-ISGX 
XSIE?=XSTEF+XINC 
GO TO ICC 
IX= IXf!SGX 
XSirP=.<SlEP + .«INC 
IY= IY-*I SGY 
Y31FF=YITEP+YINC 
GO TC ICC 
KGRS=KGRS+NGPSC 

SCUARE list ACV« COMPLETE IN IGX, IGY WITH NGRSO ENTRIES 


FINE WHICF COVER ELLIPSES TOUCH THE A TO S LINE, 
CHECK ELEVATIONS AT SI AND S2 FOR EACH SUCH ELLIPSE 
NELS=0 
CHTNAX=C. 

IF(NCV£LS.eO.C) GOTO 270 
DO 2£C X=1,AGFSC 
IX = IGX(K ) 

I Y=IGY(K ) 

N = NC(IX, lY ) 

IF(N.EC.O) GC TO 26C 
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L£=LSTCiIXTlYI 
LEN0=LS4fx-l 
CQ 25C L=LS,LENC 
.‘<ELL=KELL+ I 
IC=LI£TC(L) 

IF (KCPEP (Id . EC.KTPcP ) 3C TO 250 
KCREP(IC)=KTREC 
PX=X;^-CXC ( !C) 

RY = >4-CYC( IC) 

PPX)i = CPXX( IC) 

PFYY=CPYY ( IC ) 

PFXY = CPXY( IC) 

A/1 = PPXX^X8 ^SQ + PPYY^ Y9 ^SC+PPXY*XY8A 
BE=FFXX'-^TWOX EA*RX + PPYY^ TWCY BA^RY + po XY* ( RX’^YSA+R Y’^'XBA ) 
CC = FP>X + R X*RX + PFYY ‘KY*RY-^ FPXY^RX ^ RY-1.0 
APG=eE’*e0--.O^^AA'*CC 
IF (ARC.L5.C.) GC TO 2 50 
SC = SCPT (ARC » 

Sl=-(Be^SO)/(2cC*AA» 

£2=(SQ-F3)/(2.0^AAI 
IF( Sl.GE. 1. ) GO TO 25C 
IF(S2.LE.O^ ) GO TO 25C 
IF(S1. LE.O* ) CO TC 510 
IF(S2,GE.l. ) GO TO 51C 
C*** CHECK LOS AT SI ANC S2 
K INT = KINT + 1 
CFK=CFEAK(IC) 

XS=>A+S2^>eA 
YS=>A-»>S2’»'r0A 
CALL ELEV (XSfYStHTS ) 

H7 S=HT$-*-C FK 
ZS = 2A4S2-* Z0A 

IF( lATOB.EC.C) GC TC 210 

CALL KCVER(ZA,TMAC3,S I2c6 tTS t$2•HT$,ZS,VISFR6) 
IF(VISFRe.L£.Od CO TC 510 
21C IF( IBTOA,EC.C) GO TC 220 

S=l.-S2 


22C 


23C 


2AC 


25C 

ItC 

270 




CALL KOVER(2 5.TMA(A,3IZEAtZA,SfHTStZS» VISFRA) 
IF(VISFRA.L£.0.) CQ TC 510 
XS = >A+Sl*'^xeA 
YS=YA+S1^Y6A 
CALL ELE\(>S,YS,HTSi 
HTS=HTS^CFK 
ZS=2A+S1^2eA 

IF(LATO0.EO.C) GO TQ 230 

CALL KCVE?(2A,TMACe,S I2E3.20 »S 1 oHTSiZSiVISFR 0 ) 
IF(VISFR0.LE.O.) GC tq 510 
IF(LETOA,EG.C) GO TQ 240 
5 = 1.0-S1 

CALL KOVER (Z3 , TMACA, STZEA ,ZA ,S ,HTS, ZS, VI SFRA)- 
IF(V ISFRA.LE.O-) GO TO 510 
NELS = NELS^1 
lEL(NELS)=IC 
CSl (NELE .)^Si 
CS2(NELS)=S2 

IF (CPK,GT^CHTMAX) CHTMaX = CPK 
CCNfTINUE 
CONTINUE 

ALL ELL IPSES CHECKED 
START CN THE h^ILLS 
CO ECO K=1 tNGRSQ 
IX=IG> (K ) 

IY= IG>r ( K) 

IF(NHL(IX,IY)*EC.C) GC TO ECO 
LS=LST (IX, lY) 

LENO=lS-^NHL(IX,IY) -1 
DO 5CC L=LS,LENC 
I = L I ST H ( L ) 

I F (KHREF (I ) •EC*KTRE F ) GO TO 50 0 
KHF5P( I )=KTREP 

FFOCESSING FOR FILL I STARTS HEFE 


KH=KH+1 
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c 


28 C 

3CC 

C*** 


35C 


38C 

3SC 


C* + * 

ACC 


COMPUTE W =TOP OF THIS HI LI ALONG 0-T LINE 

TR>=XA-:XC( I ) 

TP>=YA-YC(I) 

TF>;( = PXX (I ) 

TF^Y=PYY(I ) 

T F>Y = PXY (I ) 

FO=TWOXBAvTFXX«TRX + 'rwCYe 3^TPYY ’»TRY>TPXY=^(TRX^Y2/ ^ 
GQ = TPXX+ >8ASQ+TPYY* >BASC>TPXY* XY6A 
IF (GC .EO.O.) GO TO 500 
W = -FC/(2.'«'GC) 

IF(A0S(^^}.GT.5.) GO TO 500 
FSC=FQ«FC 

EQ = TPXX^TPX* TR X+TPY V'^TP Y*TRY+T PXY^c TRX^TRY 

FGWEP=EG-FS0/(A.*GC) 

IF(PCWER .LT. -3*) GO TO 500 
HHW = FEAK (I )*eXF(POW iR) 

W=Km W-f 1 

IF(HHW.LE.BASE) go to 5CC 
ZW = 2A+U*2eA 

IF((W.LT.0.).aR.(V^.GT.l.)) GO TC300 

IFIHHW •6E. ZU) GO' 10 5 1C 

CVH1V=0. 

IF(N£LS.6C.0J GO TC 30C 
DC 280 M=1,NEL$ 

IF( (CSl(M).GE.W ).CR.(CS2(M) .LE.W)) GO TO 280 
IC=IEL(M) 

IF(CVHTW*LT.CPEAK(IC) ) CVHTW=CPEAK{IC) 
CCMINLE 

IFKHHW^CVHT;.) .GE.ZM GC TC 510 

I FIHFW + CHTMAX.LT.AM ink ZA-SI2E a, 23-SIZEB) ) GO TO 
IF WE GET TC HERE THEN NEEC TC FIND LOWEST SIGHT LINE C\ 
NEWTCN ITERATION A TC 6 GIVING VI5fR8 
IF(LATOB.£C.O) GC T0 400 
KV=KV+1 
V = W 

HHV=HH W 
NCT = 0 
FV = FC*V 
TW0GV=2.*GC^V 

FCNV = Z A+HHV’*(TWQGV+ \ + FV- 1. ) 

KN=KN+1 

FACTCR=(TKOGV^ThCGV ^2.*(GC+TWOGV*FG ) +FSQ) 
OFCNV=hHV*V=*FACTCR 

IF <AES ( OFCNV ) .LT.l.E-lC) GC TC 330 

V=V-FCNV/DFCNV 

FV=fC*V 

TW0GV=2.^GC*V 

POVER = EQ^FV+GQ^V^V 

IF(PCWER .LT. -3.) GO TO 400 

HHV=PEAK{I)=*EXF(POW iR) 

DHHV=HHV+(FO+TVaGVJ 

ELV=2A>CHHV^V 

IF (ABS(HHV-ELV) .Ll.l.) GO TO 350 
N01 = NCT> 1 

I F (NCT.LT.IO ) GO TO 33C 
IF((V.LT.C.).CR.(V.GT.l.)) GO TG400 
CVHTV=0. 

IF (NELS. EC.O ) GO TO 390 
DC 360 M=1,NELS 

IF( (CSKM l.GE.V ) .OR. (CS2(H) .LE.V ) )G0 TO 360 
IC=IEL(N) 

IF(CVHTV.LT.CPEAK(IC)) CVHTV=CPEAK(IC) 

CONTINUE 
HTV=HHV^-CVHT V 
ZV=ZA+V^ZBA 

CALL KCVER(Z/tTMACe »SIZEe,ZB, V tHTVtZVtVISFRB) 

IF (VISFPB.LE.C.) GC TC 510 
NEWTON ITERATION B TC A GIVING VISFRA 
IF ( AES(V ).GT.5.)G0 TC 4CC 
IF(L3TOA.EC.O) GC TC 5CC 


TR f^XBA) 


5 00 

3R HILL 
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KV = KV+ 1 

VM1=V-K 
HHV=hHW 
NC1=0 
fV=FQ*V 
TW0GV=2.^GC^V 

^3C FCNV=26 + I-HV» < (FO+TWOGV . ) 

KN = KN+ 1 

FACTCR=(TWGGV^TWCGV^2. -*< GQ+TWOGV>^FQ )+FSC) 

OFCNV=hHV-?^VMl^FACTC P 

IF(ABS<CFCNV) •LT,1.E-10) GG TC 450 

V=V-FCKV/CFCNV 

IF(AeS(V).GT.5.)G0 TO 5CC 

VM 3 = V-!• 

FV-FC*V 

TWCGV=2. ♦GC’^V 

PCWER = EQ^FV + GQ*V*\/ 

IF{FCWEP .LT. -3,) CO TC 50C 
hhv=feak(I)^e;<p(fgvn ir) 

OHHV=HHV’*( FG+TXGGV ) 

El V=23+C^hV=>►V^^l 

IF (ABS(HHV-ELV) .LT.l.) GO TQ 450 
NC1=NCT41 

IF (NCT.LT.IC) GC TC 43C 

45C IF((V. LT,0, ).CP. (V.GT.l. )) GO TC500 

CVFTV=0. 

IF (NEL S.eC.C) GC TC 490 
DC 4S0 f-=ltNELS 

IF( <CS1 ).GE.V ) .GR. (CS2<M) .LE.V ))G0 TO 480 
IC=IEL<K) 

1F(CVHTV.lt.CPEAK(I C) ) CVHTV = CPEAK(IC) 

<ec CC^TI^^UE 

4SC HT\=HHV+CVHTV 

ZV = 2A^-V»2aA 
S = -VM 

CALL K0VER(2BtTMACA ,Sl2EA»ZAt$ ,HTV,2V,VI5FRA) 

IF (VISFPA.ie.O.) GC TC 510 
'CC CONTINUE 

60C CONTINUE 
RETURN 

51C VI5FRA=0. 

V I$FR0 = C . 

FETUPN 

END 
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APPENDIX C 


COMPLETE INPUT DATA SET 

for the 

SMALL-UNIT AMPHIBIOUS OPERATION COMBAT MODEL 

The small-unit amphibious operation combat model consists of two 
phases of combat, ship-to-shore and land combat, and requires data 
input for each of these phases. The data set that follows is divided 
into two parts: the first part consists of all data used as input for 
the ship-to-shore phase of combat, and the second part consists of all 
data used as input for the land combat phase of combat. The input data 
set was designed to be self-documenting in that the input variable 
names or descriptive phrases are listed alongside the data being used 
as input to the model. The purpose of this documentation was to assist 
the user in associating the input data with their respective input 
variables. A complete input data set follows. 
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IPRINT 

= 1 

C FOR EACH 

TIME STEP S 

1 FOR END 

OF BATTLE 

. 

LVA‘5 

5P0MAX = 

40.00 5 

POMIN =10. 

5 HTMAX = 1. 7 

HIM IN = 

0.6 V.IOTH 

= 3.533 


TT5 = 

10. 

LENGTH OF 

EACH TIME STEP IM SECONCS. 


TANK 

MAX RANGE = 

150C. 

ATGM MAX RANGE = 2000. 

ATGM MIN 

RANGE = 

200. 

ISIGV 

TARIM = 

15. 

SARTM = 30 

. TVEL = 350 

. SV5L = 

35C. 


- c 

- ^ • 


500. 

1000 . 

20 00 . 

5000. 

10000. 

0. 




5. 

20. 

25. 

25. 



TS IGH 

’ 2 5. 


500. 

1000 . 

20 00. 

5000. 

ICOOO. 

0. 


c . 


c 

^ . 

20 . 

25. 

25. 



TMENH 

^ 25. 


5CC. 

ICOC. 

20C0. 

5000. 

lOCOO. 

0. 


I. 


5. 

10 . 

15. 

15. 



5SIGV 

^ 2f. 


25C. 

50 0. 

1000. 

2500. 

5000. 

10000. 


c . 


5. 

7.5 

14. 

15. 5 

17. 

20. 

SSIGH 

- 2 5. 


250. 

500 . 

1000. 

2500. 

5000. 

10000. 


C. 


^ . 

7.5 

14. 

15. 5 

17. 

20. 

CEF. 

WEIGHTS 

ASSIGNEE TO 

WAVE CNE 

= 2. AND 

WAVE TWO 

= 1. 



NO. OF 

LVA 

PER ^.A VE = 25 . 

20. 

15. 

10. 

5. 

SHOP 

E OEF. 

TANK 

ASSETS 

= 

10. ATGM 

ASSETS = 

10. 



^ITkIT CC‘EF* fc? 
iSllKli Cwhh« FC^ 
AI^EC FIFE 
GAIN 
GAMMA 
NO. CF HI 
BA 

HILL CAT/.tXC = 


AGGPESSOF FORCE FIRE: ALPFA1= 0.00006 ALPHA2 = O.CCCOS 

JtFENSlvE FCRCt FIRE: BtlAi == 0.0007 &ETA2 = 0.0009 

ATTRIT CCEF.: weETA( 1 ) = C.0005 Vi8ETA(2) = 0.0007 


LST(5,4) X 


11 = 0.32 

= 50.0 DELTA 

LLS IN TERRA IN MO CEL 
.<E = G.C 

2C00. YC= 1100. PEAK 

= IOC. 0 
= 46 

= 170. ANGH= 

0. 1 

SPRO = 

999.9 

3 SCO. 

2200. 

150. 

30. 


350. 

:cco. 

1900. 

150. 

130. 


3C0. 

2AC0. 

1400. 

150. 

0. 1 


30C. 

2450. 

1700. 

130 . 

80. 


SCO. 

27CO. 

leco . 

138. 

90. 


500. 

3 200 . 

1650. 

140. 

150. 


6CC. 

K300 . 

1300. 

130. 

160. 


400. 

3 750. 

17 50. 

150. 

0. 1 


660. 

^ 150. 

16C0. 

150. 

160. 


550. 

3 2CO. 

2150. 

130 . 

25. 


500. 

^ ECO. 

1700. 

170. 

45. 


5C0. 

^800 . 

1500. 

170. 

0.1 


3CC. 

d dCO • 

2600. 

170. 

90. 


350. 

2 AOO . 

2850. 

150. 

120. 


3C0. 

3100. 

2700. 

150. 

150. 


350. 

2 500. 

24CC. 

150. 

0. 1 


250. 

2 e 50 . 

2850. 

150. 

160.0 


4C0. 

27CO. 

26C0. 

150 . 

130. 


370. 

3eco. 

2 2 C 0. 

150. 

0. 1 


230. 

^500. 

2600. 

150 . 

90. 


280. 

3eco. 

26C0. 

150. 

145. 


500. 

2 700. 

3300. 

190. 

25. 


350. 

3000 . 

33CO. 

170 . 

15. 


400. 

3150. 

37 5C. 

130. 

0.1 


350. 

3 750. 

32C0. 

150. 

10. 


850. 

3eco. 

28CO. 

150. 

0.1 


650. 

3 ECO. 

36C0. 

350. 

160. 


320. 

^150 . 

3950 . 

170. 

30. 


220. 

1650. 

21C0. 

150. 

30. 


300. 

2250 . 

2100. 

180. 

150. 


220. 

22C0. 

1700. 

180 . 

80. 


150. 

26CO. 

1400. 

150. 

165. 


500. 

3 600. 

2650. 

150. 

20. 


400. 

^150. 

3CC0. 

130. 

160. 


300. 

^^50. 

3150. 

130. 

140. 


350. 

^850 . 

3600. 

150. 

45. 


400. 

^500. 

40C0. 

150. 

60. 


450. 

A850. 

3450. 

150. 

150. 


440. 

A900. 

3000. 

130. 

45. 


250. 

A7C0. 

24C0. 

150. 

50. 


400. 

4900. 

20C0. 

130. 

0. 1 


25C. 

25CO. 

1100. 

170. 

0. 1 


950. 

5350. 

43C0. 

150. 

45. 


450. 

4C00. 

2200. 

150. 

45 • 


280. 

39C0. 

22C0. 

150. 

0. 1 


300. 

= 0 

C 0 

COO 

1 

7 

18 


ECC = 


27 


8.0 

2. C 

2 . 

2. 5 
2.2 
2.2 

3. 


1 . 

2 . 

2 . 

1 . 

1. 

2 . 

1 . 

3. 

1. 

1. 

1. 

2 . 

2 . 

2 . 

2 . 


2 . 
2 . 
2 . 
1. 
1 . 
2 . 
2 . 
1. 
2 . 
1 . 
2 . 
1. 
1. 
2 . 
1 . 
8 . 

2.5 

2 . 

3.5 
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VJiuiinoo CDUi“>iui-^vnfvjOf\jO ov-ninouiUJUicDOO ooa>v.iiuiv.n ovn 














0 

33 

29 

NHL(5,A) ^ C 

C 

0 

0 


lA 

NO. OF HILLS TOTAL 
LISIH(I) = i 

= 


2 

3 

6 

c2 

33 

8 

G 

10 

8 

A2 


16 

17 

18 

2 

31 

11 

A6 

2 C 

^ 1 

A2 

A5 

46 

26 

lA 

25 

35 

A A 

26 

AO 




NCVHLS = 


62 

0 

7A 

77 

82 

93 

0 

0 

6 

11 

9 

6 

12 

C 

3 

6 

10 

9 

101 

A 

A3 

1 

3 

A 

5 

11 

31 

A3 

1 

6 

7 

33 

A3 

10 

12 

12 

9 

30 

23 

15 

3 

14 

15 

20 

*3 

6 

23 

11 

7 

20 

22 

34 

35 

AA 

45 

12 

3A 

35 

36 

AO 

41 

22 

15 

23 

24 

25 

15 

27 

28 

29 

2A 

22 

22 

29 

35 

36 

37 

38 

39 
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AT7RIT 

VAR = 

0 CSEED 


IA3257.0 


AGGPESSOP 

FORCES 1 

BETA CIET* 

INPUT 

PARAMETERS: 

PP 


7. c 

00 ~ 

21.0 

ATTP IT 

VAR s 

0 CSEED 


123457.0 

(F0« 

DEFENSIVE 

FORCES) 

BETA CIST* 

INPUT 

PARAMETEPS: 

PD 


21.0 

00 = 

7.0 


NC. CEf UNITS = 03 NO. iTK UMTS = 03 

PMINTK = OCCC.O PMXTK = 250C.0 PMNTW 
type CF route = 1 VEHICLE IPEEC = 4 


= 0500,0 


RMXTH = 4000.0 


XIC(1,1> = 2C00.0 

YIC(1,1) = 

19C0.0 







ISCO.C 


2400.0 







15C0.C 


2 ICO.0 







NO. CF NCCE5 FOP 

ROUTE 1 = 

0 1 







XlOC(ltll= fCCO.C 

YLOCd *i) ^ 

25C0. 0 

FCR 

NODE 

1 

OF 

FCUTE 

1 

NO. CF NODES FOR 

ROUTE 2 = 

01 







XLCC(2,1)= ASCO.O 

YLGC (2,1) = 

2150.0 

FOR 

NODE 

1 

OF 

ROUTE 

2 

NO. CF NODES FOR 

PCUTE 2 = 

02 







XLCC(3fl)= S2C0.a 

YLQC(3»1 )= 

ITCO.0 

FOR 

NHOE 


CF 

FCUTE 

3 

>L0C(3f2)= ^eco.c 

YLOC (3 ,2 ) = 

1750.0 

FOR 

MODE 

2 

OF 

FOUTE 

3 


DEF UNIT X = 

t 

LOCATICN 

ALT. POS. 


5EC0.C Y = 2'rCC. 0 FCPCEIO.O OIP. IFC WIDTH 120 

2800.0 22C0.0 LEVEL 5.U C'F 19C OF 120 

EtCO.O ITCO.O 10.0 FIRE 180 SkCH 120 

VAR. 0 BREAK PI CECO.C NC. TIME STEPS FCR MOVE 4 

FCP ALT. POS. 1 
= = ^2 


XA(1) 

= 4500 

.0 

YA(1) = 

3600.0 

XA(2) 

= 45CO 

.C 

YA(2) - 

27 CO.0 

X 

= 46C0 

.0 

Y = 

leco .0 

Q.6C 

0. 

70 

C.65 

0.85 

0.65 

C. 

90 

C. 65 

C.90 

0.60 

0. 

85 

0.85 

C.80 

0.75 

C. 

8C 

C.75 

C.70 

0.60 

c. 

7C 

C.65 

C.65 

0.40 

0. 

45 

0.40 

C.50 

C.F5 

PHH C. 

65 

FHM C.75 

PKH C.70 

0.60 

C. 

80 

0. 75 

C.70 

C.15 

0. 

75 

0. 70 

C.60 

0.60 

C. 

65 

C. 60 

C.55 

C.45 

0. 

50 

0.50 

C.35 

C.2C 

C. 

2C 

C.20 

C.20 
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APPENDIX D 


BLANK INPUT DATA SET 

for the 

SMALL-UNIT AMPHIBIOUS OPERATION COMBAT MODEL 

The blank data set provided with the small-unit amphibious opera¬ 
tion combat model was designed to assist the more familiar user of the 
model in the development of a new input data set to be analyzed by the 
model. It is patterned after the complete input data set listed in 
Appendix C providing input variable names or descriptive phrases to 
identify the locations of required input parameters. Underlining of 
the spaces following these descriptors is intended to serve as a guide 
for inputing values for the input variables in order that they will be 
compatible with the formatted read statements of the program. The 
blank input data set follows. 
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IPRIM = 

CVA'S 5FCMA)( = 

TTS = _ 
TANK MAX RANGE = ' 

TARTM = _ 


0 FOR E/CCH TIME STEP S 1 FOR END CF BATTLE, 
SPOMIN HTMAX = .. HTMIN = . WIDTH = 

LENGTH CF 1/iCH TIME STcP IN ScCCNCS. 

ATGM MA> range = _• ATGM MIN FANGE = 


SAPTM 


TV EL =“ 


SVEL =_, 


TSIGH = _; 


„ ^ • 

—• 




TMENh = JZl _I 


, # 

—• 


—! 


SSIGV = _I 


_ ^ . 

-T- . 

— • 


— 1 


SSIGH = JZl „Z 



— • 


—* 

- 

DEF. WEIGHTS ASSIGN? 

C TO W#VE 

It 

fUJ 

1 M 

1 


AND 

WA7?* 

Two 

NC. OF LVA PER 

WVE = _ 

_ • 





SHORE CEF. TANK ASS 

ETS = 

_ 

ATGM ASSETS 

ATTRIT CCEP. POP AGGRtSSCP FG 

RC E FIRE: ALP KA1= 



ATTRIT CCEF. FGR OEF 

ENSIVE PC 

RCE FIRE: BttAi = 



AIMEC F :rE ATTRIT COEF. 
GAIN! = 


WBETAC1) 


• . __ 


GAMMA = 


DELTA 

* -- - - • ^ 




NO. CF HILLS JK 

TERRAIN 

MO CEL 





BASE = C 

.0 






HILL DATA:XC= 2000. 

YC= 1100. 

PEAK = 

170. ANGH= 

0. 1 

SPRI 

18C0. 

2200. 


150. 


30. 


2CC0. 

19C0. 


150. 


130. 


24CO. 

14C0. 


150. 


0.1 


2450. 

1700. 


130. 


80. 


27C0. 

18C0 

. 

138. 


90. 


3 200. 

16fO. 


140. 


150. 


4300. 

13C0. 


130. 


160. 


3750. 

17fO. 


150. 


0. 1 


4150. 

1600. 


150. 


160. 


32CO. 

2150. 


130. 


25. 


46C0. 

17C0. 


170. 


45. 


4800. 

1500. 


170. 


0.1 


22C0. 

26CO. 


170. 


90. 


2400. 

2850. 


150. 


120. 


3100. 

27CC. 


150. 


150. 


2500. 

24C0. 


150. 


0. 1 


2650. 

2850. 


150. 


160.0 


27CO. 

26CC. 


150. 


130. 


3aco. 

2 200. 


150. 


0. 1 


4500. 

2600. 


150. 


90. 


36CO. 

28C0. 


150. 


145. 


2700. 

33C0. 


190. 


25. 



= . 

2ETA2 = I.; 
WBETA{2) = 


LST(5,4) s 


3CCO, 
3150, 
3750, 
36CO, 
3 600, 
4150, 
1650, 
2250, 
^ CO , 
2600, 

3 800, 
M50, 
4450, 
4850 , 
4500, 
4850, 
4<:C0, 
47C0, 
4900, 
25CO, 
3350, 

4 000, 
39CO, 


3300. 
3750. 
3200. 
3aco. 
2600. 
3 «^ 50 . 
21CC. 
2100 . 
1700. 
14C0. 
2650. 
3CC0. 
2150. 
3600. 
40CO. 
2450. 
30C0. 
24C0. 
2000 . 
1100 . 
4300. 
2200 . 
22C0. 
0 ( 


170. 
130. 
150. 
150. 
150. 
170. 
150. 
130. 
180. 
150. 
150 • 
130. 
130. 
150. 
150. 
150 . 
130. 
150. 
130. 
170. 
150. 
150. 
150. 
0 


15. 

0.1 

10 . 

0.1 

160. 

30. 

30. 

150. 

80. 

165. 

20 . 

160. 

140. 

45. 

60. 

150. 

45. 

50. 

0.1 

0.1 

45. 

45. 

0.1 

1 


<=99.9 

350. 

300. 

300. 

5C0. 

500. 

600. 

400. 

660. 

550. 

500. 

200 . 

3CC. 

350. 

300. 

350. 

250. 

400. 

370. 

230. 

280. 

500. 

350. 

400. 

350. 

850. 

650. 

320. 

220 . 

200 . 

2 20 . 
150. 
500. 
4C0. 
300. 
350. 
400. 
450. 
4^0. 
250. 
400. 
250. 
950. 
450. 
280. 
300. 
18 


ECC = 


27 
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0 

33 

39 

NHl(5»4) = C 

C 

0 

C 

6 

14 

NO. CF HILLS TCT/M 

= 


LIMHin = 1 

2 


6 

32 

33 

3 

9 

10 

a 

42 

2 

16 

17 

18 

2 

21 

11 

A6 

2C 

21 

42 

45 

46 

26 

14 

25 

35 

44 

26 

40 




NCVELS = 


62 

0 

74 

77 

83 

93 

0 

0 

6 

11 

9 

6 

12 

0 

3 

6 

10 

9 

101 

4 

43 

1 

3 

4 

5 

11 

31 

43 

1 

6 

7 

33 

43 

10 

12 

13 

9 

30 

23 

15 

3 

14 

15 

20 

■3 

6 

23 

11 

7 

20 

22 

34 

35 

44 

45 

12 

34 

35 

36 

40 

41 

23 

15 

23 

24 

25 

15 

27 

28 

29 

24 

22 

23 

29 

35 

36 

37 

38 

39 


0 


53 

0 

S 

30 

7 

11 

lA 

IS 

16 

22 

1 ^ 

26 

27 
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JlTTPIT VAR = ^ nSEEO _ 

BETA CIST, INPUT p■APA^£TEPS: FP = 

ATTFIT VAP r CSEED = 

BETA CIST. INPUT PARAVETEPS: PO = 

NC. CEF UNITS ^ _ NC. ATK JMTS = ^ 

PMINTK =_.. RN;<T!< =_. FRIRTW 

TYPE CF ROUTE = _ VEHICLE SPEED 
Xicdil) =_.. YlC(lfl) ^_ 


<FCR AGGRESSCR FORCES) 
_•_ QQ ” 

IfZR DEFENSIVE fQPZzSl 
00 = 

=_.. RMXTW =_ 


NC. CF NODE? F03 POUTS , = 
XLCC(ltl)= VlCCTltl 

NC. CF NODE? 73R ROUTE = 
XL0C(2tl)= VLCCl2tl 

NO. CF NODE? r3R ROL'^E 


XLCC(2 tl) = 
>LCC(2 ,2) = 
CEF UNIT X = 

I 

LCCATICN 

ALT. POS. 
XA (1) = 
XA(2) = 

X = 


— _ 

varT*I 


AZZ15.Z 


PhP 


>LCCT3f1 
yLCC(3,2 
Y 


BREAK 
YA(1) 
YA( 2) 
Y 


F r R 


I = 

) = 

) = 
) = 

PT 


FCR NODE 1 CF ROUTE 1 

FCR NODE 1 OF ROUTE 2 

FCR NODE 1 OF ROUTE 3 
FCR NODE 2 CF ROUTE 3 

FCRCE_.. DIR._WIDTH 

LEVEL„._ _OF 

FIRl SRCH 

NC. tTme steps f315 move 

FCR ALT. POS. 1 

= = = 2 

= = = 3 


PKH I 
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APPENDIX £ 


EXECUTIVE PROGRAM 
for the 

SMALL-UNIT AjMPHIBIOUS OPERATION COMBAT MODEL 


The combat model presented in this thesis has been provided with an 
EXEC program which is designed to set uo and execute all of the necessary 
CMS commands for the running of the model. The EXEC program will auto¬ 
matically BROWSE the output listing of the model [AMPHIBl LISTING) 
allowing the user to review immediately the results of the battle. 

A listing of the EXEC follows. 

GLOBAL TXTLIB F0RTM0D2 M0D2EEH IMSLSP NONIMSL CMSLIB 

FILEDEF 05 DISK SEA DATA 

FILEDEF 09 DISK LAND DATA 

FILEDEF 06 DISK AMPHIBl LISTING 

LOAD AMPHIB (START) 

BROWSE AMPHIB 1 LISTING 
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APPENDIX F 


COMPUTER OUTPUT 
for the 

SMALL-UNIT AMPHIBIOUS OPERATION COMBAT MODEL 

The computer output for the small-unit amphibious operation combat 
model was designed to be clear, concise, and identifiable to the user 
of the program. The combat model conducts two phases of combat: ship- 
to-shore and land combat. Therefore, the computer output was designee 
to report on each phase of combat. The computer output for each phase 
of combat begins with an initial information page which lists the 
input data provided by the user of the model. The initial inforir.atior 
page serves as a record of the battle scenario analyzed by tne ;T.odel, 
as well as a check for the user to insure that the input data provided 
were read correctly by the model. In addition, a battle summary 
report is provided reporting on the status of both the aggressor and 
defender forces throughout both phases of combat. The computer output 
based upon the input data listed in Appendix C is as follovvs. 
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*♦ INITIAL SHIF-T0-5H0FE PHASE INFOFMATION 

INITIAL FCFCE STRENGTH 
VAVE 1 2 3 A 

LVA 25*0 20.0 15. 0 10.0 


CtF. TANK ASSETS 


10.0 


5 

5.0 

OEF . ATOM ASSETS = 10.0 


LVA ENGR SPECS 

SPDMAX SPOMIN HTHAX HTMIN WIC 

AC.00 10.5C l.TC C.6C 1.52 

CEFENSIVE TACTICAL PARAHETERS 

RANGE AIH-RELCAO PROJECTILE 
MAX MIN TIME VELOCITY 

TANK 1 500.0 15.LO 350.UO 

ATOM 2CCC.0 2CC.0 2C.CC 25C.C0 

DEFENSIVE TACTICAL ALLCCATICN WEIGHTS: 

KAVE 1 = 2.0C WAVE 2 = l.OC 

CEFENSIVE FORCE ATTRITION COEFFICIENTS 
ALPHA'^A BETA^A 

OT 0.0000^ C.00070 

DS O.OOOCE 0.C0090 

AIMED FIFE ATTFITICN RATE COEFFICIENTS FOR 
CEFENSIVE TANK AND ATOM AESETS 

WEETAd ) = O.OOC50 WBETA ( 2 ) = C. COO 70 

BREAKPOINT assumption: 0.3'^(TCTAL OEF FORCE) 

CEFENDEF ATTRITICN LEVEL ALLOWING FOR LAND COMBAT 
0.22^(TCTAL DEFENDER FORCE) 

ARTM SUP FACTOR= 50.0 ERROR SUP FACT0R = 10C.0 

DISPERSION DATA 


RANGE 

TSIGV 

RANGE 

TSIGH 

R4NGE 

tmeanm 

25.0 

0. C 

25.C 

C.O 

25. 0 

0. 0 

500.C 

2.0 

5CC .C 

2.0 

500.a 

1.0 

ICOC.C 

5. C 

lOCC.C 

5.0 

1COO.C 

5.0 

2CC0.0 

20. C 

20C0.C 

20.0 

2 COO.0 

10. 0 

5C00.0 

25. C 

5000.0 

25.0 

5 COO.C 

15.0 

lOCCC.C 

25. C 

ICOCC.C 

25.0 

lOCOO.C 

15.0 

RANGE 

SSIGV 

RANGE 

SS IGH 



25.C 

0. 0 

25.C 

0.0 



250.C 

5. 0 

25C .0 

5.0 



500.C 

7.5 

500.C 

7.5 



ICCC.C 

14. C 

lOCC.C 

14.0 



2500.C 

15.5 

2500 .C 

15.5 



5C0C.C 

17. C 

5000 .0 

17.0 



lOCOC.C 

20. C 

lOOCC.C 

2C.0 



CURRENT 

STATUS 

OF WAVE 

I VARIABLE CEFINITIOf 


C - NOT ENGAGING 

1 - LANCED 

2 - UNDER FIRE BY ATGM 

3 - UNCER FIRE EY TANK 

A - UNCER FIRE EY BOTH ATGM C TANK 
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<♦♦♦♦ the ship-to-shcre phase begins 

BREAKPOINT REACHED AT TIME = 502.5 SECONDS 


Tiy€= 5C2.5 SECONDS 


W/VE 

FCFCE LEVEL 

STATUS 

LOST-PCT TOTAL SURVIVING 

1 

16.0000 

1 

0.360 

2 

12.OC00 

1 

0.35C 

*3 

11.0000 

1 

0.267 

A 

10.0000 

1 

0 .0 

c 

TANK 

5.OC00 

0.0 

C 

0.0 55.00 

1 .000 

ATGM 

C.O 


1 .000 0.0 

final 

LVA SURVIVORS 

ASHORE 

= 55.CCC 

LAND 

COMBAT STARTS 

WHILE SHCBE CCMBAT IS GOING ON 

LANC 

COMBAT ATTACK 

TIME = 

A95.C SECCNDS 
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*< INITIAL LAND CCH6AT INFORMATION 


LOCATION 


UN IT 

X 

Y 

FORCE LEV! 

i 

2COO.O 

190C.0 

18.0 

2 

1^00.0 

2400 .0 

1 8.0 

a 

l^CG.O 

21C0.0 

1 6.0 

4 

3600.0 

2700.0 

1 0.0 

c 

36CO.C 

2300.0 

5.0 

6 

3600.0 

i7CC.0 

10.0 


ATTRITICN IS STOCHASTIC 

ROUTES DETERMINED EY USER 

ATTACK \iEHICLE SPEED IS 18.0 M.P.H. 

EREAKPOINT DISTANCE IS 500.0 METERS 

DEFENCE? WILL MOVE TO ALTERNATE POSITIONS 
ALTERNATE POSITIC.NS ARE: 

UNIT X Y 

A AfOO.C 3800.0 

i A5C0.C 2700.0 

6 AECO.O 180C.C 

ATK ;:ILL I'P 08 A9 IL IT I E S 


PAfsGE 

P 

PHH 

PHM 

PKH 

50.0 

C.60 

0.70 

0.65 

0.85 

1 ecG 

0.65 

C .50 

0 .85 

C.90 

i5CC 

0.80 

0.65 

C.85 

0. 8C 

2C00 

C.75 

C .80 

0.75 

C.70 

2 500 

C.6C 

0.7C 

C.65 

0.65 

3 COO 

0.40 

0.^5 

C.40 

C. 5C 

DEF 

. KILL 

FFCBABILITI 

ES 

RANGE 

P 

PHH 


PKH 

5G0 

C.85 

0.85 

0 .75 

0.70 

ICCC 

0.60 

0.8C 

C.75 

C.7C 

150C 

C.75 

C.75 

C.7C 

C.6C 

2aoc 

o.eo 

C.65 

0 .60 

0 .55 

2 50C 

0.^5 

0.5C 

0.5C 

C.35 

3C00 

0.20 

C.2C 

0 .2C 

0.20 


CURRENT STATUS OF UNIT I VARIABLE DEFINITIONS 

0 - ALIVe NOT FIPING 

1 - AL<VE AND FIRING 

2 - K a I E C 

3 - MOVING 


VEHICLE SPEED VARIABLE DEFINITIONS 


1 - 

5 

HPH 

2 - 

12 

NFH 

3 - 

15 

MFH 

4 - 

18 

MFH 
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***** THE LiNC COMEAl PHiJE BEGINS ***** 


**** DEFENSIVE FORCE IS ELl.MUATEC. END OF BATTLE. 
1IHE= 7A5 SECONDS 


^GGRESOR UNIT 
UN IT X 

INFORMATION 

Y FORCE LEVEL 

STATUS 

LGST-PCT 

TARGETS 

1 2394.5 

1978.9 0 .0 

2 

1.000 


2 31^4.2 

2246.3 4.8 

1 

0.731 

5 

3 332^.7 

iT21.6 X7.9 

0 

0.005 


DEFENSIVE UNIT 
UNIT X 

INFCPMA7I0N 

Y FORCE LEVEL 

STATUS 

LOST-PCT 

TARGETS 

^ 45CO.O 

3300.0 0 .0 

2 

1.000 


5 4fCO. C 

27CC.C G.O 

2 

1.000 

2 

6 ^cCO.O 

1800.0 0 ,0 

2 

1.000 
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